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Résumé. — On étudie de façon détaillée les processus d’absorption et d'émission de photons 
de résonance optique par un ensemble d’atomes. Le rayonnement est traité quantiquement ; 
l’ensemble des atomes, représenté par une matrice densité dans l’état excité et l’état fondamental. 
On montre que sous l'influence de l’excitation lumineuse, l’état fondamental acquiert une durée 
de vie finie et une self-énergie. Des équations différentielles sont obtenues qui décrivent l'effet 
des différentes étapes du cycle de pompage optique sur la matrice densité dars l’état fondamental. 


Abstract. — À detailed treatment is given for resonance absorption and emission radiation 


by an ensemble of atoms. 


Radiation is treated quantum-mechanically. The density matrix for- 
malism is used to describe atoms in the ground and in the excited states. 


It is shown that the 


effect of optical excitation is to give the ground state a finite lifetime and a self energy. Differen- 
tial equations are derived which describe the effect of the different steps of the optical pumping 
cycle on the density matrix in the ground state. 


Introduction. — Nous nous proposons d'étudier 
l’évolution d’un ensemble d’atomes, possédant une 
structure Zeeman dans l’état fondamental, et éclai- 
rés par une source qui émet leur raie de résonance 
optique. Les atomes passent à certains instants 
dans l’état excité en absorbant un photon, puis 
retournent à l’état fondamental en en émettant un 
autre. C’est le phénomène bien connu de la réso- 
nance optique. Il en résulte une modification des 
populations des sous-niveaux de l’état fondamental 
et de la cohérence qui existe entre eux. Le premier 
effet, dans le cas particulier d’une excitation en 
lumière o, ou 6, n’est autre que le « pompage 
optique », pouvant créer un important taux d’orien- 
tation dans l’état fondamental, déjà souvent étu- 
dié [1]. Le deuxième effet a déjà été mis en évi- 
dence, en particulier dans l’étude de l'influence de 
l'intensité lumineuse créant le pompage optique sur 
la largeur des courbes de résonance magnétique de 
l’état fondamental [2], ou dans l’étude des phéno- 
mènes transitoires qui se produisent lors de l’éta- 


(*) Un second article, faisant suite à celui-ci sera publié 
dans le prochain numéro du Journal de Physique. 
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blissement d’un champ de radiofréquence [3]. Notre 
étude généralise les résultats obtenus, en précisant 
l’évolution dans le temps, sous l’influence du pom- 
poage optique et du champ de radiofréquence 
créant la résonance magnétique, de la matrice den- 
sité qui représente l’ensemble des atomes. Les élé- 
ments diagonaux de cette matrice représentent les 
populations des sous-niveaux, les éléments non dia- 
gonaux la cohérence [4] qui existe entre eux. Nous 
obtenons certains résultats nouveaux, tels qu’une 
possibilité de conservation partielle de la cohé- 
rence lors de la résonance optique, la prévision de 
certains déplacements de la raie de résonance ma- 
gnétique sous linfluence du pompage optique, et la 
possibilité de calculer simplement de manière géné- 
rale les effets observés lors des expériences de 
« faisceau croisé » [5]. 

Le plan de cette étude est le suivant : nous pré- 
ciserons d’abord le formalisme mathématique (1) 
puis nous étudierons successivement l’excitation 
lumineuse (II) la retombée de l’état excité (III) et 
enfin le bilan global qui en résulte pour l’état fon- 
damental (IV) ; l’effet d’un champ de haute fré- 
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quence sera étudié en dernier lieu (V). Moyennant 
certaines approximations que nous préciserons, 
nous aboutissons à un système d’équations très 
voisin du système des équations de relaxation de 
Kubo-Tomita-Bloch-Ayant [6]. Nous montrerons 
que cette analogie est plus que mathématique, et 
que, malgré certaines différences, l'excitation lumi- 
neuse peut être considérée comme un processus de 
relaxation auquel sont associés plusieurs temps de 
relaxation longitudinaux et transversaux. La mé- 
thode utilisée pour traiter le champ de rayonne- 
ment et les notations sont très voisines de celles 
déjà utilisées par l’un de nous [7, 8]. Les notations 
non expliquées sont les mêmes que dans les publi- 
cations précédentes. 

Les parties (1), (IT), (111) du plan précédent sont 
traitées ci-dessous. Les autres le seront dans un 
second article. 

Les résultats de cette étude ont été publiés sépa- 
rément aux Comptes Rendus de l’Académie des 
Sciences [9]. 


I. Formalisme mathématique. — 1. ExciTrATIoN 
LUMINEUSE. — L’action de la source lumineuse sur 
un atome sera décrite par l’action d’un ensemble 
de N ondes planes, correspondant à WV photons de 
vecteurs d’onde k;, k: ... ky. Tous les k; sont 
supposés parallèles, la direction du vecteur élec- 
trique étant la même pour toutes les ondes planes, 
parallèle au vecteur unitaire ex,. Le nombre N est 
supposé très grand et on peut admettre une répar- 
tition des nombres d’ondes avec une densité u(k) dk 
dans l'intervalle (k, k + dk).u(k) n’a de valeur 
notable que dans une bande de largeur À entou- 
rant la valeur k, correspondant à la transition 
optique considérée. À, largeur de raie de la source 
excitatrice, est supposé grand devant l', largeur 
naturelle du niveau excité, w, et w;, écarts Zeeman 
dans l’état excité et dans l’état fondamental : 


AE, Ge, Of. (E54) 


Dans la réalité physique, l’atome est soumis à 
une succession de trains d’ondes et le modèle 
adopté peut en sembler éloigné. Nous n’avons 
cependant pas pu trouver de modèle plus réaliste, 
permettant de conduire les calculs jusqu’au 
bout (*). Nous discutons dans l’appendice 1 la 
validité du modèle. Nous montrons qu'il n’a de 
sens qu’en calculant des moyennes sur un grand 
nombre d’atomes ce qui correspond bien aux expé- 
riences envisagées dans cet exposé. 


2. Érars pu sysrème. — Les états du système, 
composé d’un atome et du champ de rayonnement 
que nous considérerons, sont de 3 types [7]. 

a) Les états où l’atome est dans le sous-niveau 
Zeeman de nombre quantique magnétique & de 


(*) Une des difficultés rencontrée est par exemple l’appa- 
rition d’effets d’interférence entre trains d’onde successifs. 
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l’état fondamental, en présence des V photons 
représentant l’onde incidente : 


lu ; k;, ko Or KN, ex >. 


ou plus simplement |[u >. 

b) Les états où l’atome a absorbé un photon de 
vecteur d’ondes k; dans l’onde incidente et se 
trouve dans le sous-niveau Zeeman de nombre 
quantique magnétique m de l’état excité : 


SU RNIEN Le 


Im ; k;, PAT ki—, Ki+:, 


ou plus simplement [m ; — ki >. 

c) Les états où l’atome a absorbé le photon k; et 
est retombé dans le sous-niveau Zeeman p de l’état 
fondamental en réémettant un photon de vecteur 
d'ondes k et de polarisation caractérisée par le 
vecteur e; 


I ; K;, er Ki—1 Kipi ÉCTKNR ESS k, ex > 


ou plus simplement [u ; — k; ; k, À >. 

Si l’on prend comme zéro d’énergie, l’énergie des 
N photons de l’onde incidente, les énergies de ces 
3 types d’états sont respectivement : 


Do, Ko — 3 + mo, k— kj + Lo 


(dans tout ce qui suit, on a fait % — 1, c — 1). 
Le vecteur d’état du système à l’instant £ s’écrit : 


[F() > -È ault)le > +2 am: kilt) [nm ; —Kk > 


+ D aui-k:kRA()U; —k; k AD. (I, 2) 
u,1,k,À 


L'état initial, à l’instant t — t,, est représenté 
par : 


Ft) > = : a(to)lu > (I, 3) 


où seuls les a,(£,) sont non nuls. Le but du forma- 
lisme développé ci-dessous est de calculer les a, (t) et 
Qui —k;: KA(E) qui décrivent l’état de l'atome avant 
et après le phénomène de résonance optique, ou plus 
précisément les éléments de la matrice densité 
au(t) (et) et 2, us ki: RAGE) QG: = k,: KA(E), calculés 
en faisant la moyenne sur les V atomes identiques 
d’une vapeur. Les éléments diagonaux (4 = y’) de 
cette matrice représentent les populations des ni- 
veaux pu, les éléments non diagonaux (4 pu’) la 
«cohérence » entre les niveaux. Nous supposerons 
qu’une telle cohérence existe à l’état initial 


(au (to) 4 (to) À 0) 
et nous étudierons son évolution au cours du temps 
sous linfluence du pompage optique. 


3. Équarions Dp’ÉvoLurTIon. — L’Hamiltonien 
du système comprend deux termes : 
Je = Jo + Jr. (1, 4) 


N° 6 


ais Sens 6e ét EEE 


Sr ot dtapton fn de duos utile items és. nie 


bu. dt - 


N°6 


J, représente l’énergie de l’atome et du champ 
de rayonnement. IT est diagonal et ses valeurs 
propres sont les énergies des états données ci- 
dessus. Jr représente le couplage entre le rayon- 
nement et l’atome. Ses seuls éléments de matrice 
non nuls sont donnés par (cf. [7], éq. (I, 13) et 
(I, 14)) 


< ufi@rim ; —k; > = Ame kiR < ule,, .D|m > 
Cu; —k;; k, ANrIm ; —k > 
= Age kR < ple).D\m > 


(L, 5) 


D est (à un facteur près) le moment dipolaire de 
l'atome. A4 dépend des fonctions d’onde de l’atome 


e 1 : x 
et est proportionnel à he R détermine la posi- 


tion de l’atome, qui dépend de sa vitesse v : 


R = Ro —- vt 
| (0) > satisfait à : 


LP > = (fe + PU) >. (6) 


Passons en représentation d'interaction. On pose: 
[D(4) > = eid@et [P(r) > 


rÈ blt)lu > +2 Om ; =ki(t)|m ; — ki > 


Fe À dei ki ka (tilu : —k:k, x > 
Ho 
Kilt) — ne, ) e— ot. (12) 
Il vient : 
ID 7 = KP) 7 (I, 8) 
Les coefficients b satisfont à : 
bu =D < ulEr(t)im ; — ki > bmi k; (T, 9a) 
mi 
ibm: k; = 2 < mi — kite > by 
nm 
+ EE < melilu ; KA bui-k:ka (I, 9) 
u’,k,A 
bu: —ki:kA=D Cu; k, AIZ{(t)|m > bm:—k; (I, 9e) 
\ m 
avec 
D u|3Cr(t)|m ; — ki > — Aye KR. 
< ule10.DIm > eëlki— ko + boy — mo)t (I, 104) 
Ou; k, AI (t)Im = 
Run ki ki AE ()Ims — k; > = Age tk.R 
< p'le.Dim > el ko Fu — mot (1 108) 


(I, 9a) décrit l'absorption d’un photon et l’exei- 
tation de l’atome ; (1, 9b) décrit l’évolution de 
V’état excité ; (1, 9c) décrit la réémission d’un pho- 
ton et le retour dans l’état fondamental. 
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IT. Étude de l’excitation optique. — A. Évoru- 
TION DE L'ÉTAT FONDAMENTAL. — Le principe du 
calcul consiste à éliminer by/:_K,:Kkx entre (1, Qc) 


et (1, 96), puis bn: -Kk, entre (T, 9b) et (TI, 9a). 


1. Elimination de bs:_K,: ka. — On intègre (I, 9c) 
et on porte dans (1, 9b). [1 vient, avec la condition 
initiale du’: —ky: k(Lo) — 0 
10m ; — k; —= 5; cs ML, 


; — kildr(t lu > bu 


. Fe : , 
—i E, J <b> 


<OWRkRAIES (1) 1m" > bm’: —k;(t) de. (1,1) 


Le deuxième terme représente les effets de l’émis- 
sion spontanée sur l’évolution de l’état excité. Le 
calcul de ce terme est identique à celui que l’on 
rencontre dans l’étude de l’émission de rayon- 
nement par un atome isolé. On trouve que (IE, 1) 
est équivalent à [7]: 

im: —k, = À TE 
[3 


— kr (t)lu 7 by 


= 1e 1 iAE) bn _k(t) (IL, 2) 


F'et AE sont la largeur naturelle et la self énergie 
du niveau excité. Nous poserons : 


PAPE RON 


2. Elimination de bn: -k, — On intègre de 
même (II, 2) et on porte dans (I, 9a) avec la condi- 
tion initiale + —k;(to) = 0 


bu = — >» 


mi,’ su 


Es m ; 


E 
ta HiAE)4+ < Jr (t)[m ; — ki > 


— kil (r)lu" © but) dé. (II:3) 


Utilisant (1, 10a) et remplaçant la sommation 
sur & par une intégration sur À pondérée par la 
forme de la raie optique excitatrice u(k) (cf. $ I, 1) 
on obtient 


ZX: 


< ulex, .Dim > < mJlex,.Dlu’ > eitu—uot 


-+ co LOTS LE , : 
AE de Rate nes 
— CO b 


bul#) dt” (II,4) 


L’exponentelle qui figure sous l’intégrale est une. 
fonction de X très rapidement oscillante, et dès 
que  —t" > 1 JA où À est la largeur de u(k), l’inté- 
grale sur # du produit de cette fonction par 

k)|Axl? qui varie très lentement, est nulle. Les 
seules valeurs de {’ qui contribuent à un résultat 
non nul sont telles que t{ > # > t — (1/JA). 

Si b,(t) varie peu sur un intervalle de temps de 
l’ordre de 1/A (ce que nous justifierons par la 
suite, Cf. $ IV, 3), on peut remplacer b,.(#) par 
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b(t). L'intégrale sur t’ se calcule alors aisément et 
il reste à calculer l’intégrale sur 4 : 


so uk) [Ael? dk 


1 el —h—mo +0, —k y +iT/2 (te) 


. 


— ik — PASS moe + mo — k.v + iT/2) 


Pour t — 1, > 1/1, seul cas intéressant puisque 
la vitesse de variation de b,(t) est lente devant 
(cf. $ IV, 3) on peut négliger le terme exponentiel 
et 1l reste pratiquement 


ET ARE de —@ — © (5 

—co RÉ ko) + il f2 

et 

D 2) CAPE A ee LT (IT, 6) 
u” à 

en posant 


Ayy = D < ule»,.Dim > < mle .Dlu’ > (11,7) 
m 


EME = E 


dk 


FE ER (ko + kov)]e + (2/4) 
ITS, 8, a) 
AE = [ atylAal? LE La dk. 
me LR — (ko + ko.v)]? + (72/4) 
(IL, 8b) 


Dans ces expressions, nous avons négligé mu, et 
u'ey, qui sont très petits devant la largeur A 
de u(k) qui détermine la vitesse de variation de 
u(k)|A:/2. Nous avons remplacé de même k,.v 
par k,v (effet Doppler qui n’est pas très petit 
devant À, mais devant 4,, n’introduisant donc 
ainsi qu’une erreur négiigeable.) 


3. Interprétation physique de |”. — La partie 
réelle de (IT, 5) est une courbe de Lorentz, de lar- 
geur l' centrée au point 4, + k,.v, dont l’aire est 
égale à 7x. [étant très petit devant À, l’intégrale 
(II, 8a) est peu différente de 


'f2 [7 ut) VAuf® môUk — (fo + ko.v)] 


— ru(ko + kov) |Aÿ|? (IL, 9) 
F” s’interprète comme proportionnel à la proba- 
bilité par unité de temps que l’atome quitte l’état 
fondamental en absorbant un photon. Cette proba- 
bilité, calculée au moyen de la théorie des pertur- 
bations, est en effet [”4,, pour le niveau u. 
1/1"A,, s’interprète donc comme une durée de 
vie 7, du sous-niveau 4. 7, dépend de u par 4, 
c’est-à-dire par l’état de polarisation de la lumière, 
et est inversement proportionnelle à l’intensité 


N°6 # 
lumineuse excitatrice à la fréquence d’absorption 


déplacée par effet Doppler ko + Kko.v. 
Nous poserons I” — 1/T; (II, 10), pour nous 
conformer à des notations déjà utilisées [3]. 


4. Interprétation physique de AE". — La partie 
imaginaire de (II, 5) est une courbe de dispersion 


de largeur l', centrée au point 4, + k,.v. L’inté- 
grale (IT, 8b) est nulle si le centre k, de la raie 


excitatrice est confondu avec ko + k,.v. En pre- 
nant pour u(k) une forme simple telle qu’un créneau 
ou une courbe de Lorentz, on peut préciser l’allure 


générale de variation de AE" avec k; — (£, + ko. v) : 
AËE' commence à croître avec k; — (%, LKR 


Un maximum est atteint pour #; — G, + ko.v) de 
l’ordre de A 7/2 où A est la largeur de u(k). L'ordre 
de grandeur de ce maximum est 16/2 où I est 


la valeur de [" lorsque k, = k) + kov. ; AE” dé- 
A 


croît ensuite et varie comme — 


£ : ki —(k + ko.v) 

pour À « ka — (ho + ko.v). 
L'interprétation physique de A£’ est la sui- 
vante : Considérons un atome dans l’état fonda- 
mental en présence d’un photon dont l’énergie k#, ne 


coïncide pas avec l’énergie résonnante k, + k,.v. 
L’absorption réelle du photon par l’atome ne peut 
avoir lieu : en effet l’énergie ne serait pas conservée. 
Mais cette transition peut se faire de façon « vir- 


tuelle » à condition de durer un temps inférieur 
1 


RE AN UE de 

he CREER) 
certitude. L'effet de ces transitions virtuelles sur 
l’atome est de ramener dans l’état fondamental 
une partie de l’énergie du niveau excité, donc de 
déplacer le niveau fondamental. AZ’ est donc la 
self-énergie de l’état fondamental en présence du 
rayonnement excitateur. (Remarquons également 
que par suite de ces transitions virtuelles, les pho- 
tons se propagent moins vite au sein de la vapeur 
atomique ; le déplacement qui apparaît dans nos 
calculs est donc étroitement lié à la dispersion de 
cette vapeur atomique.) 

Notre formalisme, tenant compte uniquement du 
niveau de résonance, ne peut décrire l’effet des 
transitions virtuelles vers tous les niveaux excités. 
Nous montrerons cependant plus loin (chap. IV) 
que si nous ne pouvons atteindre la self-énergie 
absolue de l’état fondamental, nous pouvons appré- 
cier la différence des selfs-énergies des sous- 
niveaux . C'est cette différence qui est intéres- 
sante parce qu’elle se manifeste sous forme d’un 
déplacement de la raie de résonance magnétique. 
Signalons qu’un déplacement de ce genre vient 
d’être observé expérimentalement par l’un d’entre 
nous [10]. 


en vertu de la 4€ relation d’in- 


. N°6 


5. Introduction de la matrice densité. — T’équa- 
tion (II, 6) permet d’écrire : 


1 ; LAS 

bu = — (57, + AP) E Anne ee (I, 14) 
1 Se 

Posons 


ouu’ = by br (IT, 12) 
Puu eSt un élément de la matrice densité repré- 
sentant l’état de l’atome dans l’état fondamental. 


A partir de (II, 11) on peut calculer la vitesse de 
{1 


variation de la matrice densité “Ar Pau” due au pro- 


cessus d’excitation. 


du) 1 re D nr 


(Eh ar) 3 Au ed qe. 
2Ty u” 


(IL, 13) 


Les équations (II, 13) sont strictement équiva- 
lentes aux équations (II, 11). Effectuons mainte- 
nant la moyenne des 2 membres de l’équa- 
tion (II, 13) par rapport à l’ensemble des .N atomes 
de la vapeur. Nous voyons apparaître au 2° membre 


des termes tels que ee + AE") Eu. L'effet 
2T, 


Doppler varie d’un atome à l’autre. Il en est donc 
de même pour 1/7, et AE". Mais, par suite des 
chocs, la valeur de 0, pour un atome déterminé 
et sa vitesse n’ont pas de corrélation précise de 
sorte que l’on peut calculer la moyenne indépen- 
damment pour les deux facteurs (1/27T,) + 1AE' 
et pw du produit précédent. Nous redéfinissons 
donc 1/27, et AE’ comme étant les moyennes des 
expressions (II, 8a) et (IT, 8b) par rapport à la 
fonction de distribution de k,.v qui n’est autre 
que la courbe d'absorption Doppler. 

Dans ces conditions l'équation (II, 13) valable 
pour chaque atome, est valable pour la matrice 
densité qui représente un ensemble de .N atomes. 
Elle nous permet de décrire l’évolution des popu- 
lations des niveaux (,,) et de la cohérence entre 
eux (0,4) Sous l’influence de l’excitation lumineuse. 


6. Etude de cas particuliers. — a) Excitation en 
lumière o+ Ou o— ou ñ. 

Dans ce cas, les éléments de matrice 

< ue, . Dim > 

ne diffèrent de 0 que si m—u — + 1, — 1 ou 0 
respectivement. Dans ce cas A, n’est différent 
de 0 que pour u — un’. Une teile excitation est 
dite non cohérente. C’est le cas usuel des expé- 
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riences de pompage optique. Dans ce cas (IT, 13) 
s'écrit 
dat 1 
pra | Met ee L (Auu + Au’y’) 

+ AE" (Au — Auvt)] Euu’ (LT, 14) 


Cette équation met bien en évidence, la cons- 
tante de temps T\,,, telle que 
1 il 
Ta La 27. Ce NL Au): 
avec laquelle la cohérence disparait, la variation de 
la distance énergétique A£'(A,, — A) entre les 
deux niveaux — et la «durée de vie» 7, telle que 
LIT, = AIT, =" 4Ajs/To du niveau y. 

b) Excitation « cohérente ». — C’est le cas général 
où Ayy 7 0 même pour u Æu’. C’est en parti- 
culier le cas des expériences de «faisceau croisé »[5]. 
Notre formalisme permet aisément de calculer 
l'absorption d’un faisceau lumineux. En effet, le 
nombre de photons absorbés par unités de temps 
est égal au nombre d’atomes qui quittent l’état 
fondamental par unité de temps. L’absorption est 
donc donnée par — X (dt ,.. /dt). 

[TA 


B. ÉVOLUTION DE L'ÉTAT EXGITÉ. — L’intégra- 
ton de (IT, 2) donne : 


” NO ff? —(dT/2+iAEN—+ 
ibm k4 = À le (T/2+iABt—t 
[0 0 


< m; — kiR£(r}lu > bult) dé. (II. 15) 


Prenant le complexe conjugué, l’indice m étant 
remplacé par m',et multipliant membre à membre 
l’équation ainsi obtenue par (IT, 2) on trouve 


; iR : 
bre Ke Eee ha (; 1 ia) Lea (sr exp 
î î 2 î i 
= S eil(m—mo,—{(u—u1o;lt mlex, .Dlu = 


uu” 
< u'lex,.D|m’ > |Ax,[? 


t ñ 7 L [114 3 9 LA 
à Con AE amat dr ge CAT TON 
to 
(IL, 16) 
Nous nous intéressons à la cohérence dans l’état 
excité et aux populations des sous-niveaux M, 


mais non au photon k; absorbé. Nous étudions 
donc la quantité 


s' K Q » 
nn ete Joux dt x uUR) dx 


| à 
dont la dérivée qi Pr" est donnée par (II, 16) et 


son complexe conjugué. On rencontre alors les 
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mêmes intégrales qu'au $ II À 2 et on trouve 
(cf. réf. 8) 


: Aer Les = 
D Om: k, Pi k, — — (5 D AE) 2 bnik, di, 


LA 


G ss ia") Qil—-mo,—(u—uo;1t 


< mjex, .Dlu > 2 pau ijex Dim (II, 17) 
et 
US 
dt To uu’ 
rex :Dlu > pu Siufex Dim >. : {IL 148) 


Le 1% terme du 2° membre de (II, 18) repré- 
sente l’émission spontanée, le 2° terme l’effet du 
pompage à partir de l’état fondamental. 


HI. Étude de la réémission de rayonnement. — 
1. MATRICE DENSITÉ APRÈS RÉÉMISSION DE RAYON- 
NEMENT EXPRIMÉE EN FONCTION DE L'ÉTAT EXCITÉ. 
— L'intégration de l'équation (I, 9c) donne : 


| re 
bu: — ki: kAÏt) — E | < uk, AE (21m > bmi— k(t') dt 
m 0 


(IH, 4) 


Comme au $ II-B, nous prenons l’équation com- 
plexe conjuguée, avec l’indice 4’ remplaçant u, et 
multiplions membre à membre cette équation et 


(I, 9, c) : 


Z emo tu-mojt |Ax|à 
mm 


Du: k,: Kk,X Bi k: RNA 


< pile] Dim > < m/le1.Dlu’ > 


DRE ’ 
L eTro—k+k.y—uo,+mo l(t—# A0) bete) dt’. 
(IL, 2) 


Nous nous intéressons là encore à la cohérence 
globale qui existe dans l’état fondamental de 
atom: après réémission, ainsi qu'aux populations 
des niveaux |u =>, et non aux photons absorbés ou 
réémis. Nous devons donc sommer les 2 membres 
de (III, 2) par rapport à k;, k, À La sommation 
sur & revient en fait à supprimer l’indic? — k; dans 
les deux membres. La sommation sur (k, À) com- 
porte une sommation sur [k| puis sur la direction 
de k et la polarisation à. La sommation sur |k| 
conduit à évaluer 


>> |Ax/? 


14 fs 


ho —+ ke vue, + mo (1) 


| UE DER 
Z 5% k (0 bus ke dt = (5 ++ AE) pl) 


O I 


(III: 3) 
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le calcul étant déjà fait dans la théorie de l'émission 
spontanée [7] ; on a donc : 


Rnb tee 
È ui ki: KA DE ki: KA 


8x \2 
< wle1.D\m > < mlex.Dlu > Emm. 


= #ÿ G iAE) D, etlm—mo,—(u—He,lt 
um 
(III, 4) 


Cette équation permet de calculer la quantité de 
lumière de fluorescence émise par unité de temps 
dans une direction déterminée avec un état de 
polarisation donné, qui n’est autre que 


2 Chuikska bg, R; ka] (Cf: réf. 8) 

Pour achever la sommation sur k, il convient 
d’expliciter la nature des éléments de matrice 
< ule,.DIm >. Nous nous limiterons au cas d’un 
état fondamental du moment angulaire électro- 
nique J — 0 avec spin nucléaire Z, avec comme 
état excité le sous-niveau hyperfin # d’un niveau 
de moment angulaire électronique J = 1 (cas des 
isotopes impairs du mercure, pour la raie 2 537 À, 
615, —6%P;,). Avec les notations de la réfé- 
rence 7, On à : 


< ue Dim > = CirlF,m; m—u,w 


< pu = Ole1. Dlmz = m—uw>r-06 (II, 5) 


CutF,m;m —4u,u)est un coefficient de Clebsch- 


Gordan. L2 principe du calcul est évidemment le 


même pour le cas des atomes alcalins pour lesquels 
la transposition des résultats est aisée. 
La sommation sur les angles donne alors [7] : 


8T 


— Om, m'y 
] 

et 1l reste : 

2 ARE k;: k:2 CHE k;:k, à 


ik 
T : 
Æ Ë QE (AE) » 
m,m 
m—m =u—u’ 


etu—u)lo —e,) 


# 


v 


CuiF,m; m—u,u) Cf, m°; m°— 4", 4°) Emme’ 


et finalement, la vitesse de variation de la matrice 


densité de l’état fondamental d® 4, /dt due au 
processus de retombée est donnée par : 


. = d 
dEpun [dt = 2 (RER À ri kiik:X) 


Er 


mm 
m—m'=u—u 


eitu—u)(o,—0,)1 


Cu(F,m; m—4u,u) Cir(F, n°; m — 4,4) pmm. 


On voit qu’on ne peut avoir d® o,, {dt Æ 0 
pour & =£ L' Que Si pmm 7 0 au moins pour certains 
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couples (m, m') tels que m -< m'. En d’autres 
termes, il ne peut y avoir de cohérence dans l’état 
fondamental après réémission de rayonnement que 
s’il y en a dans l’état excité avant l’émission. 


2. MATRICE DENSITÉ APRÈS RÉÉMISSION DE 
RAYONNEMENT EXPRIMÉE EN FONCTION DE L'ÉTAT 
FONDAMENTAL. — Nous allons chercher à exprimer 
la vitesse de variation de P,, due à la réémission 
d’un photon en fonction des valeurs de pb, avant 
Pabsorption en éliminant pmm entre (11, 18) et 
(III, 7). L'intégration de (II, 18) donne 


e—Trt f: t , NT ee / 
Enn' = freres ue, —iTit 
To un Jt 
< mlex, .Dlu > puu’(t) < u’ler, Dim’ > dé. (III, 8) 


Si l’on reporte (IIT, 8) dans (III, 7) on obtient 


dû) Quu’ se 1% eu (u—u/"lo;t 
dt Th up 
mm, < mlex,.Dlu” > < u"lex, .D|m’ > 
MN—M = 


CrfF,m; m—u,u CulE,m'; 


m'—u", u) 


t î (ES F2 2 7 17 / A 
Î ele, + (ue, HT) ont) dt. (III, 9) 
to 


L'interprétation physique de ce résultat est la 
suivante : la retombée dans l’état fondamental à 
l'instant t{ dépend de toute |’ « histoire » de l’état 
fondamental entre t, et t. Mais le facteur eFé#) 
fait-que les valeurs de pur, (t') pour des instants #’, 
très antérieurs à { — (1/[') contribuent très peu. 
En effet, les atomes qui réémettent et retombent 
dans l’état fondamental à l’instant t l’ont quitté à 
un instant {', et t — {ne dépasse guère 1/[", durée 


_de vie de l’état excité. 


Si nous admettons que la vitesse de variation 
de p,, est lente devant 1/1 (ce que nous justi- 
fions plus loin, cf. $ IV, 6) on peut donc remplacer 
dans (III, 9) pyry(t) par pymu(t) et l'intégrale 
sur {’, pour des instants f tels que t — 1, > 1/F, 
seul cas intéressant en pratique,"se réduit à 


Puu”(1) 


Do, du 4) @ 


Donc, la variation de la matrice densité due à la 
retombée de l’état excité est donnée par : 


1 0) — (nu) ’ 
dt?) Ouu’ si : 4Lu—0)— (00 Not bth, 


LS Qu’ e 1 
dt DT or we] 
(III, 10) 
avec 
Bo — Die < mlex,.Dlu” > < u"'er, .D|m’ > 
ones 
CulF,m;,m—quu) Cu(Fm;m vw) (I, 11) 
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APPENDICE I 


Justification du modèle décrivant l’excitation 
umineuse. — Un état tel que 
HÉRNEN AREN SE (1) 
est un état stationnaire du champ de rayonnement 
lorsqu'il n’interagit pas avec les atomes. Chaque 
photon de vecteur d’onde ki, est associé à une 
onde plane indéfinie et, ce qui signifie physi- 
quement que la probabilité pour que le photon se 
manifeste dans une interaction en un point quel- 
conque de l’espace est la même dans tout l’espace 
et est constante dans le temps. 
Un état tel que 


M DR (RTC NÉE REC) 
Lake 

peut, par un choix convenable des «, représenter 
un paquet d’ondes qui se propage. Il représente 
donc de manière beaucoup plus réaliste la situa- 
tion physique : une source émet des trains d’ondes 
qui passent l’un après l’autre au point où se trouve 
un atome. La matrice densité représentant |[Ÿ > 
par rapport aux états de base du type (1) du champ 
de rayonnement n’est pas diagonale. [F > n’est 
en général pas un état propre de l’'Hamiltonien du 
champ de rayonnement en l’absence d’interaction. 
Il est beaucoup plus difficile de mener les calculs en 
prenant comme état initial un état du type (2) qu’un 
état du type (1). 

Une difficulté analogue se présente en théorie de 
la relaxation spin-réseau dans les solides pour la 
représentation des ondes de phonons dues à l’agita- 
tion thermique (réseau). On est amené à représenter 
le réseau par un mélange statistique d'états, en 
équilibre thermique, donc par une matrice densité 
diagonale, et on suppose que l’énergie du réseau 
est tellement grande devant celle des spins que les 
échanges d’énergie entre les spins et le réseau ne 
modifient pas la matrice densité de celui-ci. La 
première hypothèse se justifie par le fait que nous 
ne disposons d’aucune donnée précise sur le réseau, 
et que les effets liés à la présence d'éléments de 
matrice non diagonaux par rapport au réseau ont 
une moyenne nulle quand on étudie l’action du 
réseau sur un très grand nombre de spins, c’est-à- 
dire si on répète un très grand nombre de fois une 
même expérience. 

Guidés par cette analogie, nous représenterons 
l'excitation lumineuse par une matrice densité dia- 
gonale. Comme nous étudions toujours en pratique 
l’action du rayonnement sur un très grand nombre 
d’atomes, cela revient à faire en même temps un 
très grand nombre d’expériences, et l’on peut alors 
négliger les effets des termes non diagonaux, dont 
la moyenne est nulle parce que leur phase varie 
aléatoirement d’une expérience à une autre. Or, 
nos équations sont linéaires, par rapport à la mas 
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trice densité, et l’on peut admettre que tous les 
états présents dans le mélange statistique initial 
correspondent aux mêmes valeurs. des para- 
mètres T,, AE", etc. (c’est-à-dire au même flux 
lumineux excitateur). Il est alors suffisant de faire 
le calcul, comme nous l’avons fait, pour un seul état 
du type (1). On admet en outre que l’absorption du 


N°6 


rayonnement par les atomes ne modifie prati- 
quement pas ses caractéristiques, c’est-à-dire que la « 
matrice densité représentant l’état de l’ensemble du 
système en ce qui concerne le rayonnement (la 
trace étant prise sur les états de l’atome) reste 
toujours la même. 


Manuscrit reçu le 16 janvier 1961. 
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ÉTUDE DE L'ANTIFERROMAGNÉTISME HELICOIDAL DE MnAu, PAR DIFFRACTION DE NEUTRONS 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


Par ANDRÉ HERPIN et Prerre MERIEL, 


Service de Physique du Solide et de Résonance Magnétique, Centre d'Études Nucléaires de Saclay. 


Résumé. — Une élude par diffraction de neutrons a montré que la structure de MnAu, quadra- 
tique centrée, peut être considérée comme un empilement, normalement à l’axe c, de feuillets de 
manganèse séparés par deux feuillets d’or. MnAu, est d’autre part un hélimagnétique : tous les 
atomes de manganèse d’un même feuillet ont leurs moments parallèles entre eux et normaux à 
l'axe ec, ces moments tournant d’un angle constant lorsqu'on passe d’un feuillet au suivant. Un 
champ parallèle à ec influe peu sur cette structure magnétique. Par contre un champ normal à € 
modifie brusquement cette structure dès qu’il dépasse une valeur seuil de 9,6 kOe. Une théorie où 
l’on tient compte des interactions d’échange entre atomes de feuillets premiers et seconds voisins 
Re de rendre compte quantitativement de toutes les propriétés de cet alliage qui ont été 
observées. 


Abstract. — A study by neutron diffraction has shown that the tetragonal centred structure 


* of MnAu, can be considered as a packing normally to the c axis of layers of manganese separated 


by two layers of gold. Moreover MnAu, is a helimagnetic : all atoms of manganese of the 
same layer have their moments parallel between themselves and normal to the c axis, these moments 
turning by a constant angle when passing from one layer to the next one. A field parallel to c 
has little influence on this magnetic structure. On the other hand a field normal to c modifies 
abruptly this structure as soon as it goes beyond a threshold of 9.6 kOe. A theory, where one 
takes into account exchange interactions between atoms of first and second neighbouring layers, 
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qantitatively accounts for all the properties of this alloy which have been observed. 


I. Introduction. — Les alliages or-manganèse 
ont fait l’objet de nombreuses études. En 1953, 
Raub, Zwicker et Baur [1] donnent le diagramme 
du système Au-Mn et mettent en évidence le com- 
posé défini MnAu,. En 1956 Meyer et Taglang [2] 
font une étude magnétique très complète de cet 
alliage à l’état ordonné. Pour une température infé- 
rieure à 90 0C, ce composé se comporte comme un 
paramagnétique fort dans les champs faibles, puis à 
partir d’un champ de l’ordre de 10 kOe, son aiman- 
tation croît rapidement et tend, difficilement d’ail- 


M uern/g 
3 


leurs, vers une saturation au-delà de 20 kOe (fig. 1). 
Ces propriétés sont analogues à celles des composés 
métamagnétiques étudiés par Néel [3]. Néanmoins, 


ni la structure cristalline du composé ordonné, ni sa 
structure magnétique n’étaient connues. Une étude 
par diffraction de neutrons a été entreprise afin de 
déterminer l’une et l’autre. 


1) DisPosiTIF EXPÉRIMENTAL. — L'or ayant une 
forte section efficace d'absorption (60 barns pour 
les neutrons de 1,143 À de longueur d’onde utilisés), 
l’échantillon est une plaquette polycristalline 
d’épaisseur optimum placée en position de trans- 
mission symétrique sur l’axe du spectromètre. Cette 
plaquette a été préparée par A. J. P. Meyer (1) qui 
a vérifié que ses propriétés magnétiques étaient 
bien identiques à celles de l’échantillon ayant fait 
l’objet de son étude [2]. La plaquette peut être 
entourée d’un four permettant de la porter au- 
dessus du point de Néel. Un électroaimant permet 
d'appliquer à l’échantillon un champ magnétique 
dont la direction est soit parallèle au vecteur de 
diffusion, soit perpendiculaire au plan de diffusion. 

Le spectre de diffraction obtenu à la température 
ordinaire est représenté sur la figure 24. Un spectre 
réalisé à une température supérieure à la tempéra- 
ture de Néel permet, par différence (?), d’obtenir le 


{4} Nous remercions bien vivement A. J. P. Meyer de 
l’aide précieuse qu’il nous a apportée. 

(2) Pour obtenir le spectre « magnétique » par différence, 
il est nécessaire de tenir compte de la dilatation thermique 
qui déplace les raies vers les petits angles quand la tempé- 
rature croît. Nous avons pu déduire de la correction ainsi 
faite le coefficient moyen de dilatation & = (25 + 4).10—5, 
mais la précision n’est pas suffisante pour permettre de 
déterminer séparément œa et es 
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spectre de diffraction « magnétique » (fig. 2c). 
L'examen de ces deux spectres montre immédia- 
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(a) Spectre à 293°K 
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tement que les raies du spectre à température ordi- 
naire sont soit purement nucléaires, soit purement 
magnétiques. 


2) STRUCTURE CRISTALLINE DE MnAu,. — Les 
raies nucléaires sont bien expliquées en adoptant la 
maille quadratique centrée proposée par Smith et 
Street [4] pour laquelle 


a = 3,369 À et cfa — 2,599. 


Cette maille est à peu près triple de celle anté- 
rieurement proposée pour le composé désordonné 
par Raub, Zwicker et Baur [1]. En accord avec la 
densité, chaque maille comprend deux molécules 
MnAu,. La seule structure compatible avec les 
valeurs des rayons atomiques de Mn et Au et avec 
les intensités observées correspond aux positions 
suivantes des atomes dans la maille : 


Mn :0,0,0 


Æ!: Ë SAT l 
ATIE 0, 0, Z 0, 0, Z 9° 92 + Z 2? 

Dans une première étude [5] nous avions trouvé z 
pratiquement égal à 1/3. P. Michel [6] a simulta- 
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nément déterminé cette structure par diffraction 
électronique. La valeur du paramètre z — 0,338 


qu’il obtient est certainement plus précise que celle 


que nous avions proposée, son étude portant sur un 


nombre beaucoup plus grand de réflexions. Nous. 
avons par la suite utilisé pour le calcul des inten- … 


sités la valeur ci-dessus du paramètre z. 
La structure ainsi obtenue est représentée sur la 
figure 3. On voit que les atomes de manganèse 


O Mn 


(DAu 


F16. 3. — Structure de MnAu.. 


forment des feuillets perpendiculaires à l’axe € 
séparés par des feuillets d’or. Dans le tableau I on 


TABLEAU I 


ENTOURAGE D’UN ATOME DE MANGANÈSE DANS MnAu, 


Ç 


ORDRE DISTANCE COORDONNÉES NOMBRE 
DE VOISINAGE Mn-Mn En À RELATIVES DE VOISINS . 
1 So a, 0, 0 k 
2 4,76 aa, 0 4 
$ 4,99 a/2, a/2, c/2 8 
A 6,74 2a, 0,0 A 
5 6,90 3a/2-a/2,c/2 16 
6 FE0S DRE al 8 
7 8,38 EVANS UPOE 8 
8 8,76 0,0,c 2 


trouvera les données relatives au voisinage entre 
atomes de manganèse. 


II. Structure magnétique de MnAu, en champ 
nul. — La détermination de la structure magné- 
tique s’est avérée beaucoup plus délicate. Nous 
avons vérifié l’origine magnétique des raies du 
spectre différence en étudiant la variation de leur 
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intensité en fonction de la température. La figure 4 
montre cette variation pour la première de ces raies. 
L’allure générale est correcte, l’intensité devant 
être proportionnelle au carré de l’aimantation de 


Intensité 


75 


FIG. ‘4. 


chaque sous-réseau. On peut cependant remarquer 
que l'intensité de cette raie ne s’annule pas au 
point de Néel, ce qui traduit la persistance d’un 
ordre à courte distance. 

Si l’origine magnétique des raies figurant sur le 
spectre différence est certaine, 1l n’en reste pas 
moins difficile de les indexer dans une maille mul- 
tiple simple de la maille cristalline. Même si l’on 
parvient à trouver une telle maille, encore faut-il 
expliquer les très nombreuses extinctions observées. 
Nous avons pensé qu’un tel spectre pouvait tra- 
duire l’existence d’une structure hélicoïdale, telles 
que celles décrites par Yoshimori [7] et Villain [8]. 


1) ANTIFERROMAGÉTISME HÉLICOIDAL. — Nous 
ferons seulement une application de la théorie géné- 


rale de l’antiferromagnétisme hélicoïdal que l’on 


peut trouver en [7]et [8] au cas particulier qui nous 
intéresse. Comme nous l’avons signalé, on peut 
considérer que les atomes de manganèse sont dis- 
posés en feuillets équidistants, normaux à l’axe 
quaternaire, distants de c/2. Nous supposerons que 
les interactions d'échange entre atomes d’un même 
feuillet sont fortement ferromagnétiques, de sorte 
qu’en dessous de la température de Néel, tous ces 
atomes ont leurs moments parallèles. Nous pour- 
rons donc traduire l’interaction entre un atome 
appartenant à un feuillet z et l’ensemble des atomes 
du feuillet à + q par le champ moléculaire 


H;,; — nq Mi+a 


où M; représente le moment magnétique d’un 
atome appartenant au feuillet 7. En se limitant aux 
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interactions d’un atome avec les quatre feuillets 
voisins, le champ moléculaire agissant sur l’atome 
prend la forme : 


a= +2 à 
Hi =0 3; His = 70 Mi 
1=—2 


+ (Mit + Mis) + n(Mite + Mio). 


Les coefficients de champ moléculaire ainsi intro- 
duits sont des moyennes pondérées des intégrales 
d’échange entre les différents voisins. En désignant 
par J, l'intégrale d’échange entre pièmes voisins, 
on à : 


| 
RS ES MIE CU 4 
Au 
q 
Rire + 8J; + Le + Me) ; (1) 
AUS 
ns sine + ) 
| 2 Te 8 ss.) 


L’énergie correspondante, rapportée à un atome, 
est alors : 


W = — 7 M8 — n(Mi+, + Mi) M; 
— n(Mi+s + Mio) Mi. 


Pour qu’il y ait antiferromagnétisme le champ 
moléculaire doit être parallèle au moment magné- 
tique Mi et être de module constant. On en conclut 
que, sauf cas exceptionnels, M;,,, M;_, et M; 
doivent être parallèles à un plan P, donc que tous 
les moments magnétiques doivent être parallèles 
à P. De plus, l’angle 6 de M;_, avec M; est égal à 
l’angle de M; avec M;,,. Remarquons que ces 
conclusions ne sont valables que si les différents 
atomes sont cristallographiquement identiques, 
donc sont aux nœuds d’un réseau de Bravais, et 
dans ce cas, le parallélisme des moments à un 
plan P n’est pas une hypothèse. L'énergie peut 
alors s’écrire : 


r né ’ 
W = — us — ni bo COS  — 2n9 Lo COS 2 p. 


u, étant le moment magnétique d’un atome de 
manganèse. 

En minimisant cette énergie par rapport à o on 
trouve trois solutions possibles : 

a)o — 0 : tous les moments sont parallèles. C’est 
le cas ferromagnétique. 

b) © — tx: deux feuillets successifs ont leurs 
moments antiparallèles. C’est le cas antiferro- 
magnétique simple. 

c) + donné par cos o — — n;fäns. Les moments 
magnétiques sont disposés en hélice. La résultante 
des moments est encore nulle. On a alors affaire à 
un antiferromagnétisme d’un type particulier que 
pour abréger nous appellerons hélimagnétisme. 

En fonction des coefficients n9, n3; 9, l'énergie a 
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pour expression dans ces trois cas, en faisant abs- 
La 


n 
traction du terme constant — © pu, : 


2 
Wy = — 2o(ni + n$) 
Wa = — 24o(— mi de ne) 
n°3 , 
Wi = — 2 He 1). 


La figure 5 représente la variation de l'énergie en 
fonction du rapport æ — n,/n, dans le cas où 2 


RE) 
Fire. »: 


est négatif (interactions antiferromagnétiques entre 
feuillets seconds voisins). On voit que si la solution 
hélimagnétique existe, ce qui n’est possible que 
si [x] << 4, c’est la solution la plus stable. Pour 
2.4, la solution stable est antiferromagnétique 
simple, et ferromagnétique si æ < — 4. Quand 
est positif, les seules solutions possibles sont ferro 
ou antiferromagnétiques, résultat qu'exprime la 
figure 6. 


2) DIFFRACTION DES NEUTRONS PAR UN HÉLI- 
MAGNÉTIQUE. — Les formules usuelles donnant le 
facteur de structure des raies « magnétiques » sont 
établies dans le cas où les moments magnétiques 
sont parallèles à une direction unique. Lorsqu'il 
n’en est pas ainsi, le carré du facteur de structure, 
rapporté à tout le cristal, est donné par la relation 
plus générale : 


(Fe? = 4(yro)? f(K)? 
| à (Si.K) K 
Z( K° 


où r, est le rayon classique de lélectron, y le 
rapport gyromagnétique du neutron, K le vecteur 
de diffusion, /(X) le facteur de forme magnétique 
et s le spin du neutron. La somme est étendue à 
tous les atomes 1 de spin S; dont la distance à 
l’origine est 4, 


2 


=è Si) .S exp (:K.R;) (2) 
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Pour des neutrons non polarisés cette formule … 


peut se mettre sous la forme : 


lFel? = (yro)? f(K KP 2 Vi.Vi exp (iK.(R; — R)j)) 
avec 
re (Six) ET 


Pour un hélimagnétique nous venons de voir ‘que 
tous les moments magnétiques sont parallèles à à un 
même plan P, et parallèles entre eux à l’intérieur 


d’un même feuillet perpendiculaire à l’axe quater- 
naire. Le calcul montre alors que |F.|? n’est diffé- 
rent de zéro que si : 


K = 7T=+2 


S 1-6 


20 


w étant le vecteur unitaire de la direction de l’axe 
quaternaire et æ un vecteur du réseau réciproque 
(défini à un facteur 27 près). Dans ces conditions, 
le carré du facteur de structure, rapporté cette fois 
à la maille cristalline et non plus à tout le cristal, 


s'écrit : 
»2 
ae = tros une se(1 + Éà) (3 
ne 
où XX est la projection de EÆ sur la normale au 
plan P. 


C’est dire que l'intensité ne sera différente de 
zéro que si le vecteur de diffusion a son extrémité 
en un point de l’espace réciproque se déduisant 
d’un nœud du réseau polaire par une translation 


(fig. 7). 


cléaire Akl sera donc accompagnée de deux satel- 
lites magnétiques que nous indexerons hkl* et hkt-, 
le premier correspondant, par convention, au vec- 
teur de diffusion le plus grand. Dans le cas d’un 
diagramme de poudre, chaque raie nucléairé hkl 


29 
égale à + — Lo Chaque réflexion nu- 


sera effectivement flanquée de deux satellites ma- 
gnétiques. Toutefois, les deux satellites d’une 
raie 40 seront confondus dans un tel diagramme, 


< 
1018 00 


101= 020- 


not 
110— 


n0$ 


RrGN 7: 


les vecteurs de diffusion qui leur correspondent 
étant égaux. Ce satellite unique sera noté hk0*. 
En particulier, une réflexion 000* doit correspondre 
au nœud origine du réseau polaire. 


3) STRUCTURE MAGNÉTIQUE DE MnAu,. — Le 
spectre différence représenté sur la figure 2 corres- 
pond justement à la description que nous venons de 
donner. On voit en particulier que chacune des 
deux raies nucléaires 002 et 101 est accompagné de 
deux satellites. Pour les raies plus lointaines, l’en- 
chevêtrement des raies est trop grand pour qu’il 
soit possible de les utiliser. À partir des satellites 
de 002, on calcule l’angle + dont les moments 
magnétiques tournent en passant d’un feuillet au 


suivant 
@ = 510 + 10, 


A l’aide de cette valeur on peut déterminer a 
priori la position des différents satellites. Ces posi- 
tions sont indiquées sur la figure 2b. Sur le spectre 
(fig. 2a), on ne voit pas la raie 000* qui devrait 
apparaître à 206’. Elle est beaucoup trop proche du 
faisceau direct pour pouvoir être observée sans 
précautions. Nous avons dû pour cela limiter le 
faisceau direct par un écran de cadmium ; un 
spectre différence entre 130 0C et la température 
ordinaire donne le résultat représenté sur la figure 8 
La raie observée est très intense, comme le montre 
la comparaison avec la raie 002 représentée à la 
même échelle, et apparait bien à l’angle prévu. 
Notons en outre qu'aucun harmonique de cette 
raie n’a été observé, ce qui semble exclure une 
structure à antiphase périodique qui donnerait 
également des satellites mais accompagnés d’har- 
moniques. La structure hélimagnétique de MnAu, 
est encore confirmée par la comparaison entre les 
intensités observées et calculées. Les intensités 
mesurées ont été portées en échelle absolue par 
comparaison avec la somme des intensités calculées 
des deux premières raies nucléaires, pratiquement 
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indépendante de la valeur du paramètre z. Le calcul 
de l’intensité des raies magnétiques a été fait en 
supposant, par raison de symétrie, que les moments 


000* 


| 


NTC-R9 


magnétiques étaient normaux à l’axe quaternaire, 
ce qui conduit à remplacer dans la formule (3) Æ# 
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par K, Pour déterminer la valeur du moment 
magnétique d’un atome de manganèse, nous nous 
sommes servis de la raie 000* qui est la plus intense 
et pour laquelle le facteur de forme magnétique 
est pratiquement égal à 1. On trouve u, — 3,04 ur, 
valeur un peu plus forte que celle (2,88) déduite 
par Meyer et Taglang des mesures magnétiques. 
Le tableau IT donne la comparaison entre les inten- 
sités observées et calculées en prenant pour facteur 
de forme magnétique celui de lion Mn++. Dans le 
calcul, nous n’avons pas introduit de facteur de 
Debye, ce qui n’apporterait qu’une faible correc- 
tion, les mesures n’allant pas à de grands angles. 
On voit sur le tableau IT que le rapport entre les 


TABLEAU II 


COMPARAISON ENTRE LES INTENSITÉS MAGNÉTIQUES 
CALCULÉES ET OBSERVÉES 


hkkl I cale. (barns) : I  obs. (barns) TT 
000€ 936 936 

002— 20,6 19,9 1,04 
002+ 9,8 10,3 0,95 
101— 13,0 : 

Ne 1e 24,5 34,1 0,72 
110+ 4,0 

004— 2,0%. 44,7 14,6 0,80 
103— 5,7 

112— 3,2 7,7 0,42 
103— 3,7 

004+ 1,1 ps 16,4 0,45 
112+ 2,5 


valeurs théoriques et expérimentales décroît quand 
l’angle de diffusion augmente. Ceci peut s’expliquer 
en admettant que le facteur de forme magnétique 
de l’atome de manganèse dans MnAu, décroît 
moins vite que celui de Mn++. La figure 9 donne 
la variation de ce facteur de forme, déduit de nos 
mesures, comparé à celui de Mn++. 


4) INFLUENCE DE LA TEMPÉRATURE SUR LA 
STRUCTURE MAGNÉTIQUE DE MnAu,. — Nous avons 
étudié l'influence de la température sur la structure 
magnétique. L’échantillon, placé dans un cryostat, 
était refroidi par un courant d’azote dont on pou- 
vait contrôler la température et le débit. La tempé- 
rature de l’échantillon était uniforme à 0,29 près. 
Au cours d’une expérience, la température de 
l'échantillon ne variait pas de plus de 2°. Nous 
avons vérifié que, jusqu’à 87 °K, les spectres de 
diffraction restaient qualitativement semblables, ce 
qui semble exclure une modification de la structure 
magnétique, qui reste hélimagnétique. La position 
de la raie 002 varie du fait de la dilatation ther- 
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mique (fig. 10). On en déduit la valeur du coefficient 


de dilatation & — 25 + 5 10—$ pratiquement égale 


à la valeur du coefficient de dilatation moyen déduit 
des mesures précédentes. La détermination la plus 
précise de la valeur de l’angle + entre les moments 
des atomes de deux feuillets voisins s'obtient à 
partir de la différence des angles de Bragg des 
raies 002— et 002 (?). La figure 11 montre que la 


28502 


15°05° 
15° 
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Diet 
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variation de l'angle © avec la température est 
faible : ®9 commence par décroître linéairement et 
passe par un minimum aux environs de 125 °K. 
Le fait que l’angle © varie peu avec la température 
est en accord avec le modèle hélimagnétique utilisé, 
car © ne dépend que des coefficients d’échange qui 
varient peu avec la température. L'existence d’un 
minimum pour @ montre que la variation de m; 
et #2 ne peut être linéaire en 7, ce qui n’a rien de 
surprenant, d'autant que ces coefficients sont des 


(%) La raie faible qui apparaît entre les raies 002+ et 


101— (fig. 2a) est la seule qui ne puisse s’interpréter par les: 


structures nucléaire et magnétique proposées. Son origine 
nucléaire est certaine : elle re disparaît pas au-dessus de 
90 °C. Elle doit être attribuée à ure impureté qu’il ne nous 
a pas été possible d'identifier. Malheureusement le fait 
qu’elle soit proche de la raie 002— empêche de pointer 
cette dernière avec une bonne précision. 


6 


moyennes pondérées entre plusieurs coeflicients 


d'échange (relation (1)) dont on ne connaît pas la 
variation avec la température. La mesure de l’in- 
tensité du satellite 002— permet également de déter- 
miner la variation thermique du moment magné- 


200 
Fre 12, 


tique (fig. 12). L’extrapolation au zéro absolu donne 
une aimantation de 3,6 u# par atome de manganèse. 


III. Action d’un champ sur la structure magné- 
tique de MnAu.. — Nous avons ensuite cherché à 
expliquer le mécanisme auquel MnAu, doit ses pro- 


I berns 


20 002 
10 002 + 
TEA 
ee H.K%e 
0 10 20 
Hrce19;: 


priétés métamagnétiques, mécanisme dont nous ne 
donnerons qu’une théorie approchée. Expérimen- 
talement, ce mécanisme se manifestera par l’in- 
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fluence d’un champ sur l’intensité des raies inagné- 
tiques. Une première information intéressante est 
obtenue en étudiant la variation de l’intensité des 
raies 002* et 002 en fonction d’un champ appli- 
qué parallèlement à l’axe quaternaire. La figure 13 
montre que l’intensité de ces raies varie au maxi- 
mum de quelques pour cent pour des champs allant 
jusqu’à 20 kOe. Ceci montre que pour des champs 
de cet ordre, la direction des moments s’écarte très 
peu du plan 001. Une telle rigidité doit être attri- 
buée à un champ cristallin très intense, au moins 
égal à 100 kOe. L’existence d’une aussi forte aniso- 
tropie, et la structure magnétique elle-même 
excluent un mécanisme d’aimantation analogue à 
celui d’un métamagnétique du type FeCl, [9]. 


1) THÉORIE. — Dans ces conditions, lorsqu’on 
applique un champ magnétique à un cristal de 
MnAu, d’une intensité inférieure à quelques di- 
zaines de kOe, seule la composante de ce champ 
située dans le plan perpendiculaire à l’axe quater- 
naire pourra agir sur les moments magnétiques. Il 
nous suffit donc d’étudier théoriquement l’action 
d’un champ perpendiculaire à l’axe c sur la dispo- 
sition des moments. On trouve le résultats suivant : 
selon la valeur du champ, les moments peuvent 
présenter deux configurations, l’une voisine de la 
configuration hélimagnétique, l’autre voisine de la 
configuration ferromagnétique, rigoureusement 
atteinte pour un certain champ que nous appel- 
lerons champ limite. On obtiendra la première solu- 
tion en étudiant la susceptibilité initiale, la seconde 
en étudiant le champ limite et l’approche à la 
saturation. 

Nous appellerons 0, l’angle des moments des 
atomes du piè% feuillet normal à l’axe c avec la 
direction du champ. L'expression générale de 
l’énergie est alors, à un terme constant près, pro- 
portionnelle à 


Wie Dico 10 0, icon 0,7) 
? 
— po AH cos 6» (4) 
? 
OÙ 7 = 74 LÉ 
1H Fo: 

a) Susceptibilité intiale. — Dans un champ 
faible les moments sont peu déplacés par rapport à 
la configuration hélicoïdale, ce qui s’écrit : 

Op = 00 + PP + © (5) 
où & est petit. Dans ces conditions, l'équation 
d'équilibre pour le moment pm, obtenue en mini- 
misant l’énergie par rapport à 6, s’écrit d’après (4) : 
— ND COS P(2€p — Ep—1 — Ep+1) 


— N COS 2 p(2Ep — Ep—2 — Ep+2) = Lo H Sin (05 + po) 
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équation dont la solution est 
ep — À sin (0, + p9) (6) 
avec 
A = BE Bo H 
4 (ra cos ® sin? + + n9 COS 26 sin? +) 
ne ee | sut 
72 sin & (1 + 2 cos p (1 + cos æ)) 


La valeur moyenne de la projection du moment 
dans la direction du champ est : 


BE = bo COS (09 + po + Ep) 
———— 1 
= — Ho €p Sin (0 + po) = — 3 Uo À. 


On en déduit la susceptibilité par atome : 


u0 1 
MEET 9 s : 
sin + (1 + 2. cos æ (1 + cos œ)) 


L’approximation que nous faisons, et qui se 
trouve largement justifiée par un calcul plus poussé, 
est de supposer que l’aimantation varie linéai- 
rement avec le champ, quel que soit le champ : 


My = — 2€ sin (0, + po 
0 
le pas de l’hélice gardant la même valeur. On en 
déduit lexpression de l’énergie, rapportée à un 
atome : 


W; = —n, cos ® — n, cos 2e — 5 x = Wo—xe. 


b) Champ limite. — Dans un champ très intense 
tous les moments magnétiques sont parallèles au 
champ. Si l’on fait décroître le champ, les moments 
s’écartent de la direction du champ lorsque celui-ci 
devient inférieur à un certain champ limite H, dont 
on peut calculer l'expression. 

Soit », l'angle, supposé petit, du moment magné- 
tique p, avec le champ H ; l'équation d’équilibre 
du moment p, s'écrit, d’après (4), pour les petites 
valeurs de "y : 


— Na(2p — Vw—1 — Np+1) 


— afp — Tp—a — Tpta) = Uo Hp. 


L'ensemble de ces équations forme un système 
d'équations linéaires homogènes qui n’admet que la 
solution 1, = 0, sauf si w, À est une valeur propre. 
Le champ limite , sera donc la valeur propre la 
plus grande de ce système dont les fonctions 
propres sont : 


“= C sin (pŸ + à). 


La valeur propre correspondant à à est 


Uo H(Y) = — 2n,(1 — cos Ÿ) 
| — 2n,(1 — cos 24) = — 2{ns + ne) — Wo(Ÿ). 
où W,(4) est l’énergie du système hélimagnétique, 


l'angle des moments de deux feuillets voisins 
étant 4. H,(Ÿ) est donc maximum quand W,(Ÿ) 


‘est minimum, c’est-à-dire lorsque 4 est égal à 


l’angle o défini par cos @ — —n,/4n,. On a donc : 
1672. 
He TM (8) 


c) Approche à la saturation. — Nous allons main- 
tenant calculer l’aimantation pour des valeurs du 
champ voisines de À,, c’est-à-dire la susceptibilité 
différentielle au voisinage de la saturation. Dans 
le paragraphe précédent nous avons écrit l’équa- 
tion d'équilibre du moment », sous la forme linéaire 
valable seulement pour des angles %, infiniment 
petits. Sinon cette équation s'écrit : 


— n3 (Sin (np — Np—1) + Sin (Np — Mp+1)) 


— A3 (Sin (Mp — Np—2) + Sin (np — Np+2)) = Lo Sin y. (9) 


La non-linéarité de cette équation entraîne deux 
conséquences : une limitation de l’amplitude », et 
l'apparition d’harmoniques, ce qui est habituel 
pour les systèmes non linéaires. Les termes anhar- 
moniques étant des fonctions impaires de %, seuls 
apparaîtront les harmoniques impairs et le dévelop- 
pement de , S’écrira : 

np = Qsin (po + 8) + CG sin (3pp + à’) + ... (10) 
où C1, C3, ... sont des fonctions de 7, — Æ. Nous 
nous limiterons au calcul du premier coefficient 
C.(4) qui est de degré le plus bas en 4, — H, et 
même au premier terme de ce développement. 
Pour cela, utilisons un développement de l’équa- 
tion (9) 


— Na (2%p — Tps — Vp+ 1) 
— Nan — Tp—2 — Tp+2) do Ho TD 
n , < 
Sa G [ne — Vo) + (ro — Yp+3)°] 


+ U [ro — Np—2) + (rp — p+2)°] 


Uo H 
EN = dite =) 


Si l’on remplace », par €, cos (po + à), la pre- 
mière ligne s’annule identiquement par définition 
de 11, ; la seconde ligne nous donne des termes en 
Ci cos (pp + à) et des termes en € cos 3(po + 5), 
ces derniers devant être compensés par les termes 
provenant de €; cos (3po + à’). Ceci montre tout 
d’abord que le second coefficient €;, proportionnel 
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à Ci est d’un ordre supérieur en J,— H. En éga- 
lant les termes en cos (po + 3), on trouve : 


CG f(@) = UolHo — À) 6; 


avec 
fo) = HT (2 4 (1 + 2 cos 9) (11) 
d’où : 
ta = V BotHo — H) (12) 


f(œ) 


L’aimantation correspondante est donc 


re fl il 
M = {9 COS np — bo (1 ne D) FE do (1 —; 4) 
ou , 
me y —#$ 


On voit donc qu'au voisinage de la saturation, 
Paimantation varie linéairement avec le champ 
magnétique. Le calcul plus précis montre que cette 
approximation reste bonne, même pour des valeurs 
de 4, — H qui ne sont plus petites. En nous limi- 
tant à ce premier terme, l’énergie s’écrit : 


Wa = bo Ho 
ds cos . H 2 ( 1) 


1 — 
9 / 2 
16 sint ? H5 2 + (1 + 2 cos o) 


d) Structure magnétique de MnAu, en présence 
d’un champ. — D’après ce qui précède nous pou- 
vons nous représenter l’action sur MnAu, d’un 
champ normal à l’axe c de la manière suivante : 
lorsque le champ croît à partir de zéro, la confi- 
guration hélimagnétique de départ se déforme légè- 
rement. Pour une certaine valeur du champ, que 
nous appellerons le champ seuil, la configuration 
. des moments doit changer brusquement et devenir 
voisine de la configuration ferromagnétique. Quand 
le champ continue d’augmenter, les moments 
tendent à être parallèles, et le sont dès que le champ 
atteint la valeur du champ limite. 

Le problème que nous devons résoudre mainte- 
nant est le suivant : quelles sont les configurations 
stables pour les différentes valeurs du champ — 
ce qui revient à rechercher les lois de variation 
0, = 0,(4) — et en particulier quelle est la valeur 
du champ seuil. 

Les premiers termes du développement de 0,(/7) 
pour les deux configurations possibles ont été 
donnés aux paragraphes a) (formules 5, 6, 7) et 
c) (formules 10, 11, 12) précédents, et peuvent être 
utilisés dans une première approximation. On peut 
aussi, à cette même approximation, calculer le rap- 
port 4 du champ seuil au champ limite en égalant 

0 
les valeurs des énergies correspondant aux deux 
configurations. Si l’on pousse les développements 
aux termes d’ordre supérieurs en À, on voit que 


L'ANTIFERROMAGNÉTISME HÉLICOIDAL DE MnAu, 


345 


ces termes décroissent si rapidement que la pre- 
mière approximation est largement suffisante. Nous 
nous en sommes encore assurés par un calcul numé- 
rique : la valeur de ® étant très voisine de 2x/7, 


W X \ 


Re 


sue 


o 


D2SND7E RTC NEN0S 


(2 


IVe 


il est possible, avec cette valeur, de résoudre, grâce 
à un ordinateur électronique, le système des sept 
équations 0W/d9, — 0 ({). Le résultat confirme ce 
que nous attendions. 


(4) Nous remercions la Section de Calcul Électronique 
Arithmétique du Centre d'Études Nucléaires de Saclay à 
qui nous devons ces calculs. 
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La figure 14 donne ce résultat. On voit que pour 
H, — 0,485 H, l'énergie de la deuxième configu- 
ration devient inférieure à celle de la première. Il y 
a donc pour cette valeur du champ un saut brusque 
de la première configuration à la seconde. L’arran- 
gement des moments magnétiques des sept feuillets 
successifs est représentée à la partie inférieure de la 
figure pour À = 0, H;, H,. 

La figure 15 représente l’aimantation d’un feuil- 
let en fonction du champ appliqué dans le plan 
du feuillet. La comparaison avec la courbe d’aiman- 
tation expérimentale d’un échantillon polycristallin 
(fig. 1) suggère une valeur de /, de l’ordre de 
10 kOe. 


2) ÉTUDE PAR DIFFRACTION DE NEUTRONS DE 
L'INFLUENCE DU CHAMP SUR LA STRUCTURE MAGNÉ- 
TIQUE DE MnAu,. — Si le mécanisme que nous 
venons de proposer est exact on doit s'attendre à 
modifier considérablement l’intensité d’une raie 
satellite lorsque la composante du champ dans le 
plan des moments magnétiques est importante. En 
effet quand cette composante atteint la valeur du 
champ limite la raie satellite doit disparaître, la 
saturation ferromagnétique étant atteinte. On peut 
donc penser que, quand cette composante attemdra 
la valeur du champ seuil, pour laquelle la confi- 
guration quasi-hélimagnétique passe à la confi- 
guration quasi-ferromagnétique, l'intensité de la 
raie subira une chute brutale. Ainsi que nous l’avons 
déjà signalé, la courbe d’aimantation suggère une 
valeur du champ seuil de l’ordre de 10 kOe, de 
sorte que l’on peut espérer mettre aisément en 
évidence ces variations d'intensité, pourvu que le 
champ appliqué fasse un angle « assez faible avec 
le plan des moments magnétiques. 

En disposant le champ parallèlement au vecteur 
de diffusion, on ne peut espérer observer cet effet 
sur les raies 002* et 0027, le champ étant alors 
normal au plan des moments magnétiques. C’est 
bien ce que nous avons observé. Par contre, avec la 
même disposition du champ par rapport au vecteur 
de diffusion, la situation est beaucoup plus favo- 
rable pour les satellites 101 et 101—, l’angle « 
entre le vecteur de diffusion (donc le champ), et le 
plan des moments étant petit dans les deux cas. 
La composante Æ cos « du champ dans le plan des 
moments est en effet égale à 0,964 Æ pour 101— et 
0,896 Æ pour 101*. La figure 16 donne le résultat 
de cette expérience. On voit que l’intensité de ces 
raies subit bien une décroissance rapide dès que le 
champ atteint 10,2 kOe pour 101—-eét 10,7 kOe 
pour 101*. 

Mais il est également possible de mettre en évi- 
dence le changement de configuration des moments 
dans leur plan sur les satellites 002— et 002*, si lon 
dispose le champ non parallèlement au vecteur de 
diffusion, mais normal au plan de diffusion, le 
champ étant alors parallèle au plan des moments. 
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La figure 17 donne le résultat de l’expérience sur le 
satellite 002—. Le satellite 002* donne un résultat 
absolument identique. La valeur du champ pour 
laquelle s’ amorce la chute d'intensité est égale cette 
fois à 9,5 kOe. Cette valeur doit être considérée 


O' 2-46. 8 10 2.W 16 ‘B120 -H°Kœ 
F1G. 16. 


H. KO 


—i— 


DRM NENEEn 


HG 13: 


comme la valeur du champ seuil Æ, et doit en 
principe coïincider avec les valeurs des projections 
H, cos « des champs seuils «apparents » A; obtenus 
pour les réflexions 101— et 101*, Le tableau III 
montre qu'il en est bien ainsi. La valeur moyenne 


TABLEAU III 


DÉTERMINATION DE LA VALEUR DU CHAMP SEUIL 


hkl 002% 10e To 
Hs 9,9 »10,7 10,2 
UE 9,9 9,6 9,8 


N°6. 


L’ANTIFERROMAGNÉTISME HÉLICOIDAL DE MnAu, 


bon accord avec celle que laissait prévoir la courbe 
d’aimantation. | 
On peut en outre, avec la même disposition du 
champ, avoir la preuve qu’on tend bien, dans les 
champs forts, vers une saturation ferromagnétique. 
Il suflit pour cela d'étudier l'influence du champ 
sur l'intensité de la réflexion 002, car, lorsque la 
saturation ferromagnétique est atteinte, à l’inten- 
sité nucléaire de cette réflexion doit s’ajouter une 
composante magnétique qui, le calcul le montre, 


N° 6 
. du champ seuil qu’on en déduit, 9,6 kOe, est en 
doit atteindre 46 % de la composante nucléaire. 
| À 
| 


Lien 


002 


. La figure 18 montre l’effet observé qui confirme ces 
prévisions. 

On peut rendre théoriquement compte de la va- 
riation d'intensité des satellites magnétiques avec 
le champ observé. Nous supposerons toujours que 
les moments magnétiques ne peuvent tourner que 
dans le plan des feuillets de manganèse, et que la 
périodicité de la structure magnétique correspond 
à NN feuillets. En appelant 6, l’angle entre les mo- 
moments des atomes du feuillet 7 avec la direction 
de la projection Æ/' du champ A sur le plan du 
feuillet, une application de la relation (2) donne 
aisément l’expression du facteur de structure 
pour A feuillets : 


= \K)" (yro)" 5? 


KT A . 2 Ne RCE Di 
ip cos ds cos je | ls sin din je | ë 


Æ? Liz j=1 
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D'où le rapport de l’intensité de la raie pour la 
valeur Æ du champ à l’intensité en champ nul : 


AN pe y se 
HAE AT EU 
1(0) De 

Ne +1 


avec 


N 
T =.» cos 0;/H) cos je 
1 


î— 


N 
E — }ÿ, sin 0;({) sin jo. 
j=1 


Ce rapport a été calculé numériquement dans les 
trois cas correspondant à nos expériences sur 002, 
101— et 101*. La figure 19a donne le résultat de ce 


I calculée I observée 
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calcul. La comparaison avec la figure 19b où sont 
superposées les trois courbes expérimentales montre 
que l’expérience est en bon accord avec les prévi- 
sions. Le fait que la chute d'intensité pour la 
valeur du champ seuil soit plus progressive que 
prévu peut s’interpréter par l’existence d’un champ 
local équivalent à un champ démagnétisant, qui 
serait d’après les courbes expérimentales, de l’ordre 
de 5 kOe. 

Si l’on se reporte maintenant aux grandeurs que 
nous avons introduites au cours des calculs, nous 
voyons qu’il est possible de les évaluer à la tempé- 
rature ordinaire où l’expérience nous a donné 
H; — 9,6 kOe. La résolution numérique des équa- 
tions W, — W, nous ayant conduit à H, — 0,485 H, 
on en déduit la valeur du champ limite : 
H, = 19,8 kOe. On peut ensuite, par application 
de la relation (8) obtenir la valeur des coefficients 
d'échange : 

D Hip 1:02 10 2'erg 
ICS) 
ce qui équivaut à une température de — 7,3 °K, 
et 


ny =) na COS = 0 2,06 10775 erg: 1(18,2:9K). 
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3) FORME DE LA COURBE D’AIMANTATION POUR 
UN POLYCRISTAL. — Par diffraction de neutrons 
nous avons pu obtenir les valeurs du champ seuil 
et du moment magnétique des atomes de manga- 
nèse à la température ordinaire. Il est alors possible 
de calculer en valeur absolue, avec le modèle que 
nous avons utilisé, l’aimantation d’un échantillon 
polycristallin en fonction du champ appliqué, afin 
de vérifier si le résultat est en accord avec l’expé- 
rience (fig. 1). Dans ce calcul, on suppose toujours 
que la seule composante du champ susceptible de 
faire varier l’aimantation d’un cristallite est la 
composante normale à l’axe c de ce cristallite. L’ai- 
mantation varie avec cette composante suivant la 
loi représentée sur la figure 15, et lui est constam- 
ment parallèle. Cette loi, nous l’avons déjà signalé, 
est représentée avec une très bonne approximation 
par 


ï a 1 H 
M, x si H < H; 
M : AT 
0 
M £ 
Us STREE>H6 


En supposant les cristallites orientés au hasard, 
on trouve aisément : 


sin—1(H /H) 
7x4 u sin? 0 d0 
sin—1{H,/H) 
T De 
Le in? 640 
F Le (Æo/A) Rues ; 


Les termes successifs de cette somme corres- 
pondant, dans l’ordre, aux contributions des cris- 
tallites orientés de telle sorte que la composante du 


[(4— x) H9 + x’ Hsin 0]sin?0 40 
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champ normale à leurs axes c soit inférieure à 4, 
comprise entre /7, et Æ,, et supérieure à A4. 


Le calcul numérique conduit à la courbe repré- 


sentée en trait plein sur la figure 20. La compa- 


G uem/g 


Calculé 


(ko=304p, 3 Ho=198 Ke) 


FrG020: 


raison avec la courbe expérimentale reportée sur la 


2 


i 
, 


même figure, en tirets, confirme la valeur du. 


modèle proposé. 


La croissance trop lente de la courbe expéri- 


mentale, immédiatement après le champ seuil, est à 
rapprocher de la chute d’intensité progressive 
observée pour les intensités des satellites magné- 


tiques. Le fait que la susceptibilité initiale soit 


inférieure à la susceptibilité observée peut s’expli- 
quer par l’existence d’un paramagnétisme constant 


dont nous n’avons évidemment pas tenu compte. . 


En outre, aux champs élevés, il est certain que la 
courbe théorique doit remonter, dès que le champ 
devient assez intense pour faire sortir les moments 
magnétiques du plan des feuillets. 


Manuscrit reçu le 4 janvier 1961. 
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ÉTUDE DES RÉACTIONS PRODUITES PAR DES IONS DE Li SUR LE BÉRYLLIUM 


Par C. LEMEILLE, L. MARQUEZ (*), N. SAUNIER, 
Section de Physique Nucléaire à Basse Énergie, CG. E. N., Saclay. 


et M. COSTE, 


Laboratoire de Physique Nucléaire, Faculté des Sciences de Bordeaux. 


. Résumé. — Nous avons étudié les réactions nucléaires produites par les ions de Li sur le béryl- 
lium, avec énergie variable jusqu’à 2 MeV. Nous indiquons les résultats obtenus avec des scintil- 
lateurs, des plaques nucléaires et des jonctions. Les mécanismes possibles des réactions observées 


sont discutés. 


Abstract. — We have studied the nuclear reactions produced by fLiions of energy up to 2 MeV 


on Be. 
Possible reaction mechanisms are discussed. 


I. Introduction. — Nous avons utilisé le Van de 
de Graaff de 2 MeV de Saclay pour accélérer, selon 
la technique de Allison et ses collaborateurs [1], des 
ions Li. Dans un premier travail, nous avons 
étudié les réactions nucléaires produites par ces 
ions sur une cible de °Be [2]. 

Nous avions trouvé des raies y de forte intensité 
qui correspondaient au premier état excité de 7Li 
. de 479 keV et au premier état excité de !°B de 
717 keV. Nous avions également mis en évidence 
des y de plus faible intensité. 

Dans les lignes qui suivent nous exposons les 
résultats obtenus dans la suite de nos expériences. 


IT. Observation du spectre y avee un seintil- 
lateur. — Un détecteur muni d’un scintillateur 
Nal(Ti) de 3,8 cm de diamètre et de 5 cm de hau- 
teur nous a permis d’enregistrer le spectre de 
rayons y émis sur un sélecteur R. I. D. L. à 100 ca- 
naux, précédé du dispositif électronique usuel. Le 
compteur était placé à 00 par rapport au faisceau 
et à 3 cm d’une cible épaisse de béryllium. En 
faisant varier le gain, nous avons observé diffé- 
rentes raies y d'énergies plus fortes. Les résultats 
obtenus sont notés sur les figures 1 et 2. 

La raie de 1 010 keV provient probablement de 
la transition entre le deuxième et le premier état 
excité de Be. Celle de 2 110 keV peut provenir, 
soit du premier état excité de 11B, soit du troisième 
état excité de 1°B, soit du premier état excité de 
13C de 3 090 keV, car 
hv, — 2mc? = (3 090 — 1 022) keV = 2 068 keV. 

Nous observons une raie de 2 700 keV qui pour- 
rait être un mélange de Compton et du pic de 
3 090 keV avec une perte d’un seul photon d’anni- 


(*) Chercheur de la Comissao Nacional de Energia 
Nuclear. Professeur du Gentro Brasileiro de Pesquisas 
Fisicas, en stage au Centre d'Études Nucléaires de Saclay. 


The results obtained by scintillation counters, nuclear plates and junctions are presented. 


hilation. La raie trouvée à 3 160 keV semble corres- 
pondre au pic de 3 090 keV. Les intensités des 
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Fic. 1. — Spectre de rayons y de 300 keV à 4 000 keV 
provenant de la réaction du $Li de 2,0 MeV sur ‘Be. 


autres raies y sont, compte tenu des efficacités 
relatives, au maximum quatre fois plus faibles que 
celle de 479 keV. Les intensités des raies de 479 keV 
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et 717 keV ont été mesurées en fonction de l’éner- 
gie de 6Li. Ces mesures ont été faites par inter- 
valles de 0,05 MeV dans un domaine d’énergie 
compris entre { et 2 MeV. Nous avons utilisé une 
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F1G. 2. — Spectre de rayons y de 1 000 keV à 7 000 keV 
provenant de la réaction du Li de 2,0 MeV sur ‘Be. 


cible mince de béryllium déposé sur tantale, Cette 
cible était placée dans une chambre reliée à un inté- 
grateur indiquant le nombre de microcoulombs dé- 
posés. Les courbes de la figure 3 sont tracées à 
partir des résultats obtenus. La cible utilisée étant 
mince, les ordonnées sont proportionnelles à la 
section efficace. Nos résultats sont d’autre part en 
accord avec les courbes théoriques fournies par le 
calcul du facteur de Gamow. Pour la raie de 
479 keV, dans le cas de l’irradiation d’une cible 
épaisse par le $L1 de 2 MeV, le rendement est de 
l’ordre de 5 X 104 gamma/microcoulomb de Lit, 
l'erreur commise étant de 50 %. 

Nous avons mesuré le rapport des intensités des 
raies de 479 keV et 717 keV pour différentes valeurs 
de l’énergie de 6L1. La mesure de ce rapport a été 
faite plusieurs fois pour chaque point avec une 


N°6 


bonne précision, les deux raies étant enregistrées | 
sur le même spectre et au cours de la même irra- 
diation. Le tableau I a été établi à partir des, 
résultats obtenus après correction d’efficacité, . 
celle-ci restant la même pour toutes les mesures. 
Nous remarquons que le rapport considéré est, 
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Fic. 3. — Rendement des raies de 479 keV et 717 keV en 
fonction de l’énergie de 6Li sur une cible mince de ‘Be. 


constant dans l’étroite limite des erreurs expéri- 
mentales. 

Nous avons également mesuré l’effet Doppler- 
Fizeau produit sur la raie de 479 keV. Pour effec- 


TABLEAU I 


RaProrT 


ÉNERGIE Éa e DÉ VIATION. 

pu 6Li (MeV) DATE STANDARD 
1,1 93 19 
1% 1:05 06 
1,3 99 "04 
1,4 1,01 02 
15 96 702 
1,6 ,98 ,02 
1,7 97 02 
18 198 702 
1,9 1,01 02 
2.0 1,05 04 
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tuer cette mesure, le compteur de Nal a été placé 
à 00 et à 900 par rapport au faisceau et à 4 em de 
la cible. On a trouvé un déplacement de la raie 
positif avec Br — (0,77 + 0,14) X 102. L'énergie 
moyenne des 6Li qui produisent la réaction étant 
de 1,8 MeV, la vitesse du centre de masse est donnée 
par B, = 1,02 X 10—2, Donc en moyenne, le 7li 
dans le centre de masse sort vers l’arrière avec une 
faible composante moyenne de la vitesse dans la 
direction contraire du faisceau et telle que 


Be = — (0,25 + 0,14) X 10%. 


III. Observation des particules chargées au 
moyen de plaques nucléaires. — L'étude des parti- 
cules chargées provenant de la réaction a été faite 
au moyen de plaques nucléaires et de jonctions. 
Pour déterminer la réaction qui produit le ?Li 
et le 10B, nous pensions examiner les particules 
chargées émises préférentiellement dans le demi- 
cône avant avec une énergie de quelques MeV. 
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Fic. 4, — Histogramme des particules chargées provenant 
de la réaction de SLi de 2,0 MeV sur °Be, observé à 420. 
Le premier pic est celui de 4He et le deuxième celui 
de $He de la réaction Li + 1H — 4He + °He. 


La réaction 6Li + 1H — $He + ‘He qui se 
produit à basse énergie nous a empêchés d’effec- 
tuer les observations dans cette région. L’hydro- 
vène qui provoque cette réaction pourrait pro- 
venir des matières organiques des joints et de la 
pompe à diffusion. Les intensités des pics °He et 
4He sont environ 10 fois plus élevées que celles 
des réactions recherchées. La figure 4 représente 
l’histogramme obtenu avec une plaque à 420. 
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IV. Observation des particules chargées au 
moyen de jonetions. — [1 est possible d’effectuer 
des mesures dans le domaine de 909 à 1800 sans 
être gêné par la réaction due à l’hydrogène. Des 
observations ont été faites dans cette région avec 
un compteur à jonction. 


DISPOSITIF EXPÉRIMENTAL. — Nous avions à 
notre disposition une jonction R. C. A. au silicium 
type p de 20 mm°, et le matériel électronique 
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Fc. 5. — Spectre des particules chargées provenant de la 
réaction de 6Li sur °Be, observé à 1202. Les pics d’alpha 
indiqués sur la figure proviennent de la réaction 

6Li + Be — x + 11B, 


courant : un préamplificateur gain 25 légèrement 
modifié, suivi d’un amplificateur proportionnel 
2 Me type C. E. A. La jonction était placée à 5 em 
de la cible et son fonctionnement sous vide pen- 
dant un mois ne l’a pas détériorée. La position du 
préamplificateur à l’extérieur de la chambre à 
diffusion et la géométrie de cette dernière nous 
ont imposé une liaison capacitive de 10 à 20 pF ; 
d’où une perte de résolution du détecteur que nous 
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pensons éviter au cours des prochaines expériences 
en utilisant un préampli de charge. 

Avec une source naturelle d’alpha (Th ©”) et 
pour une tension de 45 V appliquée à la jonction, 
nous avons obtenu des raies dont la largeur à mi- 
hauteur (l), indépendante de l’énergie des parti- 
cules détectées dans le cas d’une jonction, est de 
135 keV, en connectant directement la jonction au 
préamplificateur. Dans les conditions de nos expé- 
riences : [ — 300 keV. Cette résolution a été obte- 
nue sans aucune précaution spéciale. Elle est très 
loin des conditions optimum des jonctions, mais 
elle a été suffisante pour nos mesures. 

Sur la figure 5, qui représente le spectre des 
particules chargées observé à 1200, il est possible de 
séparer les pics «, et «, auxquels correspondent des 
énergies différant de 380 keV ; pour «, le F mesuré 
est de 340 keV, les raies pouvant d’ailleurs pré- 
senter un élargissement supplémentaire dû à la 
variation relativement rapide de l’énergie des par- 
ticules détectées en fonction de l’angle. Avec le 
cristal de CsI que nous possédions, la meilleure 
résolution obtenue était de l’ordre de 20 %. 


RÉSULTATS EXPÉRIMENTAUX. — Les raies mar- 
quées &op, A1, >, 3 Ont été identifiées et corres- 
pondent à l’état fondamental et aux trois premiers 
états excités de B. La réaction qui produit les 
alpha est la suivante : 

SLi + °Be = ‘He + 1B. 

Ces quatre groupes d’alpha ont déjà été trouvés 
et étudiés par Leigh et Blair [3] dont nous avons 
vérifié les résultats. 

Dans la région de basse énergie, plusieurs 
raies plus intenses, dont lidentification n’est pas 
encore définitive, apparaissent. 


V. Interprétation des résultats. — Parmi les ré- 
sultats obtenus, 1l semble intéressant de souligner 


ceux qui concernent la production intense de Li . 
et de 1°B et le rapport constant des sections effi- | 
caces de formation des premiers états excités de Li . 
et de 1°B en fonction de l’énergie du fLi incident. . 
Ces données expérimentales nous conduisent à deux 
hypothèses : le 7Li et le 1°B-proviennent, soit de la 
désintégration d’un seul état excité de B, soit 
de la désintégration des niveaux différents de 1B 
dont les rapports de désintégration doivent être 
identiques. Nous pouvons supposer, d'accord avec 
ces hypothèses, que la réaction primaire se pro- 
duit par stripping du deutéron de Li, ce dernier 
étant alors constitué d’un deutéron et d’un alpha. 
Dans ce cas, la réaction qui se produit est : 


SPiTEMBeE EEE; 


Le He emporte une énergie de 1 à 2 MeV. 
Le 11B est formé dans un état très excité et se 
désintègre en &« + 6Li ou n + 1#B. L’énergie dis- 
ponible dans notre réaction étant de 14,35 MeV, 
nous pouvons atteindre les niveaux de 4B de 
13,16 MeV, 14,0 MeV et 15,1 MeV. Ces trois niveaux 
se désexcitent par émission de « et n [4]. 

Dans notre première communication [2], nous 
avions envisagé la possibilité de l’existence des 
réactions de transfert d’un seul nucléon ; nous 
poursuivons nos expériences afin d'obtenir des ré- 
sultats complémentaires qui préciseront la contri- 
bution relative des deux mécanismes suggérés. 

Nous remercions M. le Pr Chastel pour la colla- 
boration apportée par le Laboratoire de Physique 
Nucléaire de Bordeaux, MM. A. Claes et J. Petres 
pour la conduite du Van de Graff, M. J. Quidort 
pour la mise en route du détecteur à jonction, 
Mme Bouraoui et Mlle Garrabos pour le dépouil- 
lement des plaques. 

Nous remercions particulièrement M. E. Cotton 
qui nous à encouragés au cours de nos expériences. 
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INFLUENCE DE LA TEMPÉRATURE SUR LE SPECTRE RAMAN DE LA GLYCÉRINE 


Par Ro8erT EBENSTEIN et Monique HARRAND, 


Laboratoire des Recherches Physiques, Sorbonne, Paris. 


Résumé. 


Les auteurs ont étudié l'influence de la température sur le spectre Raman de la 


glycérine. Ils n’ont pas décelé de variation de fréquences quand la température passait de 25° à 
170 °C, mais ils ont observé une diminution de largeur des deux bandes à 4 465 cm—t et à 1 066- 
1105 cmt (respectivement 26 % et 6 %), cette diminution résultant de la dépolymérisation. 


Abstract. — The authors have studied the influence of temperature on the Raman spectrum 


of glycerin. 


They have not seen any shift of frequencies when the temperature is increased from 


259 to 170 °C, but they have observed a decreasing width of the two bands at 1 465 em—t and at 
1 066-1 105 em! (respectively 26 % and 6 %) ; this decrease is a result of depolymerization. 


Pour obtenir un spectre Raman convenable de la 
glycérine il est nécessaire d’éliminer le fond continu 
de fluorescence qui se superpose aux raies Raman. 
Cette fluorescence est due pour une part aux impu- 
retés et pour une autre part à la décomposition de 
la glycérine sous l'effet des radiations ultra-violettes 
et violettes. 

Nous avons excité la glycérine par la radiation 
4 358 À du mercure et nous avons éliminé les radia- 
tions de plus courte longueur d’ondes en utilisant 
un double filtre formé d’une solution aqueuse satu- 
rée de NO,Na et d’un verre Schott GG2. Lors de 
Virradiation à température élevée, il a été néces- 
saire de renouveler la glycérine au cours d’une pose 
à cause de l’altération produite par le chauffage et 
par la décomposition photochimique. 

Nous avons utilisé la glycérine bi-distillée pour 
analyse livrée par plusieurs laboratoires. Le produit 
Prolabo, de densité 1,26 contient jusqu'à 3 % 
d’eau en masse ; la glycérine étant hygroscopique, 
cette teneur représente un minimum ; pour éviter 
la fixation de vapeur d’eau atmosphérique, les 
remplissages des cuves Raman ont été effectués en 
atmosphère anhydre, dans une enceinte constituée 
par un sac plastique transparent contenant de 
Vactigel. La glycérine Merck est d’une plus grande 
pureté, mais sa densité est 1,23 car elle contient 
jusqu'à 12 % d’eau. Un tel pourcentage était 
inacceptable car nous désirions obtenir le spectre 
de la glycérine et non celui d’un mélange ; d’autre 
part, lors des poses à chaud, l’eau en solution est 
libérée et forme des bulles qui perturbent l’image 
spectrale. Nous avons procédé à une deshydratation 
par chauffage modéré sous pression réduite : la 
température restant inférieure à 100 0C, la pression 
était abaissée à 12 mm de mercure. La teneur en 
eau’a été ainsi ramenée à 4,25 %. 

Les spectres ont été réalisés avec un montage 
Raman transversal ; la cuve était éclairée par deux 
arcs à vapeur de mercure. Le spectrographe com- 
porte deux prismes, il est ouvert à F/4,5 ; dans la 


région spectrale utilisée, la dispersion est de 
28 À/mm ; la largeur de la fente d’entrée était 
réglée à 0,1 mm. 

Nous avons effectué des poses à des températures 
comprises entre 20 0C et 200 0C. La cuve Raman 
étant disposée dans un four parallélépipédique 
étroit percé de trois fenêtres et placé entre les con- 
denseurs. Un système de réglage par thermomètre à 
contact et relais permettait de maintenir la tempé- 
rature à la valeur désirée. Pour éviter que les 
vapeurs libérées lors du chauffage ne se condensent 
sur les surfaces optiques, nous avons relié la cuve 
Raman à un tube manométrique contenant, lui 
aussi, de la glycérine. 

Les raies Raman de la glycérine sont assez larges 
et pour les spectres obtenus à température élevée le 
contraste diminue du fait de l’augmentation du 


fond continu de fluorescence. Pour cette raison les 


mesures de fréquences ont été faites sur les enregis- 
trements microphotométriques, par interpolation 
linéaire entre deux raies du fer. Nous avons évalué 
à 5 cm1 la limite supérieure d’une détermination 
de fréquence. Nous avons pris la moyenne de plu- 
sieurs mesures et pour chaque raie nous avons cal- 
culé l’écart moyen qui est, sauf dans un cas, nette- 
ment inférieur à 5 em” 1. 


Infiuence de la température sur les fréquences. — 
Les molécules de glycérine, qui comportent trois 
groupements OH, sont associées par des liaisons 
hydrogène. Pour connaître l’effet de la polyméri- 
sation sur le spectre Raman, il est nécessaire de 
trouver des conditions expérimentales pour les- 
quelles le degré d’association est inférieur. Dans ce 
but Saxena [1] a étudié les spectres des mélanges 
eau-glycérine. La diminution de la viscosité prouve 
l’abaissement du degré d’association, mais le choix 
d’un solvant polaire est critiquable. La seule modi- 
fication de fréquence observée par Saxena concerne 
la raie 674 cm1: la dilution s’accompagnerait 
d’une élévation de fréquence de 15 cm? Sundara 
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Rao [2] a confirmé ce changement et l’a obtenu en 
provoquant la dépolymérisation par élévation de la 
température de 30 °C à 100 °C. 

Par comparaison des spectres obtenus à la tem- 
pérature ambiante et à haute température, nous 
n'avons décelé aucune variation de fréquence 
consécutive à l’élévation de température. 

Le tableau suivant donne la valeur moyenne des 
fréquences pour les deux domaines de température 
avec l’écart moyen correspondant. Pour qu’une 
variation de fréquence soit certaine, il faudrait 
qu’elle ait une amplitude supérieure au double de 
l'erreur commise sur une mesure (5 cm1), ce qui 
n’est pas le cas, en particulier pour la raie de fré- 
quence 671 em 


ÉcarT Écarr 


Va50 V150-1650 V150-165° 
MOYEN MOYEN 
324,4 7 YA) D — 0,2 
413,6 2,6 412,5 3,6 — 1,1 
486,6 2) 485,8 2,0 — 0,8 
590,0 275 545,7 4,2 — 4,3 
671,0 on 671,2 39 + 0,2 
822,5 1,6 818,5 3,8 — 4,0 
848,4 229 849,2 mil + 0,8 
918,0 4,3 919,8 2,3 —+ 1,8 
977,4 29 OP 7,4 + 1,9 
1 065,8 f./# 1 062,4 2,8 — 3,4 
4 087,2 2,4 1 087,0 0,8 — 0,2 
1 105,8 3,6 1 106,1 kA + 0,3 
4 260,3 37 1 260,0 — — 0,3 
4 464,9 1,8 1 465,8 2,4 + 0,9 


S'il y a variation des fréquences Raman avec la 
température, celle-ci est trop faible pour être 


décelée dans nos mesures et n’est certainement pas 
de 15 em 1. 


Infiuence de la température sur la largeur des 
raies. — Les raies Raman de la glycérine sont 
larges et floues. Les molécules étant lourdes et 
associées par des liaisons, on ne peut attribuer la 
largeur. des raies aux rotations, mais les liaisons 
hydrogène qui unissent entre elles les molécules, 
perturbent les vibrations des liaisons voisines et 
élargissent les raies correspondantes. 

Un des effets de la dépolymérisation devrait être 
dans ces conditions d’accroitre la finesse des raies : 
pour étudier ce phénomène, nous avons mesuré sur 
les enregistrements les largeurs à mi-hauteur des 
pics présentant un profil simple aussi bien à tempé- 
rature ordinaire qu’à température élevée. Nous 
avons observé, dans un cas au moins, une nette 
diminution de largeur consécutive à l'élévation de 
température. 
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Nous avons d’abord vérifié que dans le domaine “1 
des faibles densités utilisées, la largeur d’une raie F 
ne dépendait pas de sa densité photographique. En M 
effet pour chaque température nous avons été 
amené à chercher les meilleurs conditions expéri- \ 
mentales, done à varier les temps de pose ; d'autre 
part à cause de la fluorescence, les contrastes sont 


variables d’une plaque à l’autre. Nous avons évalué 
la densité photographique relative de la raie de fré- 
quence 4 465 em 1 par rapport au fond continu de 


fluorescence. Pour sept densités croissant de 0,08 à « 
0,28, nous obtenions les largeurs respectives sui-. 
vantes (à 25 0C) : 9,5 — 9,5 — 11,0 — 9,5 — 9,0 =" 


9,5 — 10,0 mm. 


4 


+ 


Dans ces conditions il nous était possible de 
comparer les largeurs de raies obtenues à différentes 


températures. 


La précision sur la détermination des largeurs à « 


mi-hauteur est médiocre, au point d’enlever toute 


valeur aux résultats, si les raies ont un profil com-# 
plexe et changeant : cette complexité et cette insta- " 


bilité proviennent souvent du gram de la plaque, 
mais il est possible que l’élévation de température 
s'accompagne d’une variation de structure des 
raies comme si la dépolymérisation libérait des 
composantes nouvelles. Pour ces raisons, les me- 


sures de largeur n’ont été faites que pour la-raie de : 


fréquence 1 465 cmt et la bande 1 066-1 106 cm1 
qui s'inscrivent de façon satisfaisante aussi bien à 
température ordinaire qu’à température élevée. 


INFLUENCE DE LA TEMPÉRATURE SUR LA LAR- 


GEUR DE LA RAIE A 1 465 em. — Largeur obtenue 
pour divers clichés à 25 0C : 31,0 — 31,0 — 31,0 — 
32,6 — 31,0 — 29,4 — 35,9 cmt. 

Valeur moyenne : 31,7 cm1; moyenne des 
écarts par rapport à cette valeur : 1,5 cmt. 

A température élevée les résultats sont les sui- 
vants : 

Température : 140 
170 oC. 

Largeur : 19,6 — 24,8 — 23,8 — 26,4 — 24,5 — 
21,2 cm 

Valeur moyenne : 23,4 cm1; moyenne des 
écarts par rapport à cette valeur : 2,0 em. 

La diminution de largeur est de 8 em-1 alors 
que la valeur moyenne des écarts est de 2 em-t, 


150 — 155 — 155 — 160 — 


INFLUENCE DE LA TEMPÉRATURE SUR LA LAR- 
GEUR DE LA BANDE A 1 066-1 106 cm1. — Largeur 
obtenue pour divers clichés pris à 25 0C : 99,2 — 
97,8 — 89,2 — 103,2 — 111,6 — 105,1 — 102,6 — 
TUE: 

Valeur moyenne : 102,0 em-1; moyenne des 
écarts par rapport à cette valeur : 3,8 emr1. 

À température élevée les résultats sont les sui- 
vants : 

Température : 140 
160 — 160 — 165 oC. 


150 — 155 — 155 — 160 — 


N9ù6 


Largeur : 102,6 — 102,6 — 99,2 — 95,2 —. 92,6 
— 92,6 — 85,9 — 97,2 cm1. 

Valeur moyenne : 96,0 cm1; moyenne des 
écarts par rapport à cette valeur : 4,4 em. 


x 
x 
o x 
ED DÉTANEXOX LX x 
18 21 24 27 30 83 36 cm” 
Raiïe à 1465 cm”! 
x XAX XX XX x 
o H. Er (al 
| oO 
87 90 38 96 99 102 105 108 11 cm°! 


Bande à 1066-1105 cm°! 


x Largeur A 252 G © largeur à température élevée 


Free 
Raie à 1465 cm! 


CET e 
a 1Æ0°C 


c——— 


FEG. 2; 


La diminution de largeur est de 6 cmt, mais ce 
résultat est entaché d’une grande incertitude car 
cette diminution n’atteint pas le double de lécart 
moyen adopté comme valeur indicative de lerreur 
sur les mesures de largeur. 
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La figure 1 rassemble les résultats précédents et 
montre la répartition des nombres obtenus pour 
les largeurs de ces deux bandes en fonction de la 
température. On constate que pour la raie à 
1 465 cmt il y a deux groupes de valeurs nette- 
ment séparés. 

La figure 2 est obtenue par superposition de deux 
enregistrements microphotométriques de la raie à 
4 465 cm1 dans les deux domaines de température 
envisagés. La diminution de largeur, symétrique, 
est bien visible. 


Conclusion. — La raièé Raman à 1 465 cmt 
est due à la déformation plane de l’angle des 
liaisons C-H. Comme les liaisons hydrogène sont 
variables d’une molécule à l’autre et même aux 
deux bouts de la molécule, elles ont pour effet soit 
d'augmenter, soit de diminuer l’angle d’où une va- 
riation de la fréquence correspondante de la vibra- 
tion, de part et d’autre de la valeur que donnerait 
le groupement CH, isolé. Il en résulte un élargis- 
sement à peu près symétrique. 

Quand la glycérine est chauffée, les liaisons 
hydrogène sont détruites, au moins en partie, 
comme le prouve la diminution de viscosité et l’on 
doit se rapprocher alors de l’état vibratoire du radi- 
cal non perturbé : la raie devient plus fine, ainsi 
que nous l’avons trouvé pour la raie à 1 465 cm1. 

Le phénomène devrait être le même à l’autre 
bout de la liaison hydrogène, c’est-à-dire pour les 
raies de déformation de l’angle COH. Or la bande 
1 066-1 106 cmt est attribuée à la déformation 
plane de cet angle, mais cette bande est très large 
par suite de l’existence de deux sortes de fonction 
alcool dans la molécule. L’action de l'hydrogène 
sur l’oxygène est certainement plus faible que 
l’action inverse ce qui laisse prévoir une faible dimi- 
nution de largeur par dépolymérisation. Nos me- 
sures semblent confirmer ce résultat. 

A la température de 150 °C on obtient donc déjà 
une dépolymérisation importante de la glycérine. 


Manuscrit reçu le 127 septembre 1960. 
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SUR LES INTERFÉRENCES DE DEUX IMPURETÉS MÉTALLIQUES 


Par B. CAROLI, 


Physique des Solides, Faculté des Sciences, Orsay (Seine-et-Oise). 


Résumé. — On étudie comment les densités électroniques d’un métal sont perturbées par la 
présence de deux impuretés. Sur un modèle schématique, électrons libres à une dimension diffusés 
par deux puits de potentiel à, on analyse les oscillations de charge produites autour de deux impu- 
retés. On montre qu’elles comprennent des termes d’interférences, en plus des oscillations de 
charge produites par chacun des puits pris séparément. Ces termes correctifs sont d’un ordre 
supérieur en perturbation aux oscillations produites par deux puits isolés. 


Abstract. — The perturbations in the electronic density of a metal by two impurity atoms 
are studied. On a schematic model using one dimensional free electrons scattered by two à 
potential wells, the charge oscillations produced by the two impurities are analyzed. They are 
shown to be the sum of the charge oscillations due to each of the wells scattering independently, 
plus interferences terms. These corrective terms are of higher order in perturbation than the oscil- 


lations due to independent wells. 


1. Introduction. — La diffusion des électrons 
métalliques par des atomes d’impureté pris isolé- 
ment est assez bien comprise [1]. L’atome d’impu- 
reté perturbe localement la structure électronique 
et s’entoure de régions concentriques alternati- 
vement d’excès et de défaut de charge électro- 
nique [2]. Les oscillations de charge ont des lon- 
gueurs d’onde comparables à celles des électrons au 
niveau de Fermi et des amplitudes reliées au pou- 
voir diffusant des impuretés et décroissant avec la 
distance à l’impureté. Elles permettent, en parti- 
culier, d'interpréter de façon satisfaisante l’action 
des impuretés sur la résonance magnétique des 
noyaux de la matrice, par un calcul qui donne 
également une résistivité résiduelle convenable 
pour les impuretés [3]. 

Les interactions entre atomes d’impuretés dis- 
soutes dans les métaux sont beaucoup moins bien 
comprises. On sait que, dans l’approximation de 
Born, la résistivité résiduelle doit varier de façon 
parabolique avec la concentration [4] mais on 
ignore l’importance des termes d’ordre supérieur. 
L'énergie de liaison et la résistivité résiduelle d’une 
bilacune ont été estimées de façon approchée [5]. 
Mais en général les estimations d'énergie chimique 
d'interaction de deux atomes d’impureté, faites 
dans le modèle de Thomas-Fermi, négligent les 
oscillations de charge à grande distances qui jouent 
probablement un rôle prépondérant ; elles n’ont 
donc sans doute aucune valeur [6]. Dans l’inter- 
prétation des couplages magnétiques, où les oscil- 
lations de charge interviennent [7], on a supposé 
que les oscillations d’impuretés voisines super- 
posent leurs densités de charge sans les modifier. 

L'analyse des interactions de deux atomes d’im- 
pureté dissouts dans un métal présente donc un 
intérêt assez général. Dans cette étude préliminaire, 


nous avons simplifié le problème au maximum, en 
considérant un gaz d’électrons libres à une dimen- 
sion, diffusé par deux puits de potentiels à iden- 
tiques et de faible amplitude. Après avoir déter- 
miné les fonctions d’ondes et la densité d’états 
perturbés (paragraphe 2), nous étudiens comment 
les densités électroniques varient avec la distance 
entre les puits en considérant la région entre les 
impuretés puis à l’extérieur (paragraphes 3 et 4). 


2. Détermination des fonetions d’onde et de la 
densité d’états. — Nous considérons un gaz d’élec- 
trons à une dimension diffusés par deux puits carrés 


Ê 


ETC 


de potentiel, de profondeur V, et de largeur a 


(Ag. 1). 


2.1. DÉTERMINATION DES FONCTIONS D'ONDE. — 
En prenant l’origine des abscisses au milieu entre 
les deux puits, les fonctions d’onde d’énergie 
E = k?/2 s’écrivent : 


1 (a) — 2as sin kx. 


d 
AS E 
a (x) — 2448 


[(k/ko) cos kd sin ko(x — d) + sin kd cos ko(x — d)] 


+ 


r 


” 
PT, 


nv 


CRE T UA 


he Got et és 


No 6 : 


A 
LES (x) — 2aas 


[(k JA) cos kd sin ça + sin kdcos kg a] M À& — L) 


sin (a + d— I.) 
ktg ka + k, tg kd 


Ka + d— Lilas = (ko) re Lu) 
et 
LS (x) — 2ag cos kt. 
Jô(x) — 2as 
[cos kd cos kç(x — d) — (k]ks) sin kd sin kç(x — d)] 
Yslx) — 2as 


Ne Elan dan par on Me D) | 
[cos Ad cos kç a — (kJk,) sin kd sin koa] = ro Pr à 


ko — ktg kdtgk,a 
ki tg Koa + ktg kd° 


(2.1.2) 


tg kla + d— L)|s — — (k/ko) 


Les indices S et AS se rapportent respectivement 
aux fonctions symétriques et antisymétriques, les 
indices 1, 2, 3, aux trois régions I, II, III, définies 
(fig. 1), as et aus sont des coefficients de norma- 
hsation, L est une limite Paire où les fonctions 
d’onde s’annulent. 


2.2. DÉTERMINATION DE LA DENSITÉ D’'ÉTATS. 
THÉORÈME DE VON LAUE. — Introduisons un dé- 
phasage à(k) défini par : 


Km L = mr —(k). 


La densité électronique, avec deux électrons par 
état k s'écrit : 


nu 


2 2 dù 
= ! 
Km +3 — Xm 2 (£ nl 


La quantité dô/dk se calcule à l’aide des for- 
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mules (2.1.1) et (2.1.2) suivant que l’on considère 
le cas AS ou S, en posant : 


tg ka + d— Li = (kfko) (Ni(R) Di) i = 5 ou AS 


et en tenant compte de la relation de définition 
de (k),il vient alors 


m(k) = (2/7) (1/Ail?) (2.2.2) 


où À; est l’amplitude des fonctions (x) et se 
met sous la forme 


K& D?(k) + k? N?(k) 
(1 + te? kd) (1 + te? ko a) 


Lorsque L tend vers l'infini, la densité électro- 
nique à grande distance des impuretés est en 
moyenne (1/2) #|Aif? — (2/r). Elle est donc la 
même qu'en l’absence de perturbations : on re- 
trouve le théorème de von Laue [8] dans le cas 
particulier présent. 


PRET ET (2.2.3) 


3. Modification de la densité électronique pour 
deux puits 8. — Nous simplifions les formules en 
considérant désormais des puits à : 


lim A2a— y = Cte. 
ko 00 


a—0 
La formule (2.2.2), compte tenu de (2.2.3), 
devient : 
1 Ÿ Fe lee 
nas(K) = ea tte [1 EN .sin 24d + -.sim ka] 


—1 


Ar [1 + À sin 2kd + L cost ka] 


27 aë k 
L'introduction des puits perturbateurs doit 
changer localement la densité électronique. 


3.1. ÉTUDE ENTRE LES PUITS. — La variation de 
densité électronique par rapport au gaz non per- 
turbé s'écrit : 


A PU) fx 2Y? k? — [2yk cos 2kd — ÿ? sin 2kd]? — [4yk® sin 2kd + 24°? k? cos 2kd] cos 2kx 
1 2 x Jo 


où km est le vecteur d’onde du niveau de Fermi. 
Cette formule contient des termes d’addition simple 
des effets des deux puits sans interaction situés 
Pun à x = + d l’autre à x — — d et les termes 
d'interférence. Ceux-ci sont alors donnés par : 


A1 Pint.(z) = A3 p(x) — [A polx + d) + A: polz — dj]. 


où Ap,(x) est la contribution d’un puits isolé à 


AR Æ [2yk cos 2kd — y° sin 2kd]° ARE AOMIERS 
il vient : 
À; Fint.(x) (92182) 
9 k dk 
ions ro 9kd — y sin 2kd)?| === 
F-Jdr [2% (2k cos 2kd — y sin 2kd) 15 
k 
NS [k? — (2k cos 2kd — y sin 2kd)?] 
T 0 
dA 
1 re ele dima 9 }: 
(sin 2k(x + d) — sin 2k(x — d)) GE + v) EU 


l’origine. Ap,(x) s’obtient évidemment à partir de 
Ap(x) pour d = 0 et en divisant y par deux. 
En posant : 


E(k) = 4k4 + [2yk cos 2kd — y? sin 2kd]? 


— (2k cos 2kd — y sin 2kd)?] 


9 k 
—<w/f# k? 
T Te 


d) + cos 2k(x — d)) 


dk 
(cos 2k(x + (x)° 


GR + Y)E 
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_ 8.2. ÉTUDE À L'EXTÉRIEUR DES PUITS. — Les 
fonctions d’onde s’écrivent : 


VS — Aas(k) sin (kr + Das) ASP NE 


JO — Aslk) cos (kr + 8) ds = RS — RL. 
Ôs et das sont liés à la présence des impuretés par 


les relations (2.1.1) et (2.1.2) : 


fl 
te(kd | OS) notes (y/k) te kd 
le) 
et tg(kd + Ss) = (vlk) + tg kd. 


Un calcul analogue à celui du paragraphe (3.1) 
donne pour les termes d’interférence : 


ke, [fm k 
A3 Pint.(æ) = ESS À 


re 


(sin 2k(x + d) + sin 2k{(x — d)) 


— cos 2k(x 


EU) + d) dk 


(2k cos 2kd — y sin 2kd)°] k 


dk 
(4k? + y?) (ER) 


Vera j 
ee Cf (8 sin 2k{x + d) — 2*sin 2k{z — d) 


k dk 


— sin 2k(x — 3d) —— EX) 


(2k cos 2kd — + sin 2kd)?] 


M 
A fE 


(cos 2k(x + GË 


d) + cos 2k(x d)) (&k2 + y?) E(k) 


4 ke 
—T f à (cos ST ou 
TK Jo 2 


dk 
% sd) _ 


Cette formule et (3.1.2) sont exactes pour des 
puits à. 


à 
3 COS 2k(x 


4. Cas des faibles perturbations. — Nous nous 
contenterons 1c1 d’évaluer le premier ordre signi- 
ficatif pour de petites perturbations y. 


2 fk s Lkd 
| M COS 4Kc Be. A 
Ene k2 + ve ME 0a ir) 


Y? —2V2ya :y° d 
= ——<—e Le A (5 — Sin kr a) 


? TT 


A ATEN TE 


A 


y: [* cos #km d d 


TH km M 
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2 ky cos 2k(x + d) 
Aa finir = À [EE en dE + Oalrl 
y? —Veræ+d 2y° s 
Rene Fee (x + d) 


(5 — SRE R a)) 
y? cos 2km'x + d) 
T km 4 


Alors que pour un puits isolé, les deux premiers 
termes significatifs sont : 


Y L fox k sin 2kx dé 
k2 + 2]h 


-f£m cos 2kc 
.dA 
[x EST 0e 


Apolt 


en ET eye +1 Sin 2km #| 
D 210 K72 
y: 
—_—— EC -Ÿe 
k e 


y? fcos 2km x : —?)] 
2 [AE + 2x | Sin 24m x 5)|- 

Sin x est mis pour sinus intégral de x. Une inté- 
gration par partie donne : 


2 sin 2km x 
A3 Pint.(x) TA 


250 Ce LE 
et 
cos 2km x 2sin 2km x - 
Apolt) — + = 
27 kmax ET Kÿx 


SIT CO: 


5. Conelusion. — Dans le modèle très théorique 
étudié, les variations de densité électroniques, dues 
aux interférences entre impuretés, ont donc les 
propriétés suivantes : 1. Elles sont d’un ordre supé- 
rieur aux variations de densité des puits isolés. 
2. Au premier ordre significatif, elles sont cons- 
tantes entre les deux puits et décroissent à grande 
distance comme une fonction oscillante de 2km x 
sur la distance x. 

On peut donc s'attendre, dans les alliages de 
concentration finie, à des couplages non nuls par 
interférence. Les résultats à trois dimension 
doivent cependant être un peu différents : Ceux-ci 
ne seraient plus constants entre les deux puits. 

Je tiens à exprimer toute ma reconnaissance à 
M. le P° J. Friedel, qui m’a inspiré ce travail et 
m'a guidé tout au long de sa réalisation. 


Manuscrit reçu le 15 février. 
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LE FILTRAGE MATHÉMATIQUE DANS LA SPECTROSCOPIE PAR TRANSFORMATION DE FOURIER 


Par J. CONNES (*) et V. NOZAL,, 


Laboratoire Aimé-Cotton, Bellevue. 


Résumé. — On est conduit pour obtenir le rapport signal/bruit maximum dans le spectre calculé 
numériquement à partir d’un interférogramme à utiliser un nombre de points très supérieur au 
nombre d'éléments spectraux, d’où un temps de calcul très élevé. Dans la méthode du filtrage 
mathématique on calcule d’abord la convolution de l’interférogramme enregistré par la réponse 
percussionnelle du filtre idéal isolant les fréquences Fourier. Il suffit ensuite de faire la transformée 
de Fourier du nouvel interférogramme avec un nombre de points minimum d’où une réduction 
considérable du temps de caleul total. 


Abstract. — In order to obtain the maximum signal to noise ratio when computing a spectrum 
from an interferogram, the number of points that have to be measured is much greater that the num- 
ber of spectral elements ; hence the computing time is excessive. With the method of ‘‘ mathe- 
matical filtering *” one first compute the convolution of the interferogram by the impulse 
response of an ideal filter isolating the wanted domain of Fourier frequencies. Then one computes 
the Fourier transform of the new interferogram with a smaller number of points ; the computing 
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time is thus much reduced. 


Une étude précédente de la méthode de spectros- 
scopie par transformation de Fourier a montré 
dans quelles conditions le spectre pouvait être cal- 
culé en faisant la transformée de Fourier de l’inter- 
férogramme à partir de valeurs discrètes relevées 
sur cet interférogramme [1]. En l’absence de bruit, 
il suffit pour reconstituer le spectre de relever un 
nombre de points égal au nombre d’éléments spec- 
traux à étudier. Mais on montre que pratiquement 
le spectre du bruit étant plus étendu que celui des 
fréquences de Fourier, pour obtenir un rapport 
signal/bruit dans le spectre calculé voisin du rap- 
port signal/bruit maximum qu’on obtiendrait en 
calculant le spectre par une intégrale de Fourier, le 
nombre de points à utiliser est beaucoup plus 
grand. La durée et le prix de revient du calcul d’une 
transformée de Fourier étant proportionnels au 
nombre de points d’entrée sont augmentés d'autant. 
On a donc intérêt à limiter le spectre du bruit au 
domaine du spectre à étudier. Plusieurs méthodes 
d'enregistrement de l’interférogramme sont alors 
possibles pour obtenir ce résultat : emploi d’un 
filtre passe bande, méthode du changement de fré- 
quence de Mertz [2]. Mais elles présentent toutes 
une certaine complication instrumentale. Par exem- 
ple, un filtre passe bande avec une bande passante 
centrée sur la valeur moyenne des fréquences à 
étudier est difficile à réaliser aux fréquences bas- 
ses généralement utilisées. 

La méthode de filtrage mathématique utilise 
l’interférogramme Z/'() enregistré à travers un 
simple filtre passe bas. À partir de valeurs discrètes 
de Z'(5) on calcule un second interférogramme /”(à) 
en effectuant la convolution de Z(3) par la réponse 


(*) Adresse actuelle : Service de Calcul Numérique. 
Observatoire de Meudon. 


percussionnelle d’un filtre dont la bande passante 
coïncide avec le domaine des fréquences de Fourier 
à étudier. Il suffit de faire ensuite la transformée de 
Fourier de Z"(5) avec le nombre de points minimum 
pour obtenir dans le spectre un rapport signal/bruit 
maximum. Comme le temps de calcul du produit 
de composition est faible par rapport à celui de la 
transformée de Fourier, le temps de calcul est voi- 
sin du temps de calcul minimum. 


I. Principe de la méthode. -— Considérons un 
interférogramme Z'(3) enregistré à travers un 
simple filtre passe-bas. C’est la superposition de 
l’interférogramme (5) qu’on aurait en l’absence 
de bruit et d’un bruit æ(à) (fig. 1). Le spectre corres- 


{7 


Fra. 1. — Interférogramme enregistré 1’(à). 


pondant B"p(o) obtenu par calcul de l'intégrale de 
Fourier donc avec le rapport s/b maximum compa- 
tible avec le temps d’enregistrement 7 de l’inter- 
férogramme, est composé du spectre optique à 
étudier Bp(o) limité au domaine Ac, compris entre 
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61 et 6, auquel se superpose un bruit X(6) (fig. 2). 
La figure 3 indique la position relative du spectre 


B') (c) 


6: ( 
NCA 
He — 


Fic. 2. — Spectre correspondant à Z’(à). 


Fi&. 3. — Position relative du spectre du bruit 


et du spectre optique. 


optique Bp(o) et du spectre de puissance du bruit 
contenu dans Z/'(5) Dans une transformée de 
Fourier numérique la fonction d'appareil est cons- 
tituée d’une série de pics semblables distants de 1 /h 
si À est la distance entre 2 points relevés sur l’inter- 
férogramme (fig. 1). Si on choisit » le plus grand 
possible pour reconstituer le spectre sans ambi- 
guité : À — 1/26, (1), un seul pic de la fonction 
d'appareil explore le spectre à étudier, mais plu- 
sieurs pics de la fonction d'appareil explorent le 
spectre du bruit, d’où une augmentation du bruit 


Fic. 4. — Fonction d'appareil 
explorant le spectre optique et le spectre du bruit. 


dans le spectre calculé (fig. 4). Dans ce cas, avec 
apodisation (2), le nombre de points utilisés pour 


(1) Pour choisir le pas k convenant à chaque cas, il faut 
tenir compte de la largeur du spectre optique et de sa 


position par rapport à l’origine. Si o, —_kAo,, on a 
h = 1/26, ; dans les autres cas k est légèrement infé- 
rieur [1]. 


(?) La fonction d’appareil est une fonction en sin x/« 
qui présente des pieds qui peuvent être gênants. Pour 
l’apodiser il faut pondérer l’interférogramme par une fonc- 
tion A(Ô). Si c’est une fonction en triangle, pour obtenir la 
même résolution que sans apodisation, il faut employer un 
interférogramme 2 fois plus long et le nombre de points à 
relever est double. 
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faire la transformée de Fourier est n — 
M étant le nombre d'éléments spectraux à étudier. 
Ce que nous nous proposons, c’est d'obtenir à par- 
tir de Z'(3) un interférogramme 7/"(5) dont le 
spectre calculé par une intégrale de Fourier serait 


B’p(o) — B'p(6).G(o) 


FiG. 5. — Transformée de Fourier de 7”(à). 


G(o) étant une fonction créneau qui vaut 1 dans 


les domaines Ac, symétriques par rapport à l’ori- 
gine et est nulle ailleurs (fig. 6). Le bruit qu'il 


Fic. 6. — Filtre idéal. 


contient est donc bien nul pour toutes les fré- 
quences extérieures au domaine 6:, 62. Si, momen- 
tanément, on ne s'intéresse qu'aux fréquences com- 


prises entre 6, et 6,, le gain du filtre qu’on choisira 


aura la forme représentée par la figure 7a et le 


Rice FIG; 


spectre Sp(o) sera nul en dehors des fréquences 61, 
62 (fig. 7b). 


F1G. 8. — Réponse percussionnelle du filtre idéal. 


Soit F(5) la réponse percussionnelle du filtre dont 
la courbe de gain est G(o). C’est par définition la 
transformée de Fourier de G(o). 


sin tAocÈ 


A ere 


COS 2T70m Ô (2) 


Om — (01 + 62)/2 étant la fréquence moyenne du 
spectre à étudier (fig. 8). 


N°6 à 
2M, @ 


(fig. 5) (1) 
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De l’équation (1) on déduit immédiatement la 
relation intéressante entre Z”(5) et Z'(5). 


1"(S) = I'(S) A(S) + F(5) (3) 


A() étant la fonction d’apodisation utilisée pour 
calculer B’p(o) à partir de Z'(3) (éventuellement 
égale à l’unité s il n’y à pas d’apodisation). Il suffit 
donc de faire la convolution de J'(3) avec la 

réponse percussionnelle d’un filtre passe-bande 
centré sur le domaine à étudier pour obtenir un 
interférogramme dans lequel le spectre du bruit 
est nul en dehors des fréquences intéressantes. 

Pratiquement pour calculer Z”(5) on utilise une 
réponse percussionnelle Æ”(3) limitée à 3 — L/ soit 


F'(5) = F(5).D(S) … (4) 


D(S) étant une fonction créneau valant 1 pour 
[3] < L'et 0 ailleurs. 

Les courbes de gain G’ des filtres pratiquement 
réalisés ont pour expression : 


G'{s) = Go) x% d(os) avec d(o) = T[D(Ô)] (fig. 8a et 8b). 


—v20. 


 S © 
Hu "ei 
s“F 


Een 
o 


d(s) 


F1c. 8a. — Filtres non apodisés. F1. 8b. — d(o). 


Elles présentent donc des oscillations (phéno- 
mènes de Gibbs) et s’éloignent d’autant plus des 
courbes théoriques que le quotient ! — Ac,/1/2L 
est plus petit. Pour amortir ces oscillations, 1l faut 

, pondérer la réponse percussionnelle F'(5) par une 
fonction convenable, par exemple 


A) = (8/LYT. 
Alors : 
G'(o) — Go) %x als) avec a(o) — T[A'(o)] (fig. Ia et 9b). 


Pour les fréquences limites le filtre introduit une 
altération de l’amplitude. Si on appelle Ao la zone 
utilisable du filtre pour laquelle l’altération est, par 
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Fic. 9a. — Filtres apodisés. Fi. 96. — ao). 


exemple, inférieure à 1/100 (fig. 10a), on peut 
définir un facteur de qualité du filtre par : 


p = Ac/Ao4. 


Fic. 10a. 


Par exemple avec la fonction de pondération 
A'(3) envisagée : 


__ L'Aco — 2,2 
L’A6, 


p = (fig. 10b). 


012345 10 L'AG, 


F1ic. 10. 


II. Réalisation. — Pour réaliser avec un ordi- 
nateur la convolution : 


I'(S) — I'(5).A(S) — F(3).A‘(à) (5) 


on ne peut utiliser que des valeurs discrètes de 
l’interférogramme et de la réponse percussionnelle. 


24 
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Nous allons montrer que lorsqu'on fait avec des 
valeurs discrètes équidistantes la convolution de 
2 fonctions à spectre limité on obtient bien des 
valeurs discrètes du produit de composition ; mais 
si les fonctions ont un spectre illimité, on ne peut 
obtenir par cette méthode que des valeurs appro- 
chées du produit de composition, l’approximation 
étant d’autant meilleure que le pas k' est choisi 
plus petit. 


19 CONVOLUTION DE 2 FONCTIONS A SPECTRE 
LIMITÉ. —- Soient 2 fonctions Æ(3) et X(3) ayant 
des spectres limités e(o) et k(o) comme l’indiquent 
les figures 11a et 116. Leur convolution s’écrit : 


P(S) — E(S) x ÆK(ÿ). (6) 
Comparons P(3) et la convolution : 
P'(S) = h'IE(S).Ry (8) x K(S).Rr(8)] (7) 


R;,(3) étant une distribution de Dirac à support 
périodique de pas k’. La transformée de Fourier 
de ce produit de composition s’écrit : 


p'lo) =%[ls) * Rawlo)].[A(s) % Ra (oi] 


\ 


(His 14014 le): 


a) 


qe(c)* Ray, (6) 
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C’est une fonction périodique de pas k'. Suppo- 
sons que »' ait été choisi de telle sorte que les 
zones qui se recouvrent de la fonction 


e(s) x R,Jh'(0) 


ne pénètrent pas dans le domaine occupé par 
k(s) x R;fh'(o), ce qui correspond dans notre 
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cas à : 1/h! — (63 — 64) /2. On voitimmédiatement 
que dans ce cas : 


: 
po) =; lels) * k(o)].Rimio). 
On en déduit que sa transformée de Fourier a 
pour expression : 
P'(S) — P(S).Ry (6). (8) 


Les nombres calculés P/(3) sont donc bien une 
série de valeurs discrètes équidistantes du produit 
de composition cherché P(). 


29 CONVOLUTION DE DEUX FONCTIONS A SPEC- 


TRE ILLIMITÉ. — Quand les deux fonctions ont des 
spectres illimitées les relations entre P(Ô) et P'(à) 
sont beaucoup plus complexes. Nous examinerons 
en détail uniquement un cas particulier, celui de 
la convolution d’un interférogramme Z'(à) enregis- 
tré à travers un simple filtre RC, de constante de 
temps tr — À C; et de la fonction F'(Ô) définie pré- 
cédemment, puisque c’est celui qui nous intéresse. 
La fonc tion calculée a alors la forme : 


P'(5) = l’(5).Ry(S). (8 bis) 


Elle ne prend que des valeurs discrètes et a 
évidemment un spectre périodique. L’interféro- 
gramme continu /”(3) qui aurait les valeurs précé- 
demment calculées pour des points d’abscisse 
Ô — ph'a lui un spectre non périodique. 

Sans nous préoccuper de l’expression exacte de 
1"(3) ainsi défini nous allons préciser le spectre du 
bruit contenu dans /”(3) et dans sa transformée de 
Fourier. 

La courbe de gain correspondant à la réponse 
percussionnelle F’(5).R;.(3) est une fonction pério- 
dique de période Z/h' (fig. 12). Le spectre du bruit 


G (S)*R 1/7 
1 


Va —— à 


Fic. 12. — Filtre périodique. 


contenu dans l’interférogramme continu Z'(3) a 
alors pour expression : 


Boo) — Boo) (IG) Riym(o)] (fig. 13) (9) 


et l'écart moyen quadratique des fluctuations dans 
1"(5) vaut : 


—C9 


Il est évidemment inférieur à l’écart moyen qua- 
dratique des fluctuations dans (5) qui vaut 


L— 


+00 
O% Se Bp!5) do 


N°6 


N°6 


et sera d'autant plus faible que Z/h! sera plus 
grand donc que le pas k' utilisé pour faire la convo- 
lution sera plus petit. Étant donnée la méthode de 
calcul, le pas À utilisé pour faire la transformée de 
Fourier de Z"(3) qui ne dépend plus que de la 
largeur Ac du spectre optique à étudier et de sa 
position par rapport à l’origine est obligatoirement 
un multiple de k' : 


R= y. 


La distance séparant deux pics de la fonction 
d'appareil est alors : 


C’est un sous-multiple de 7 /h'. Rappelons que la 
puissance moyenne des fluctuations dans le spectre 
obtenu par transformation de Fourier de l’interfé- 
rogramme est égale à l’énergie découpée dans le 
spectre du bruit par le carré de la fonction d’appa- 
reil [1]. L'écart moyen quadratique des fluctuations 
dans le spectre calculé à partir de Z”(3) avec un 
pas À est rigoureusement le même que celui qu’on 
aurait obtenu avec un pas k' puisque les pics sup- 
plémentaires de la fonction d’appareil tombent en 
des endroits où le spectre du bruit est nul (fig. 13). 
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F1G. 13. — Fonction d'appareil de période 1/h. 
explorant le spectre du bruit contenu dans 7”(Ô). 


Le problème du calcul de l’écart moyen quadra- 
tique des fluctuations dans le spectre obtenu par 
transformée de Fourier numérique d’un interféro- 
gramme enregistré à travers un simple filtre RC a 
déjà été traité [1]. On montre que la variance dans 
le spectre calculé avec un pas À: 6 est liée à la 
variance 6% qu’on aurait dans le spectre calculé par 
une intégrale de Fourier par la relation : 


* avec 


_h  sinkhfr 
” 2rcos hhfr— cos 27h 


D,,(h[r) 


v, étant la fréquence particulière envisagée. 

Donc, dans notre cas, l’écart moyen quadratique 
des fluctuations dans le spectre final 6 sera lié à 
l’écart moyen quadratique minimum qu’on aurait 
eu en calculant directement par une intégrale la 
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transformée de Fourier de Z'(3) par la relation : 


ox? = 0ù D,,(h’[r) 


avec 
DIE h’ sin À hole 
BAR 15) - 2rcos h h’[r —cos 27 v, h’ DÉREenEt, 
| 1072) 
V2=0 
4 


0 5 


10 A 
Fic. 14. — Fonction ®,({k’[x). 


h' étant cette fois le pas choisi pour faire la convo- 
lution (1). 

Donc pour obtenir à partir de Z”(3) le même 
rapport s/b dans le spectre que lorsqu'on fait direc- 
tement la transformée de Z'(3), il faut relever: 
autant de valeurs sur Z'(5). Mais le nombre de 
valeurs de Z"(3) qu’on calculera et qu’on utilisera 
pour faire la transformée de Fourier sera nota- 
blement inférieur. Il sera déterminé uniquement 
par la largeur du spectre optique et sa position par 
rapport à l’origine. 

On peut évaluer la perte S en rapport s/b qui 
résulte du fait que k’ est fini en faisant le quotient 
du rapport s/b obtenu avec h' ayant une valeur 
déterminée à celui qu’on obtiendrait avec h' infi- 
niment petit : SC 


(fig. 14 bis). 


F1. 14 bis. — Variation durapports/b dans le spectre final 
en fonction du pas k’ avec lequel est faite la convolution. 


On voit d’après la figure 14 bis qu’il y a intérêt, 


(1) Les courbes ®,,(h/+r’) sont pratiquement indépen- 
dantes de vs. 
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comme dans le cas d’une transformée de Fourier 
directe, à choisir k’ dé l’ordre de +. 


Méthode de calcul. — Le calcul comprend 5 par- 
ties distinctes : 

19 Calcul des N' valeurs de la réponse percus- 
sionnelle F’(3) du filtre avec N' = L'Jh'. 

20 Pondération de ces V' valeurs de Æ”(à) par 
A'(S). 

30 Pondération de l’interférogramme ZJ'(à) par 
A(à). 

49 Calcul du produit de composition 

I"($) = T'(6).A(S) *x F”(8).4'(). 

5° Calcul de la transformée de Fourier en cosi- 
nus de Z”(5). 

La durée des 3 premières opérations est négli- 


geable devant celle des deux dernières (1). Le temps 


du calcul du produit de composition peut se mettre 
sous la forme : 


Ge —= kN'n secondes 


N' étant le nombre de valeurs discrètes de F'(à) 
utilisées, z le nombre de points de Z”(5) calculées 
et k un coefficient numérique qui dépend de la 
machine utilisée. Nous donnons sa valeur pour 
l’ordinateur 704 IBM utilisé dans un problème pour 
lequel le nombre de points d’entrée N (qui n’inter- 
vient pas directement dans le calcul de Ge, mais 
seulement par l'intermédiaire de n) est supérieur au 
nombre de mémoires rapides. Ceci oblige à faire 
entrer les valeurs N en plusieurs fois et augmente 
considérablement le temps de calcul. Dans ces con- 
ditions : 4 — 10%. Rappelons que le temps de 
calcul d’une transformée de Fourier classique à 
partir des N valeurs relevées sur l’interférogramme 
serait : 


Gr = 0,8.10—% NN, secondes 


N, étant le nombre de points de sortie ; on peut 
l’évaluer en fonction de n, nombre minimum de 
points utilisés pour faire la transformée de Fourier 
n est de l’ordre de 2M, M étant le nombre d’élé- 
ments spectraux. Si on calcule 3 points par élément 
spectral pour faciliter le tracé du spectre : 
N, = 2M—\{(3I2n, 
Dans ce cas : 
Gr = 1,2.10—5% Nn secondes 


ce qui conduit avec l’ordinateur 704 utilisé dans les 
conditions que nous avons précisées à : 


Gel Gr = 0,8 N’IN. 


(4) Les sinus et cosinus qui interviennent dans le calcul 
de F’(8) aussi bien que dans le calcul de la transformée 
de Fourier de Z”(ô) sont calculés par récurrence à l’aide 
des formules de Tchébicheff : 

cos (p + 1) x = 2 cos x cos px — cos (p —1) x 
sin (p + 1)x = 2 cos x sin px — cos (p — 1) x. 
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Examinons dans un cas particulier le gain de 
temps qu’on peut réaliser en faisant un filtrage 
mathématique. Supposons que le spectre à étudier 
s’étende sur un domaine 


(fig. 15) 


et que l’interférogramme soit enregistré à travers 
un simple filtre passe-bas. Ce filtre introduit des 
déphasages pour les diverses fréquences contenues 
dans le signal. Ces déphasages se traduisent dans le 
spectre reconstitué par une dissymétrie de la fonc- 
tion d'appareil [1]. Si on tolère une dissymétrie de 
1/100 entre les hauteurs de deux premiers pieds 
négatifs on peut choisir une constante de temps : 


AG = Om — Om = OM|S 


Tr = 1/2rom. 


Le spectre optique et le spectre du bruit ont 
alors la position indiquée dans la figure 15. Si on 


Fic. 15. — # = 1[2y hr = 7x S = 0,7 h = 5k. 


veut obtenir dans le spectre final un rapport s/b 
qui soit les 7/10 du rapport s/b maximum, il faut 
choisir »’ tel que [1]: 
hlr= 7" 
c’est à dire : 
k = 1/26. 


Supposons de plus que nous voulions traiter ce 
problème à la résolution 6 000 ; des relations : 


N=1/N et R= Low. 
On déduit que : 
N:="122000: 
L'élément spectral vaut : 
806 — 6/6 000 
et le spectre à étudier contient : 
M = Ao]5o = 12 000 éléments spectraux. 


Si on calcule 3 points par élément spectral, le 
calcul de la transformée de Fourier directe durerait: 


Cr — 0,8.10—5% X 12 000 X 3600 — 9 heures 36 minutes. 


Pour réaliser un filtrage numérique pratiquement 
di : op — 0,98, il faut d’après la figure 10b 
choisir 


Ac L'= 60 soit ZL’= 300/om et N’= L'[h' — 600. 


N°6 
Le temps de calcul du produit de composition 
correspondant est 
Ge —= 600 X 2 400 X 16—3s — 24 minutes. 


La nouvelle transformée de Fourier demande 
alors : 


Or — 0,810. %X 2 400:X°3 6008 = 1 h 55 minutes. 
* d’où une réduction du temps de calcul total par un 
facteur 4,15. . 


III. Vérifications expérimentales. — Elles por- 
tent sur deux points précis : calcul de l’écart 
moyen quadratique des fluctuations dans Z”(5) et 
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comparaison de deux spectres obtenus à partir 
d’un même interférogramme, l’un par transformée 
de Fourier directe avec un pas k/, l’autre après fil- 
trage numérique, le pas utilisé pour faire la convo- 
lution étant également }’. 


19 ÉGART MOYEN QUADRATIQUE DES FLUCTUA- 
TIONS DANS /”(3). — Soit un interférogramme Z”(t) 
constitué uniquement par du bruit (fig. 16a) (1). 

Constante de temps : + = 25. 

Spectre de puissance du bruit : 


Bp(v) = _. 


RTE TL (fig. 16b). 


Bp () 


b) 


Fic. 16b. 


B5 (0) 


Fic. 18a. — J,(t), Fic. 18b, Fic. 18c. 


Écart moyen quadratique des fluctuations : 
+ 
LT — fe Bo(v) dv = Bolé. 


Avec l'échelle arbitraire employée : 67 = 37 u. 


(1) Les variables sont ici t et v, proportionnelles à 8 et o 


puisqu'on a 8 = Vtet v — oV, V étant le double de la 
vitesse du miroir mobile de l’interféromètre à 2 ondes 
employé. 
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a) 1% filtrage. — Le filtre est centré sur l’origine ; 
largeur du filtre : Av, — 0,09 Hz (fig. 17b). 

Sa réponse percussionnelle est représentée par la 
figure 17c. 

Si hk' était nul, l’écart moyen quadratique des 
fluctuations dans Z(t) aurait pour valeur : 


An 

Fa dv 
LA. PES PA RPM ANRE 
OCg° — so / 1 + (anv)? 


ire +Av/2 
— B(0) [Arc tg an v| TA po. 


Avec les unités choisies on trouverait : 
ch oidi 68 ri èt-e, = 270. 


Le pas utilisé k’ étant de 1,25, le bruit dans Z;(r) 
est légèrement supérieur. D’après la courbe 14 brs 
pour k'/r — 06,5 =—10,98 donc 
6, devrait être : 


c% — 27.X 1,02 — 27,54. 


La valeur expérimentale est : 26 u (fig. 17a). 


b) 2€ filtrage. — Le filtre est centré sur l’origine ; 
largeur du filtre Av, — 0,0225 Hz (fig. 18b). Sa 
réponse percussionnelle est représentée sur la 
figure 18c. 

Écart moyen quadratique des fluctuations pré- 
vues avec h' — 1,28. 


VAE 


Valeur expérimentale : 15 u (fig. 18a). 

Les 2 interférogrammes /:(t) et Z;(t) paraissent 
lissés par rapport à /'(t). Leurs rayons de corré- 
lation sont de l’ordre de l’inverse du double des 
fréquences de coupure soit 6 s pour le premier et 
24 s pour le second. 


29 COMPARAISON D'UN SPECTRE OBTENU PAR 
TRANSFORMÉE DE FOURIER DIRECTE AVEC UN 
SPECTRE OBTENU APRÈS FILTRAGE MATHÉMATIQUE. 
— Soit un interférogramme synthétique J/({) obte- 
nu en ajoutant au bruit précédent /’(t) une cons- 
tante. Le spectre correspondant se compose d’une 
raie à la fréquence 0 à laquelle se superpose un 
bruit X(6). : 

Dans le spectre B;(v) obtenu par transformée de 
Fourier directe avec le pas k' — 1,2 s, l’écart 
moyen quadratique des fluctuations mesuré sur le 
spectre vaut 6, — 1 280 u (fig. 19a). 

La transformée de Fourier de l’interférogramme 
1:(t) obtenu avec le filtre de largeur 0,0225 Hz 
donne le spectre 19b. Elle a été effectuée avec le 
pas À maximum soit h — 18 h/. 

Le spectre B,(v) est évidemment périodique de 
période 0,045 Hz, donc il est différent dans son 
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ensemble de B,(v), mais les parties comprises entre 
v = 0et v — 0,0225 Hz sont rigoureusement iden- 
tiques et l’écart moyen quadratique des fluctua- 
tions o, vaut aussi 1 280 u. 


F16. 19b. — Spectre obtenu après filtrage numérique : 
REA NS EAST LE 


IV. Conelusion. — Pour obtenir, avec un filtrage 
mathématique, un rapport s/b donné dans le 
spectre, le pas k' qui est à utiliser pour faire la 
convolution est exactement le même que celui qui 
serait utilisé pour faire la transformée de Fourier 
directe. On doit donc choisir »’ de l’ordre de 7, 
constante de temps qui a servi à enregistrer l’inter- 
férogramme. Donc on relève autant de points sur 
l’interférogramme. Mais la transformée de Fourier 
est faite à partir du nombre de points minimum 
pour qu’il n’y ait pas recouvrement des spectres 
optiques, la valeur du pas À ne dépendant plus que 
de la largeur du spectre optique et de sa position 
par rapport à l’origine, d’où une réduction consi- 
dérable du temps de calcul (?. 

Les auteurs remercient le Comité Européen de 
Calcul Scientifique pour ses attributions d'heures 
de calcul à l’ordinateur 704 IBM. 


Manuscrit reçu le 14 janvier 1961. 


(*) La méthode du filtrage numérique a permis d’étudier 
dans un temps raisonnable, par transformation de Fourier, 
la lumière de recombinaison intrinsèque dans le germa- 
nium. 


BIBLIOGRAPHIE 


[1] Conxes (J.), Revue d’Optique, à paraître, février 1961. 


[2] Merrz (L. W.), Symposium de Stockholm, 1959. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


EFFET D'IMPURETÉS NON MAGNÉTIQUES DE SUBSTITUTION 
SUR LES COUPLAGES MAGNÉTIQUES DES MÉTAUX 
APPLICATION A LA RÉSISTIVITÉ ÉLECTRIQUE 


Par YA ROCHER (*), 


Résumé. — La dissolution d’impuretés non magnétiques dans un métal magnétique est accom- 
pagnée. de fortes variations des propriétés électriques et magnétiques. Pour tenir compte de ces 
effets, un modèle simple est étudié, dans lequel l’impureté couple magnétiquement deux atomes 
magnétiques diamétralement opposés par rapport à elle. Dans la limite des approximations de 
champ moléculaire et de spins classiques, le déplacement du point de transition, la forme des 
corrélations de spin à haute température, les modifications de structure magnétique autour d’une 
impureté à basse température, sont étudiés en fonction du signe et de la force de cette interaction 
indirecte. Les résultats obtenus sont appliqués au calcul de la résistivité électrique. Enfin, l'étude 
de résultats expérimentaux sur les alliages Mn Cu et Pu Al montre que ce modèle conviendrait 
particulièrement bien pour expliquer les propriétés de ces alliages. 


Abstract. — The presence of non magnetic impurities in a magnetic metal strongly affects its 
electrical and magnetic properties. In order to explain these effects, a simple model is studied, 
in which the spins of two magnetic atoms are coupled through an intermediate impurity. Within 
the limits of molecular field and classical spin approximations, the displacement of the transition 
temperature, the high temperature shape of spin correlations, the low temperature changes in 
magnetic structure around an impurity are studied and related to the sign and strength of the 
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indirect interaction. 


The results are used to calculate the electrical resistivity. The model 


seems to explain particularly well the properties of Mn Cu and Pu Al alloys. 


I. Introduction. — Des études récentes ont per- 
mis de relier certaines propriétés physiques des 
métaux magnétiques, la résistivité électrique [1], 
[2] et le pouvoir thermoélectrique par exemple, à la 
structure magnétique de ces métaux. Ainsi la résis- 
tivité par désordre de spin a-t-elle été mise en évi- 
dence ; dans certains cas, elle est bien supérieure à 
la résistivité due aux phonons. 

_ La substitution à un atome magnétique d’un 
atome non magnétique (Al, Cu, Zn, Ga, etc...) per- 
turbe localement la structure magnétique du métal 
et l’effet doit s’en faire sentir sur les propriétés 
magnétiques (températures de Curie ou de Néel, 
moment magnétique, susceptibilité, spectre de dif- 
fraction de neutrons) et sur les propriétés élec- 
triques. Expérimentalement, en effet, on observe un 
abaissement des points de Curie ou de Néel, une 
diminution du moment magnétique atomique, en 
général proportionnels à la concentration C d’im- 
puretés, mais avec un facteur de proportionnalité 
mférieur, égal ou supérieur à l'unité. Dans certains 
cas, on observe de très fortes résistivités résiduelles 
à 0 °K, plus ou moins proportionnelles à la concen- 
tration d’impurétés, et associées à de très grandes 
sections efficaces de diffusion par atome d’impureté. 

Ces faits suggèrent que l’impureté non magné- 
tique n’a pas un rôle purement passif ; s’il en était 
ainsi, en effet, les variations relatives du point de 
Curie et du moment magnétique seraient égales à la 


(*) Faculté des Sciences d'Orsay (Seine-et-Oise), France, 
Service de Physique des Solides, Groupe II, B. P. n° 11. 


concentration, et les sections efficaces de diffusion 
de l’ordre des sections atomiques. L’action des 
impuretés peut alors être recherchée dans deux 
directions : 1) les impuretés peuvent modifier di- 
rectement la structure électronique des atomes voi- 
sins dans la matrice ; 2) sans changer apprécia- 
blement cette structure, les impuretés peuvent sim- 
plement influer sur l’orientation des moments ma- 
gnétiques à leur voisinage, en changeant localement 
la nature des couplages magnétiques. 

Il est bien connu que, dans le nickel et le cobalt, 
la plupart des impuretés non magnétiques rédui- 
sent le moment magnétique dans des proportions 
plus fortes que s’il y avait simple dilution du métal 
ferromagnétique [3], [4]. Dans ces alliages, les impu- 
retés réduisent donc les moments magnétique des 
atomes voisins de la matrice ; un modèle du pre- 
mier type est donc vraisemblablement nécessaire. 

On sait, par contre, que beaucoup d’impuretés 
non magnétiques diluent le moment magnétique du 
fer en proportion de leur concentration, donc 
n’affectent guère la structure électronique de ce 
métal. Des considérations théoriques [5] suggèrent 
qu’il doit en être de même pour les alliages d’autres 
métaux de transition porteurs de moments perma- 
nents (Cr, Mn, …). Tel doit aussi être le cas des 
alliages des métaux de terres rares et peut-être de 
certains métaux lourds comme le plutonium, sur 
lequel nous reviendrons à la fin de cet article. Dans 
tous les cas, on peut penser un modèle du second 
type suffisant. 

C’est ce modèle simple que nous allons dévelop- 
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per. Nous montrerons qu’il permet d'expliquer les 
maigres résultats expérimentaux. Il prédit cer- 
taines relations entre les effets des impuretés sur 
diverses propriétés : ordre local, point de Curie, 
résistivité. La vérification de ces relations dans 
d’autres alliages permettrait de voir dans quelle 
‘* mesure ce modèle est valable. 

Nous supposerons donc que l’impureté ne modifie 
pas les moments magnétiques de la matrice, mais 
introduit un couplage indirect entre deux atomes 
magnétiques diamètralement opposés par rapport à 
elle. Ce couplage de superéchange peut être, sui- 
vant le cas, ferromagnétique ou antiferromagné- 
tique. Nous allons étudier les modifications locales 
de structure magnétique entrainées par cette inter- 
action : 1) le point de transition se déplacera en 
fonction de la concentration d’impuretés ; 2) au- 
dessus du point de Curie d’un ferromagnétique ou 
du point de Néel d’un antiferromagnétique, l’allure 
des corrélations de spin sera modifiée ; 3) à basse 
température, des distorsions se produiront dans le 
réseau des spins ordonnés, autour d’une impureté 
non magnétique. 

Nous étudierons ces effets en fonction du signe 
et de la force du couplage indirect introduit par 
l’impureté, dans certains cas précis d’ordre magné- 
tique : cas ferromagnétique dans les structures c. c. 
et c. Î. c. ; cas antiferromagnétique dans les struc- 
tures c. c. (tous les moments portés par les atomes 
premiers voisins d’un atome donné sont antiparal- 
lèles au moment porté par cet atome) et c. f. c. 
(les moments sont parallèles dans un plan (001) du 
cube. et tous antiparallèles aux moments situés 
dans les plans voisins). Ces structures sont en effet 
assez courantes dans les métaux et alliages anti- 
ferromagnétiques et semblent en particulier s’ap- 
pliquer à Mn Cu [6] (c. f. c.) et peut-être Pu AI [7] 
LENTC:): k 

Dans chaque domaine de température, nous relie- 
rons l’effet du couplage indirect à l’évolution de la 
résistivité électrique. 


II. DÉPLACEMENT DU POINT DE CURIE 7, 
ET DU POINT DE NÉEL 1%. 


1. Cas ferromagnétique. — Soient J l'intégrale 
d'échange direct entre atomes du métal et J’ 
l’intégrale d'échange indirect par les atomes d’im- 
pureté. Dans l’approximation du champ molécu- 
laire, le champ effectif agissant sur un spin donné 
sera : 

JZ : & 


2 
Hnett = Ho LE (1 — 


+ — C{1 = C) + Sy = (1) 
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avec 
£ — facteur gyromagnétique, 
B — magnéton de Bohr, 
Z — coordinance, 
He — champ extérieur, 
< S, > — valeur moyenne des spins. . 


J'est positif ou négatif suivant que le couplage 
indirect introduit par l’impureté est positif ou 
négatif. 

On a : 


1: 
LS ep (OS Het). (2) 


Près du point de Curie, la fonction de Brillouin 
se développe : 


SBsly) 3 (S + 1) y + Only) (3 
Il vient donc 


a 86 
LS, > =38(S +1)É 


(4) 


Pour Æ, — 0, nous avons une solution non tri- 
viale de l’équation (4) pour : 


TT = S(S + 1)Z{t— C)(J + JC). (5) 


Pour 


Il vient donc 


Mad Tee) 


. = cha]. (6) 


J 

Aiïnsi, suivant la nature de l’interaction, la varia- 
tion relative du point de Curie sera-t-elle plus ou 
moins rapide que C. Le point de Curie peut 
même croître avec la concentration si l’intégrale de 


superéchange J’ale même signe que Jet assez forte. 


2. Cas antiferromagnétique. — Nous décom- 
posons le réseau initial en p sous-réseaux, tels 
que les premiers voisins d’un atome sur l’un des 
sous-réseaux soient tous situés sur des sous-réseaux 
différents de ce sous-réseau. Pour une structure 
cubique centrée, p = 2, et pour une structure 


NE =. 
PRES 
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cubique à faces centrées, p — 4. Notons que dans … 


une telle décomposition, l’atome diamétralement 
opposé à un atome d’un des sous réseaux par rap- 
port à un premier voisin de cet atome appartient au 
même sous-réseau. 

On obtient alors pour le champ effectif Hyer, 
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agissant sur un spin S$, (pris sur le sous-réseau n) : 


Haut = Ho + 2 | H— 0 2, < Sn > 
LAUEERS (7) 
et ; 
<s,> =5s( +0 [m4 ao 
È,<Sn> + cu — ci <s, >]. (8) 


Nous obtenons en tout p équations analogues 
pour l’ensemble des p sous-réseaux. Pour A, nul, 
le système des p équations à une solution non tri- 
viale pour 


2 J ’ 
T=Tx = S(S +1) (1 —C) Aer 7 cr]. 
(9) 
Il vient 
TR(0) — Tx(C) , à J’ 
en CU cs Ai]. uo 


Nous voyons qu’une interaction de superéchange 
antiferromagnétique accélérera la diminution de Tx 
avec la concentration d’impuretés. Une interaction 
ferromagnétique produira l’effet opposé et même, 
si elle est assez forte, une augmentation de TX avec 
la concentration d’impuretés. Notons enfin que 
les effets d’une impureté seront trois fois plus forts 
dans la structure c. f. c. que dans la structure €. c. 


_IIT. CORRÉLATIONS DE SPIN 
AU-DESSUS DU POINT DE TRANSITION. 


1. Corrélation de spin. — Soit un site 0 magné- 
tique. Si l’on y fixe un spin S, suivant une cer- 
taine direction, les autres spins prennent, aux 
sites À, une orientation moyenne Sx, liés à S, par 
la fonction de corrélation : 


1 LS; Sr > 
PÉTER à 


Au voisinage de la température de transition, on 
relie facilement Yr à Yr, Où Yr est la fonction de 
corrélation pour les sites À”, premiers voisins de /?. 
Aiïnsi, pour un ferromagnétique, on a, dans l’appro- 
ximation du champ moléculaire et en se limitant à 
un développement au deuxième ordre de la fonc- 
tion de Brillouin (éq. 3) : 


2 ve + OT). (11) 


nm = 72 
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On montre facilement que pour un métal antifer- 
romagnétique 


Y en YE". (12) 
E 

Un calcul de proche en proche montre que les 
corrélations y, intéressant les premiers voisins de 0 
sont en 1/7, que les corrélations y, intéressant les 
deuxièmes voisins sont en (1/7)?, que les corré- 
lations y, intéressant les voisins d’ordre nr sont 
en (1/7)", si l’on ne tient compte que des inter- 
actions avec des voisins d'ordre (n — 1) au plus. 

La substitution à un atome magnétique d’un 
atome non magnétique couplant entre eux les 
moments de deux atomes non magnétiques situés 
de part et d’autre modifie l’aspect des corrélations 
de spin. On voit immédiatement que si l’un des 
sites premiers voisins de 0 est occupé par une impu- 
reté, la fonction de corrélation y, liant le spin de 0 
au spin de l’atome diamétralement opposé sera du 
premier ordre en 1/7. 

Étant données les approximations faites pour 
obtenir les relations (11) et (12), seuls les termes de 
corrélation en 1/T et (1/T)? sont significatifs : 

Les termes en 1/T proviennent des corrélations 
Y1 et Yz:; les termes en (1/T)? proviennent des 
corrélations Yo, Yæ1, Ye. Ces deux termes sont liés à 
l’environnement des sites x : si les atomes premiers 
voisins de + sont magnétiques, les corrélations Y:1 
liant les spins de ces atomes au site central sont 
en (1/7)?, et si certains sites premiers voisins sont 
occupés par des atomes non magnétiques, l’inter- 
action indirecte introduira des corrélations y» en 
(LT) 

Si l’on suppose la tendance à l’ordre loca 
isotrope, ce qui n’est strictement vrai qu'à haute 
température, on détermine aisément à l’aide des 
relations (11) ou (12) les corrélations Y:, Yz, Ya Yen 
Yæ, Correspondant à une concentration C d’impu- 
retés : (voir tableau p. 370) 


2. Application à la résistivité électrique. — La 
présence d’ordre local dans un métal magnétique 
crée un terme de résistivité additionnel qui a été 
analysé par de Gennes et Friedel [1]. Si pom est 
la résistivité due au désordre complet des moments 
magnétiques, la résistivité pn(7) s’exprime de la 
manière suivante : 


AIT, : 
En lT) = pom f. sin 0(1—cos 0) d0 D yr exp (i.Æ.R) 
À R 

(13) 

où > indique une sommation sur tous les sites 
R 

magnétiques. En prenant une moyenne sur toutes 
les orientations possibles du vecteur k, et en nous 


limitant aux corrélations Y3, Ye, Vos Yen Yew 1l 
vient : 


Em(T, C) = Pom(C) [1 + Z1 Y1 gr R1) + Za Ye g(kr Ro) 
+ Ze Ye 8(kr Re) + EYa (KE Ro) + EYar g(Kp Rox)] (14) 
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où kr est le vecteur d'onde au niveau de Fermi et 


ie ; 0 
1 sin (2v sin J 
g(y) = 3 sin O{1 — cos Br een d0 
ë 2y LUS 
un 1 sin 2y 1 — cos 2y 
— de (cos 2y — un PE 
STRUCTURE C. €. 
FERROMAGNÉTIQUE | ANTIFERROMAGNÉTIQUE 
1 Te(0) 1 TN) 
16 SRE SET 
TES 1 J'Tx'0) 
es SÉTEAT NIET 
S AC LU RE FT 
Y2 16 T 16 T 
PNR ÉeeUl UT 
ANT FPE 64 J | ji 
a 1J on 
Ci DOTE QU 64 J | Ti 


o(C) au premier ordre en C, qui se met sous la 


forme 
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e(C) = u(kr D) + C(J'IJ) v(kr D). 


Le calcul des fonctions w(kr D) et v(kr D) est 


]. 15) 


par 


Pm[TelC)1— pom 

Fom 
où T}(C) est la température de Curie ou de Néel de 
l’alliage. On obtient. ainsi un développement de 
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Il est intéressant de connaître l'amplitude maxi- 0 
| male de l’anomalie de résistivité. Elle est donnée r 


= (C) 


STRUCTURE C, F. C. 


ANTIFERROMAGNÉTIQUE 


FERROMAGNÉTIQUE 
Tel0) 4 Tx(0) 
12724 ICT (0) ET EC) T0) 
1 J'Tel0) 3 J'Tx!0) 
LT CAD TT 
10 T&!0) FRAC Tx°(0) 
37- THAT EUES IC) TON 7 TT EUE = CET 
ce Ë | Li È ur 
188 LJ OT AGE JT 
7 [Te(0)]? 3 [Tx(0)]? 
LAN TAIT SA CNT AE 3 


pour 4C (1 — C)?sites 


DE JE [Te(0)]? 


TÎT +54 — 0 Txto)] 
pour 4C (1 — C)? sites 
3 J' [Tx'0)]? 


12 HN T(121 APE C) T0) 
pour 7C (1 — C}? sites 


16 J TIT + (1 —C) Tw0)] 
pour 7C (1 — C}?sites 


développée en annexe [ pour les diverses struc- 


tures envisagées et une représentation graphique 


(16) 


de ces fonctions est donnée (fig. la, 1b, lc 1d). 
Nous avons fait varier kr D entre les valeurs 2 et 5 
ce qui correspond aux cas les plus usuels. On peut 


F1G. 1. — Représentation graphique des fonctions u(kr D) et v(kr D). — (a) Structure c. c. ferromagnétique. — (b) Struc- 
ture c. c. antiferromagnétique. — (c) Structure c. f. c. ferromagnétique. — (d) Structure c. f. c. antiferromagnétique. 


Saut 15 
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noter qu’excepté pour de très petits intervalles de 
variations de *r D les fonctions x et v sont de 
même signe pour les substances antiferromagné- 
tiques à basse température et de signes opposés 
. pour les substances ferromagnétiques. En général 
donc et si J’/J est positif on devra noter une 
augmentation de ®(C) avec la concentration dans 
les substances antiferromagnétiques ; l'inverse sera 
observé dans les substances ferromagnétiques. Ce 
raisonnement s'applique si kr D reste constant 
sinon la situation est plus compliquée. 

D'une manière générale il ressort que les termes 
en (1/7)?et en (1/7) sont comparables au point de 
. transition. On doit donc s’attendre à une décrois- 
sance de l’anomalie de résistivité parfois plus rapide 
que (1/7) au-dessus du point de transition. 


IV. ROTATION DES MOMENTS MAGNÉTIQUES 
VOISINS D’UNE IMPURETÉS 
DANS LA STRUCTURE MAGNÉTIQUE 
ORDONNÉE. 


1. Caleul des rotations — Nous substituons tou- 
jours à un atome magnétique une impureté non 
magnétique qui introduit une interaction indirecte 
entre deux atomes diamétralement opposés par 
rapport à elle. Dans les états magnétiques ordon- 
nés que nous nous proposons d'étudier ces deux 
atomes ont des spins parallèles. Par conséquent si 
Vinteraction indirecte provoque des rotations des 
spins voisins de limpureté on voit immédiatement 
que ces rotations n'auront lieu que si l’impureté 
introduit une interaction d'échange antiferroma- 
gnétique (J" << 0). Nous nous bornerons donc à 
l’étude de ce cas. 

Nous prendrons deux approximations de base : 
champ moléculaire et spins classiques. L’opérateur 
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spin quantique sera remplacé par un vecteur clas- 
sique de longueur donnée #$ mais d’orientation 
libre. 

Par raison de symétrie les rotations des spins 
ont lieu dans le plan défini par les deux atomes et la 
direction non perturbée de leurs moments: dans 
des sens opposés pour chacun des deux atomes. On 
montre aisément que dans tous les cas que nous 
considérons l’axe passant par l’atome non magné- 
tique et portant la direction des moments non per- 
turbés est axe de symétrie de l’ensemble pour la 
configuration résultante stable. La distorsion du 
réseau de spins s’étend en principe assez loin mais 
dans un cas un peu différent de Gennes [8] a 
montré que l'effet d’une impureté diminue comme 
l'inverse du carré de la distance pour deux sites 
premiers voisins initialement perturbés. Comme 
nous avons huit ou douze premiers voisins per- 
turbés nous pouvons penser que l'effet total di- 
minue très rapidement avec la distance. Aiïnsi 
pouvons-nous considérer comme nulles les rota- 
tions de tous les moments portés par des atomes 
plus éloignés que les premiers et deuxièmes voisins 
de l’impureté centrale. 

Nous allons indiquer maintenant les configu- 
rations pour les divers cas proposés. Le calcul est 
facile mais long ; aussi pour simplifier la lecture 
donnerons-nous en annexe II le développement 
complet de ce calcul pour une structure c. f. c. 
antiferromagnétique par exemple. Les distorsions 
sont caractérisées par les rotations €, et e, des 
moments premiers et deuxièmes voisins respec- 
tivement. 


1. CAS ANTIFERROMAGNÉTIQUE (J ET J' NÉGA- 
TIFS). — a) Structure c. c. — Autour d’une impu- 


reté À les moments initialement parallèles à l’axe 
Az prennent l’une des configurations représentées 


(fig. 


2). Elles ont toutes la même énergie. 


Fi. 2. — Structure c. c. antiferromagnétique. Rotation des moments autour de l’impureté À, aux sites premiers et 
deuxièmes voisins, dans les trois configurations possibles de même énergie. 


Les rotations €, sont liées aux rotations €, par : 


foie — (V2/2) sin e,f(1 + cos €) (17) 


ou 


€ (Value, si e, est petit. 
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Outre les solutions &, — kx €, est solution d’une 
équation du quatrième degré en cos e, dont les 
racines ne sont pas immédiates. Si e, est petit, 
alors 

103 
32 [J’[J] —1 


il y a une solution si |[J'[J| > 3,25. 

On vérifie que la configuration stable correspond 
s é ue 103 
de Osi|J|1J7| < 3,25 et à cos Ci 32 J'1J] —1 
si [J'{[J] > 3,25. La solution €, — x n’est jamais 
stable. 

b) Structure c. f. ce. — Il y a deux configurations 
possibles, d’égales énergies, représentées figure 3. 


COS € — (18) 


Fic. 3. — Structure c. f. c. antiferromagnétique. Rotation 
des moments autour de l’impureté À, aux sites premiers 
et deuxièmes voisins, dans les deux configurations pos- 
sibles de même énergie. 
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Nous obtenons 
eee ET sin 6e (19) 
On voit immédiatement que e, est faible, puisque . 
tg es < 2 ss — 0.147. 


Outre les solutions €, — Xx, nous obtenons une 
solution approchée pour cos e; 


1 


———— (20) 
à IJ'1JL+ 0,255 


cos <= 


Cette solution existe si [J’/J| > 0,87. 

On vérifie encore que la solution e = 7, n’est 
jamais stable, que la solution £ — 0 est stable si 
HJ'IJ] < 0,87, et que la solution 

4 
(677) [II] + 0,255 
est stable si [J'[J] > 0,87. 


Ey =" AFC'COS 


2. CAS FERROMAGNÉTIQUES (J POSITIF, J' NÉGA- 
TIF). — à) Structure c. c. — Trois configurations de 
même énergie sont possibles (fig. 4), e, et e, pren- ” 
nent les valeurs calculées en fonction de |J'/J] 
dans le cas antiferromagnétique. La condition de 
rotation des spins reste la même. 


J'1J1 >.3,95. 


F1G. 4. — Structure c. c. ferromagnétique. Rotation des moments autour de l’impureté À, aux sites premiers et deu- 
xièmes voisins, dans les trois configurations possibles de même énergie. 


b) Structure c. f. c. — Deux solutions d’égales 
énergies sont possibles (fig. 5). On obtient : 


DEV EE ( sin €, 
4 2 + cos «, 
et en outre les solutions €, — Xx, la racine d’une 


équation compliquée en cos €, dont la solution 
pour €, petit est : 


Îg Es — (21) 


20,4 
61/17] — 5,9 
La solution existe si [J'[J| > 4,4. 
On vérifie que la solution stable est e, — 0 pour 


20,4 
<< TE — 0 RS ET E Ses 
[JT] < 4,4 et e, — Arc cos GI [IT — 5,9 
TNAL. 


COS e, (22) 


pour 


La condition pour qu'il y ait rotation des 
moments autour d’une impureté semble difficile à 


F1G. 5. — Structure c. f. c. ferromagnétique. Rotation des 
moments autour de l’impureté À, aux sites premiers et 
deuxièmes voisins, dans les deux configurations possibles 
de même énergie. 


N°6. 


N°6 


réaliser dans les structures €. ©. ([J'/J1 > 3,25 pour 
les ferromagnétiques et les antiferromagnétiques) 
et dans les structures ©. f. ©. ferromagnétiques 
([J'1J[ > 4,4). Par contre, dans une structure c. f. e. 
antiferromagnétique type Mn Cu, des rotations de 
spins devraient avoir lieu assez fréquemment 
([J'1J1 > 0,87). 


2. Résistivité électrique aux basses tempéra- 
tures. — Aux basses températures et en particulier 
au 0 °K, si les conditions de rotation des moments 
autour d’une impureté sont réalisées, nous obser- 
verons de fortes résistivités résiduelles, d’autant 
plus fortes que la perturbation s’étendra plus loin. 
Si la perturbation ne s’étend de façon notable 
qu'aux premiers voisins d’une impureté, nous 
aurons C(Z + 1) centres diffusants pour un élec- 
tron de conductibilhté. Le potentiel diffusant est 
proportionnel aux cosinus des angles € de rotation 
des moments, et si l’on suppose que lapproxi- 
mation de Born est valable, c’est-à-dire que ce 
potentiel est faible vis-à-vis de l’énergie de Fermi 
des électrons, la section efficace de diffusion sera 
donc porportionnelle à (1 — cos £)2. Nous pouvons 
relier la résistivité 9, au 00 absolu à la résistivité pom 
de haute température où le désordre des spins est 
total et nous obtenons, en supposant que seuls les 
premiers voisins sont perturbés : 


Pa __ C(z + 1)%5 (1 — cos €)? 
Let NOT 
À cos? 0 d6 

7% Jo 


Ceci n’est valable qu’à très faible dilution, 
lorsque les centres perturbateurs n’interagissent pas 
entre eux. On conçoit que si CZ est voisin de 
l'unité, la résistivité au 0 ©K soit égale à celle de 
haute température. Dans un tel modèle, la résis- 
tivité résiduelle ne sera donc proportionnelle à la 
concentration d’impuretés que pour une très grande 
dilution. 

Nous pouvons prévoir un intéressant effet dans 
les substances métamagnétiques c. f. ©. Si 

0,87 < [J'[JI < &,4 

il y a rotation des spins dans l’état antiferroma- 
gnétique et cette rotation n’a pas lieu dans l’état 
ferromagnétique. On devrait observer alors une 
sensible anomalie de résistivité à la transition. Ce 
point pourrait être facilement vérifié dans les solu- 
tions solides des terres rares métamagnétiques de 
structure c. f. c. La condition de rotation ne doit 
pas être très différente pour une structure antiferro- 
magnétique d’un type autre que celui que nous 
avons étudié. 


= 2C{Z + 1)2/3 (1 — cose)? 


Pom 


V. APPLICATION } 
A QUELQUES STRUCTURES MAGNÉTIQUES. 


1. Les alliages Mn Cu quadratiques à faces cen- 
trées. — La phase y du manganèse, cubique à faces 
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centrées, peut être stabilisée à basse température 
par des additions de cuivre. Elle prend une défor- 
mation quadratique à faces centrées que l’on a pu 
associer par des mesures de diffraction de neutrons 
à une structure antiferromagnétique, où les spins 
sont portés par Oz et ont des directions opposées 
d’un’plan (001) à l’autre. La température de Néel de 
ces alliages est, pour C faible, une fonction linéaire 
de la concentration ; si T'x(0) est l’extrapolation du 
manganèse pur : 


TR(0) — Tac) _ 
TT 


En appliquant notre analyse, on en déduit :, 


4 — 3J'IN] = 3 et donc J° = — (2/3) [J|. 
Un ion Cu introduirait donc un couplage antifer- 
romagnétique entre deux ions Mn. La condition 
pour qu’il y ait rotation des spins à basse tempé- 
rature n’est pas réalisée (|J’/J] > 0,87). Effecti- 
vement, la résistivité résiduelle du cuivre dans le 
manganèse reste proportionnelle à la concentration 
jusqu’à des teneurs assez élevées en cuivre (20 %) 
et la section efficace de diffusion d’un atome de 
cuivre est de l’ordre de l’aire atomique (x (1,9 À)?). 


2. Les alliages plutonium-aluminium eubiques à 
faces centrées. — La structure magnétique de ces 
alliages n’a pas été obtenue avec certitude. Mais les 
mesures de résistivité [9] et de chaleur spéci- 
fique [10] semblent indiquer une transformation 
ordre-désordre en dessous de 200 °K. Les mesures 
de susceptibilité magnétique n’excluent pas un 
antiferromagnétisme à basse température. Nous 
désignerons par 7% la température de transfor- 
mation (extrapolée au plutonium pur : Tx(0)) 


(TN(0) — Tx(C))ITx(0) = CH + 8(J'[J) (1 —C)]. 
On en déduit, si notre analyse s’applique : 
TERRES 


Les courbes de résistivité présentent une ano- 
malle qui a été associée à un antiferromagnétisme 
possible à basse température [7], [10]. Un 
Em(Tn)/6om à été calculée [7], et il varie de la 
manière suivante avec la concentration d’alumi- 
nium, 


04 A1 2,25 3 5 [10 
T 

EmÜTN) 408 41107 113 414 
Pom 


Si l’on ne tient pas compte de l’interaction indi- 
recte, augmentation du rapport pn(7x)/6om doit 
être liée à une diminution de (kr D) de 3,5 à 3,3 
lorsque la concentration d’aluminium croît, c’est-à- 
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dire à une diminution du nombre d’électrons de 
conductibilité, voisin de l’unité pour le plutonium 
pur [7]. C’est cette interprétation que nous avons 
donnée précédemment. Si par contre on tient 
compte d’une interaction indirecte par l’intermé- 
diaire des atomes d'aluminium, il apparaît, d’après 
les résultats de la section IIT, et en utilisant la 
valeur de J'/J — 1,8, que kr D reste sensiblement 
constant lorsque la concentration d’aluminium 
varie ; on trouve alors : 


9 AI 2,25 3 5 10 


kr D 95 3475 3425, 3478 

La valeur moyenne de kr D — 3,475 indique la 
présence de 1 électron de conductibilité par atome 
dans le plutonium. Ce nombre dépendrait peu de 
la teneur en aluminium. 

La décroissance de l’anomalie au-dessus de la 
température de transition est sensiblement plus 
rapide que (1/7). Ce phénomène est bien en accord 
avec le fait que les coefficients des termes en (1/7) 
et (1/7)?dans le développement de px(7%) /pom Sont 
de même signe et du même ordre de grandeur pour 
la valeur de kr D — 3,475. 

Un ion aluminium introduirait donc un couplage 
antiferromagnétique entre deux ions plutonium et 
la valeur J'/J — 1,8 trouvée satisfait à la condition 
pour qu'il y ait rotation des spins à l’état ordonné. 
Effectivement, on observe expérimentalement de 
très grandes résistivités résiduelles au 0 °K : à 100% 
d'aluminium (et sans doute dès 80/) la résistivité 
au 0 °K est égale à la résistivité de haute tempé- 
rature. Pour l’alliage à 5 °/, d'aluminium, on trouve 
une section efficace de diffusion par centre pertur- 
bateur égale à x (3,6 À )2. Comme rs — 1,8 À, ceci 
montre bien que la perturbation intéresse les pre- 
miers voisins de chaque atome d'aluminium. De la 
valeur de |J’/J| calculée à partir du déplacement 
du point de transition, on déduit effectivement une 
déviation forte des spins sur les premiers voisins 
seulement : 

E, © 9590 


En = 200) 


3. Alliages de terres rares. — Il n’y a pas eu de 
recherches systématiques des déplacements des 
points de Curie ou de Néel par des impuretés dans 
les terres rares. Cependant les résistivités rési- 
duelles au 0 0K, même pour les métaux purifiés, 
sont élevées, ce qui suggère que notre modèle 
pourrait s’appliquer. Il serait intéressant de vérifier 
ce point. 


VI. Conclusion. — Une interaction indirecte 
ferromagnétique ou antiferromagnétique par l’in- 
termédiaire d’un ion non magnétique couplant 
entre eux deux ions magnétiques diamétralement 
opposés par rapport à lui, semble pouvoir expli- 
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quer les propriétés de certains alliages riches en 
ions magnétiques : abaissement des points de Curie 
ou de Néel, variation de la résistivité au-dessus du 
point critique, fortes résistivités résiduelles au 0 °K. 

Une étude systématique des résistivités rési- 
duelles au 0 °K de métaux magnétiques dilués en 
ions non magnétiques permettrait de préciser dans 
quelles conditions notre modèle s’applique. Des 
métaux à magnétisme de couche, ou voisins de ce 
type, comme les terres rares ou le manganèse, 
semblent particulièrement indiqués pour ces vérifi- 
cations, qui, si elles étaient en accord avec notre 
modèle, permettraient de préciser la nature de cet 
échange indirect. Uné interaction de ce type a été 
envisagée dans les alliages Fe AI à forte teneur en 
aluminium, compris entre Fe; Al et Fe AI [11]. 

Ajoutons que l’analyse des fonds continus des 
spectres de diffraction de neutrons à basse tempé- 
rature permettrait éventuellement de mettre en 
évidence les rotations de spin dont nous avons 
calculé l'amplitude à 00K. 

L'auteur tient à exprimer sa gratitude à M. le 
P* J. Friedel qui l’a constamment guidé pendant la 
préparation de ce travail. 


Annexe I : Calcul du terme de résistivité due à 
l’ordre local. — Nous devons expliciter les coeffi- 
cients de l’équation (14). Les valeurs calculées 
pour Z et À sont portées sur le tableau suivant, 
où D est la distance interatomique. 

En utilisant les notations 


(1 — C) Te(0)[T = 1/0 

A — CT OT = 18 
g(kr D) = g(1) 

g(2kr D( = gf2), ete. 


Nous obtenons immédiatement les résultats sui- 
vants pour les divers cas étudiés : 


1.-STRUCTURE €. ©. 
— Cas ferromagnétique. 
1 


Em(7', C) li JS 9 
— À Le o(4\ - — p\2 =. 
ds Hg [et + See +2 


IV) + Let) + SVT) + SAV) 
C2 : ee 

+ WIN Ee() + 8e(4lV3) + ag(a V2F3] I] (24) 

glC) = 8h) + (23) + cir in) | 


— ali) —Ÿ el2lV5) + (1/64) 


Cats) 34(VTI7S) + (VTT (5) 
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STRUCTURE C. C. STRUCTURE C, F. C. 
Z R 2 R 
Zn HEC) D 111 0) D 
FN SU 0) 2D 42C{1 0) 2D 
air PO 0) 21V5 D Bic V2 D 
Ze fact CVs DU | EC CVS D 
3C(1 — C}? V19/3D 3CH.—C)? IMVE.D 
LG" 3D SOU CYR MY LD 
20(1 — C)? | 2V3 D 
| Zen, Re | 3C4— Ch | É]V3D |aci—c| 2D 
3C2{1 — C)? | 4/V21V3 D | 20241 — c)? | 2V2D 
&C?{1 — C)? | 2V3D 
C?{1 — C)? 4D C?{1 — C)? 4D 
— Cas antiferromagnétique. — Cas antiferromagnétique. 
ÉRCENRRER JÉPE TT PaiT CSS PE Pre Peas 1 
Eu | «1 EST Éelss a O ee mnt doit 04) 
3 NC ee - : = 
(Le 2V3) + Les) + 3e(V1178) + 3(V1973)] |Se(v5) + CT De(V5) + 2e(V7) +28(2 V3)] 
C2JJ\2 4 5 3 J’ = 
Re sa 6e) +0 (VA) damien ce 
a ER 
3 f/2 " de lGre 
he el} sl 1\+0c 1 g([— 1607 
(C) g(1) ue T g() + Fa 4 


met + | et8) + 8e (JA) + 8e 


2. STRUCTURÉ C. F. C. 
— Cas ferromagnétique. 


Em(T, C) … 3 Û 
Pom(C) SE 


0e(300 — 1) 


[3e(V5) + 2e(V3) + 2e(2 va)]} 


1 M ré < 
Foto 1 03 V3) 


CAE re ea Arr 1 


1 1 = sl 
+ nn 3 V0) AC 


A 
CET 


“ re CR [ag(2) + 2g(2 V2) 


144 J 6e 


+ 4g(2 V3) + g(4)] 


gl0) = À alt) +7 e(V2) 
" cr | ral) — 2 8(V5) + . g(V3) + (2 
+2 (V5) +2 (eV?) + el Va). 


(28) 


Cag(2) + 2g(2 V2) + 4g(2 V3) + g(4)] (30) 


(VAE Fa o(C) = —© sl) + r g(v3) + FRE 


F-; g(1) + a(V3) + 5 (V3) — 9g(2) 


3 (V5) + (V7) + Eee V5)]. (31) 


Annexe II : Calcul de la rotation des moments 
autour d’une impureté. — Nous calculons les rota- 
tions e, et e, des moments portés par les premiers 
et deuxièmes voisins de l’impureté, en supposant 
que les rotation de tous les autres moments sont 

nulles. Il y a deux étapes dans le calcul : 

1. Les rotations e, des deuxièmes voisins sont 
déterminées en fonction des rotations €, des pre- 
miers voisins, en minimisant l'énergie d'échange de 
chaque deuxième voisin avec tous ses voisins. 

2. On calcule la variation d’énergie d'échange 
pour l’ensemble de tous les moments perturbés ; en 
minimisant celle-ci par rapport à e;, on obtient la 
valeur de &;, puis de e. 

Nous allons effectuer ce calcul pour une sub- 
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stance antiferromagnétique de structure c. f. c., 
type Mn Cu. Dans l’état parfaitement ordonné, les 
moments sont portés par la direction [001], sont 
parallèles dans un plan (001) et tous antiparallèles 
aux moments portés par les atomes situés sur les 
plans voisins. Notons par (+) et (—) les deux types 
de moments. A un site (+), substituons un 
atome À non magnétique. La configuration la plus 
stable correspond à la fermeture des moments (+) 
vers l’axe Az et à l’ouverture des moments (—), 
soit à la solution inverse. 

Calculons les rotations des deuxièmes voisins 
de À. Ce sont des moments (+) et seuls les quatre 
moments situés dans le plan (001) passant par À 
pivotent. L’énergie d'échange d’un deuxième voi- 
sin (+) avec ses premiers voisins (+) et (—) est : 


Eex =2|J|S(S +1) [2 (2 e sin €, sin €2 + COS €, cos «.) 
— 2(F sin €, sin €, + cos €, cos €;) — 4 cos «|. (32) 


La configuration de plus basse énergie est obte- 
nue pour une fermeture des moments (+) des deu- 
xièmes voisins si les moments (+) des premiers 
voisins se sont fermés, 
premiers voisins se sont ouverts. Alors : 


Eex = 2 [J|S(S + 1) [—sin e, sin (2 Va) = rcose|. 
(33) 
On minimise par rapport à & 
dEez 
ae = 2MIS(S + 1) 
[— sin € COS (2 — V2) + 4 sin cl 007) 
d’où 
OR 
L 


tg £a = sin £2. (35) 


L’énergie d’échange pour l’ensemble des mo- 
ments perturbés s'écrit : 
Eez tot = 12 [J”| S(S + 1) cos 2e, — 56 [JIS(S + 1) cos €, 
— 16 [JIS(S + 1) (V2 sin? e, + cos? €) 
— 16 [JFS(S +1) (: sin €, sin € 


— 16 [ZI S(S + 1) cos &. (36) 


ou une ouverture si les 


N° 6 


On remplace €, par sa valeur en fonction de & 
et on minimise par rapport à e.. Il vient 


dEez tot 
ns Er 

de; ] 

S(S + 1) sin €, cos €, + 56 [JIS(S + 1) sin e 
Ê ee =} 
z sin €; COS €1 
9\2 1/2 

[1 + C _ sin? «| 


3\2 
) sin €; COS & 
(37) 


— [48 |J1 + 32 [J1(V2 


8217] SE) 


16 |J|.5(S + 1) É 


k 
SNS 3/2 
E ne sin? «| 
Nous trouvons d’abord les solutions €, = kr. 
L’autre solution est donnée par : 


—= 10. 


cos €; i8| 7] 1 39(V9 1) 


Rd 
| ve 41 fe) sin? | ‘+ 


En négligeant [(2 — 4/2)/4]° sin? e, devant 1, 
on a une solution très bien approchée en cos €; 
1 
(6/7) [J’JJI + 0,255° 


Cette solution existe |[J'/J] > 0,87. 
On vérifie facilement que la solution stable est 


[J'[J| < 0,87 


€, = Arc cos Re. si JT 087 
; IJ'I1 + 0,255 


Le calcul s'effectue de la même manière dans les 
autres Cas, 


COS € (39) 


€ = 0, si 


Manuscrit reçu le 11 janvier 1961. 
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DÉTERMINATION THÉORIQUE ET EXPÉRIMENTALE DU RAYONNEMENT DE FREINAGE 
ÉMIS PAR LES g@ DU PHOSPHORE 32 DANS L'OR 


Par R. QUIVY (*) et J-FRANEAU (*#), 


Résumé. — Les formules théoriques, relativistes et non relativistes, donnant la section efficace 
d'émission de rayonnement de freinage, pour des électrons monoénergétiques, sont appliquées à la 
prédétermination du rayonnement émis lors du freinage des 8 du phosphore 32 dans l'or. 

Pour étudier expérimentalement ce rayonnement de freinage, on établit une méthode d’inter- 
prétation des spectres relevés par un compteur à scintillations. 

L'expérience montre que les prévisions théoriques se vérifient de façon satisfaisante pour des 
radiateurs de faible épaisseur et que, dans le cas de radiateur épais, le rayonnement de faible 


énergie augmente considérablement d’intensité. 


On peut en déduire quelques considérations relatives à l’utilisation pratique du rayonnement 


de freinage. 


Abstract. — Relativistic and non relativistic theoretical formulas, giving the bremsstrahlung 
cross-section for monoenergetic electrons, are used to compute the bremsstrahlung produced by 


P 32 £ rays in gold. 


In order to study experimentally this radiation, a method is proposed to interpret experimental 


spectra obtained with a scintillation counter. 


Experiments agree satisfactorily with theoretical predictions for thin targets. 


With thick 


targets, the intensity of the low-energy radiation increases considerably. 
The practical use of this bremsstrahlung is discussed. 


I. Introduction. — On sait qu’un électron accé- 
léré par le champ d’un noyau dans le voisinage 
duquel il passe, émet de l’énergie sous forme de 
photons au détriment de son énergie cinétique. Ce 
rayonnement de freinage (bremsstrahlung) a été 
l’objet de nombreuses études théoriques. 

Nous nous sommes proposé d'étudier expérimen- 
talement le rayonnement de freinage émis par les 
électrons d’une source £ à spectre continu, lors de 
leur ralentissement dans un radiateur, en vue de 
comparer nos résultats aux prévisions théoriques 
et de voir dans quelle mesure la théorie permet de 
prédéterminer le spectre émis. 


II. Considérations théoriques. — Comme émet- 
teur B, nous avons utilisé le phosphore 32. Le fait 
que l’énergie des f varie de façon continue Jusque 
1,7 MeV nous place dans une situation assez déli- 
cate au point de vue théorique. En effet, la théorie 
ne donne la probabilité d’émission d’un photon Av, 
sous une forme assez simple pour se prêter au 
calcul, que dans le cas d’énergies nettement rela- 
tivistes ou dans le cas d’énergies nettement non 
relativistes. Pour les énergies relativistes, nous 
avons utilisé la relation de Bethe et Heïtler [1], 


(*) Chercheur au Centre de Mons de l’Institut Interuni- 
versitaire des Sciences Nucléaires (Belgique). 

(**) Professeur à la Faculté Polytechnique de Mons 
(Belgique). 


dérivée du travail de Sauter [2], et donnant la 
section efficace différentielle : 


TR NE E \° 
COCA SE 137 Fe V [ (x 


_— 
dE, 
9 
2E, E 1 ü) 


à (18 me2.hv 9 
énergie totale initiale de l’électron inci- 
dent, 
Z : nombre atomique du radiateur, 
Æ  : énergie totale résiduelle de l’électron 
après émission du photon Av: 
E =E,— h. 
Pour les énergies non relativistes, nous avons 
appliqué la formule de Kirkpatrick et Wied- 
mann [3] : 


avec : E;: 


To = Iy sin?0 + Zy(1 + cos°0) (2) 


qui donne l’intensité totale de rayonnement de 
freinage émis dans la direction 0 (fig. 1) en fonction 
des composantes Z;(v, E,) et Z,(v, E:). Ces deux 
composantes sont d’expression très compliquée 
mais Kirpatrick et Wiedmann les ont tabulées en y 
faisant intervenir des variables réduites qui en 
facilitent l’application ; le calcul numérique des 
composantes /, et Z, a conduit les auteurs précités 
à faire des hypothèses qui ne sont correctes que 
pour des énergies £, non relativistes. 
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Electron incident 
A —————————————— 


Hicele 


D’autre part, nous avons comparé, dans le do- 
maine non relativiste, la relation de Kirkpatrick et 
Wiedmann à celle de Schiff [4], dans la forme 
donnée par Kock et Motz [5]. 

Z2yne?2\ 20% 167? E 

Ga SE 


do(E, v, () = 739 Fa Te y? + 1) He 


a 
+08 [RE WE 


x log. My) | y.ay (3) 


avec : Yu Ed 
et 


1 hy \? zaB  \2 
M{y) ee Cri + ï) ‘ 

Cette dernière relation (3) intégrée à tout l’angle 
solide est en bon accord avec la relation (2). Nous 
avons utilisé la relation (2) parce qu’elle est d’une 
application beaucoup plus aisée à un cas concret. 
Il nous reste à déterminer quelle partie du spectre £ 
doit être considérée comme non relativiste. A cette 
fin, posons : 

O(E,, v) — ER 

et traçons les courbes donnant ®(Æ,, v) en fonction 
de v pour des valeurs croissantes de Æ£;, par appli- 
cation des formules (1) et (2). On voit sur la 
figure 2 que pour E, = 1,5 m,c°, les courbes (1) 
et (2) sont assez voisines. On appliquera donc la 
formule (2) pour moe? < E, < 1,5 mc? et la for- 
mule (1) pour Æ, > 1,5 mc? ; la valeur 


Es = 4,5 me 
correspond à une énergie cinétique de m,c?/2 
(c’est-à-dire 256 keV), ce qui est bien la limite géné- 
ralement admise du domaine non relativiste. 


IIT. Mesures expérimentales. — 1. Disposirrr 
DE MESURE. — La figure 3 schématise le dispositif 
expérimental, La source de %?P, sous forme de 
phosphate de sodium, est déposée par évaporation 


dans un petit cylindre creux en aluminium. Ce dis- . 
positif donne un faisceau 8 dont le demi-angle | 
d'ouverture est inférieur à 80. La nature et les di- 
mensions du support de source et du collimateur 
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_— 
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N 
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256 canaux 
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ont été choisies de façon à rendre aussi faible que 
possible le rayonnement de freinage émis par les 8 
absorbés. 

La feuille de métal servant de radiateur est main- 
tenue dans le faisceau B, juste à la sortie du colli- 
mateur, Le dispositif source-collimateur et le détec-. 
teur sont entourés d’un blindage de plomb de 10 em 
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d'épaisseur afin de réduire au minimum la détec- 
tion des photons provenant directement de la 
source et du collimateur. Le système de détection 
est constitué d’un scintillateur en iodure de sodium. 


2. CORRECTIONS ET INTERPRÉTATION DU SPECTRE 
RELEVÉ. — La détection de photons monoéner- 
gétiques au moyen d’un scintillateur ne présente 
guère de difficulté ; en effet, l’examen de la raie 
photoélectrique fournit les renseignements néces- 
saires à des mesures quantitatives, connaissant la 
géométrie du montage, les dimensions et les coef- 
ficients d'absorption du scintillateur. Appelons : 
N, : le nombre de photons d'énergie Av incidents 
au scintillateur et supposés être de même 
direction ; 

Nn : le nombre de coups enregistrés dans la raie 
photoélectrique ; 

uw : le coefficient total d'absorption du matériau 
constituant le scintillateur (effet photo- 
électrique + effet Compton + création de 
paires), pour les photons d’énergie À ; 

Ua : le coefficient d’absorption totale de l'énergie hv 
par le scintillateur ; 

L : la longueur du scintillateur (parallèlement au 
flux de photons incidents). 

On établit facilement : 


ANph = — (Upnfut)-AN = ph: No-e Ut. 
Non = (uphfut).(1 — el). Ni. 
On pose : 


1/Foù = (uphlut).(1 — et) 
d’où : 


No — Dh IVph- 


Il est bien évident que pour appliquer ces rela- 
tions, il faut disposer des valeurs de 4 relatives au 
type de scintillateur utilisé (nature et dimensions 
géométriques). 

A partir des courbes du N. B.S. [6], nous avons 
tracé la courbe de la figure 4 donnant Fa en fonc- 
tion de l’énergie hv des photcns incidents. 

Cependant, lors de la détection par scintillateur 
d’un flux de photons monoénergétiques, 1l y a, en 
decà de la raie photoélectrique, un spectre continu 
dû essentiellement à l’effet Compton : d’une part, 
effet Compton à l’intérieur du scintillateur avec 
absorption par celui-ci de l’énergie de l’électron 
Compton et échappement du photon diffusé; 
d’autre part, effet Compton à l’extérieur du scin- 
tillateur avec détection du photon diffusé mais non 
de l’électron Compton. 

Dans le cas de photons d’énergies différentes, 
aux impulsions données par absorption totale d’une 
même énergie, se superposeront celles provenant de 
l’effet Compton des photons d’énergie plus élevée. 

On pourrait espérer prédéterminer théori- 
quement la contribution, à la formation du spectre 
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Foh, 


Energie (MeV ) 


FIG 


mesuré, des électrons Compton créés à l’intérieur 
du scintillateur. Malheureusement, il est impossible 
de prévoir l'influence des effets Compton extérieurs 
au scintillateur. Ceci nous a conduit à procéder de 
la façon suivante. 

Nous avons déterminé expérimentalement la 
réponse du scintillateur pour des y monoénergé- 
tiques de différentes énergies et nous avons inter- 
polé pour prédéterminer le spectre donné par des 
y d’énergie intermédiaire. 


Cs 


Epaisseur support de source 


0 680 
(mg /Æmè) 


RG: 


Pour faire ces déterminations, nous devions dis- 
poser de sources y aussi monoénergétiques que 
possible. Il fallait donc que les influences de l’épais- 
seur de la source et du support de source soient 
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négligeables. Les hautes activités spécifiques que 
l’on peut facilement obtenir permettent de résoudre 
le problème en ce qui concerne l’épaisseur de la 
source, D’autre part, étant donné la distance qui 
sépare la source du scintillateur, on devait 
s'attendre à ce que seule, la partie du spectre relevé 
correspondant à la diffusion y à 1809 soit affectée 
par l’épaisseur du support de source, ce qui a été 
confirmé par l’expérience. Nous avons donc relevé 
des spectres, avec une même source y, pour des 
épaisseurs de support de source croissantes. De 
façon à pouvoir comparer nos résultats, nous avons 
calculé le rapport À de l’aire de la raie de diffusion 
vers l’arrière à l’aire de la raie photoélectrique, nor- 
malisé à l’unité pour une épaisseur de support de 
source de 1 mg/em?. La figure 5 donne les courbes 
de À en fonction de l’épaisseur du support de 
source pour le cobalt 60 et le césium 137. On voit, 
à partir de ces résultats, que l'emploi de support de 
source dont l’épaisseur est de l’ordre de 5 mgr/cm? 
évitera toute déformation de spectre due à la rétro- 
diffusion, le seintillateur étant placé à une distance 
de la source grande par rapport au diamètre du 
scintillateur. 

Pour déterminer le spectre Compton de y mono- 
énergétiques, nous avons utilisé le mercure 203, le 
chrome 51, l’or 198, l’iode 128 et le césium 137, 
sous forme de source d'épaisseur négligeable, dé- 
posée par évaporation sur film de formvar de 
5 mg/cm?. 


0,66 MeV 


9% de l’énergie de la raie photo-électrique 
F1. 6. 


Les différents spectres relevés ne sont pas sem- 
blables, ce qui nous a amenés à les normaliser en 
fonction d’une énergie réduite obtenue en divisant 
les abscisses d’un même spectre par l’énergie de la 
raie photoélectrique. 

La figure 6 donne ces spectres ; en abscisse, on a 
porté l’énergie réduite, la position de la raie photo- 
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électrique correspondant donc à 100 % ; en ordon- 
née, le facteur de réponse Fr de notre dispositif, 
c’est-à-dire, les ordonnées du spectre relevé divi- 
sées par l’aire de la raie photoélectrique. Il faut 
signaler que, une fois les spectres normalisés, les 
raies photoélectriques coïncident pratiquement. 

À partir de la figure 6, on peut tracer un réseau 
de courbes donnant le facteur de réponse F; en 
fonction de l’énergie du photon incident. 

L’interpolation graphique permet de compléter 
le réseau de la figure 6 et finalement de tracer 
l’abaque de la figure 7 où nous portons le facteur 


F. 10% 


0 20 #40 60 80 100 120 KO 160 180 200 


de réponse F; en fonction de la hauteur d’impul- 
sion (exprimée en nombre de canaux de notre ana- 
lyseur R. C. L. à 256 canaux) pour une position 
déterminée de la raie photoélectrique notée en para- 
mètre sur les courbes. 

En outre, pour étudier un spectre continu, nous 
devrons le découper en tranches d’égale largeur. 
Or, la largeur de raie photoélectrique variant avec 
l'énergie, il est bien évident qu’une tranche de 
largeur constante axée sur le sommet d’une raie 


Nraie 


N0 canaux réels centrés sur sommet raie 


hY (MeV) 
de la raie 


dis: 
: 
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photoélectrique n’embrassera pas la même portion 
relative de celle-ci selon qu’elle est à basse ou à 
haute énergie, ce qui nécessite l’apparition d’un 
nouveau facteur de correction r (fig. 8) : 


ne Nraie 
Naans la tranche de largeur AZ centrée sur la raie . 


D'où : 
NAE = Nraie [r. (4) 


Un exemple illustrera la façon de procéder pour 
corriger un spectre relevé. 
Considérons un spectre relevé quelconque (fig. 9). 


N 


! Energie h? (canaux) 


FiG. 9, 


On divise le spectre en tranches de largeur cons- 
tante AE (dans notre cas, AE — 4 %, de l’échelle 
totale d’énergie de l’analyseur). Nous pouvons 
admettre que la dernière tranche AE, centrée sur Æ; 
est due uniquement à l’effet photoélectrique ; soit 
N, le nombre de coups enregistrés dans cette 
bande. 

La figure 4 nous donne Fn = F;, pour kv — E:. 
Il suffit donc de multiplier le nombre de coups NW; 
de la tranche Æ, par F, pour obtenir le nombre de 
photons incidents d'énergie comprise entre 

E; — (AE 2) et E, + (AE /2). 

A ces W, coups de la bande AE centrée sur E;, 
correspond un certain nombre de coups dans la 
bande AE centrée sur £, (E, << E;). A la figure 7, 
le point d’abscisse Æ, situé sur la courbe portant 
en paramètre Æ,, a comme ordonnée #,,, ; d'autre 
part, la figure 8 donne r — r, pour £ — F;. La 
contribution de la bande AE centrée sur Æ; est 
donc : 

N,.Friofr, en vertu de (4). 


Il faudra donc soustraire cette contribution du 
nombre de coups W,, enregistrés dans la bande AE 
centrée sur Æ,. Il suffira d'appliquer le même calcul 
aux tranches successives d’énergie décroissante à 
partir de l’énergie maximum Æ£;, chaque tranche 
ayant été au préalable corrigée de la contribution 
des photons d’énergie supérieure à son énergie 
moyenne. 
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A la figure 10, nous avons représenté un spectre 
relevé et le même spectre corrigé de la réponse du 
scntillateur et des phénomènes parasites. 
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200 
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La méthode d'interprétation que nous avons 
employée tient compte de toutes les interactions de 
photons avec le dispositif expérimental (support de 
source, collimateur, air, blindage, coiffe de scintil- 
lateur, etc). 


3. MESURE DU RAYONNEMENT DE FREINAGE. — 
Grâce au dispositif de la figure 3, nous avons étudié 
le rayonnement de freinage émis par les 6 du *?P 
dans l’or en fonction de 0. 

2s mesures ont été faite pour 0 — 390, 450, 600, 
759, notre dispositif expérimental ne nous per- 
mettant pas d'explorer le rayonnement sous des 
angles supérieurs à 750, 

Nous avons utilisé comme radiateur, des feuilles 
d’or de 19, 315%39,%44, 63, 88/1165, 455,074h%°et 
918 mg/cm? d'épaisseur. 

Tous les spectres relevés ont été corrigés de la 
réponse du dispositif expérimental, de l'absorption 
dans le radiateur et de la décroissance de la source 6 


IV. Comparaison et discussion des résultats. — 
1. PRÉDÉTERMINATION DU SPECTRE THÉORIQUE DE 
RAYONNEMENT DE FREINAGE. — L’épaisseur de la 
source de *?P utilisée provoque une importante 
déformation du spectre B. Aussi, nous avons dû 
relever, au spectrographe 8, le spectre 8 réellement 
émis par notre source. Il est représenté à la 
figure 11. | 

Nous avons appliqué à ce spectre 8 les équations 


- de Sauter et de Kirkpatrick et Wiedmann selon le 


principe que nous avons exposé plus haut, ce qui 
nous donne le spectre théorique du rayonnement de 
freinage relatif à une géométrie de 4x (fig. 12). 


2. COMPARAISON AVEC LES COURBES EXPÉRI- 
MENTALES, — Nous avons intégré graphiquement 
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Fc. 12. — Spectre théorique du rayonnement de freinage 
émis par les 8 du *P. 


nos résultats expérimentaux ; aux figures 13a, 13b 
et 13c, nous portons 


6= 82030" 
1 hv.N.dQ, 
6 = 22030" 


en fonction de l’énergie. 
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Sur chacun de ces trois graphiques, nous avons 


tracé la courbe théorique en la normalisant avec la 
courbe expérimentale pour une énergie corres- 
pondant au 140° canal de l’analyseur. 


82030" 
,hNdû 
at fa a) 72P +39 mg Jr d'or 
ET 
100.1 


82930" 


l2=22°30" bd by 2P+ 88 mg /em? d'or 
I 


200. 
I — : Sauter et Kirkpatrick 
I Due I : Expérimental 
EI I 
100.10? 
0 Energie h3 (canaux ) 
50 100 150 200 250 
tn 
da 
6= 22°30° I 
400 
c) 32p +455 mg /em? d'or 
h> (keV) 
300 I 600 


Hr6212: 


Dans le cas des faibles épaisseurs de radiateur, on 
voit que les courbes expérimentales et théoriques 
ont bien la même forme. En outre, l’examen des 
courbes montre que l’augmentation de l’épaisseur 
du radiateur ne favorise que très peu l’émission de 
rayonnement de freinage d’énergie assez élevée, 
celle-ci ayant surtout lieu dans les toutes premières 
tranches du matériau (fig. 14, 15 et 16). Le décolle- 
ment, aux basses énergies, de la courbe expérimen- 
tale par rapport à la courbe théorique est d’autant 
plus important que le radiateur est épais, ce qui 
était prévisible puisque nos équations théoriques de 
départ sont relatives à une seule interaction par 
électron, donc rigoureusement applicables dans le 
cas d’un radiateur d'épaisseur monoatomique seu- 
lement. En effet, au fur et à mesure que l’électron 
ralentit, il émet des photons d’énergie de plus en 
plus faible, ce qui enrichit fortement le rayonne- 
ment mou du spectre expérimental d’un radiateur 


az minister 
NL 


ble 
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épais par rapport à la prédétermination relative à 
une seule interaction par électron. 

Pour vérifier la théorie, il était donc nécessaire 
d'utiliser des radiateurs aussi minces que possible 
tout en ayant un nombre de photons suffisant, 

Wyard [7]et plus récemment Goodrich et coll. [8] 
ont appliqué la théorie de Sauter en l’étendant à 

toute la trajectoire de la particule 8 et ont établi : 
RENE mc) [: (: > ] il 


2 2 
0 — Mo C Es — moc 


x Log re Le a £ ] 
où Z(hv) représente l’intensité du rayonnement de 
freinage de fréquence v émis tout le long de la 
trajectoire d’un électron d’énergie totale initiale £,. 

Nous avons intégré cette formule de Wyard à 
tout le spectre B de notre source de #?P (fig. 13c). 

On voit que le spectre prédéterminé par la for- 
mule de Wyard présente une grande divergence 
par rapport au spectre expérimental. Une augmen- 
tation de l’épaisseur de radiateur jusqu’à une valeur 
proche du parcours maximum des B entraînera évi- 
demment une augmentation du rayonnement mou 
mais par contre le rayonnement dur ne sera que 
très peu affecté ; en effet, il faut remarquer que 
l'épaisseur de radiateur de 455 mg/cm? (fig. 130) 
correspond à 60 % environ du parcours maximum 
des 8 du phosphore 32. D’autre part, il ne faut pas 
oublier que, à cause de la rétrodiffusion, le spectre 8 
est fortement modifié, au fur et à mesure de la 
pénétration des électrons dans le radiateur, ce qui 
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entraîne certaines réserves au sujet de l’application 
de la formule de Wyard au cas où le radiateur 
n’entoure pas complètement la source £. 
Toutefois, nos mesures expérimentales ont été 
limitées à des valeurs de 0 comprises entre 300 et 
759, pour des raisons d’encombrement de blindage 
et de précision statistique. Le spectre de rayon- 
nement de freinage émis dans tout l’angle solide 
(6 variant de 00 à 1800) sera différent du nôtre. 
Cependant, le rayonnement émis sous des angles 0 
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grands est riche en faibles énergies, mais son inten- 
sité diminue rapidement au fur et à mesure que 0 
croît, de sorte que le spectre émis dans l’angle 
solide de 47 ne sera pas très différent de celui que 
nous avons relevé [9]. 
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En conclusion, on peut dire que l’application de 
l’équation de Kirkaptrick pour £, < 1,5 mc? et de 
l’équation de Sauter pour Æ£, > 1,5 m,c?, étendue 
à un spectre continu $ permet de prédéterminer le 
spectre de rayonnement de freinage émis dans un 
radiateur mince (d’épaisseur de l’ordre de 1/20° du 
parcours maximum des $ les plus énergétiques) 
avec une précision meilleure que 5 %. 
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‘Pour des épaisseurs de radiateur plus impor- 
tantes, le spectre s’enrichit fortement en rayon- 
nement mou, alors que le rayonnement dur n’aug- 
mente que très faiblement. 


V. Possibilités d’utilisation du rayonnement de 
freinage. — Si, du point de vue théorique, il est 
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intéressant d’étudier le rayonnement de freinage 
émis par des 8 dans un radiateur, du point de vue 
pratique, il importe surtout de connaître le spectre 
de rayonnement émergeant du radiateur. 

Les figures 17, 18 et 19 donnent quelques 
spectres émergeant du radiateur. De l’examen de 
ces résultats, on peut conclure que, pour une épais- 
seur faible de radiateur (39 mg/cm?), le rayon- 
nement mou est favorisé par rapport au rayon- 
nement dur au fur et à mesure que l’angle 0 aug- 
mente (fig. 17). Une faible augmentation de l’épais- 
seur de radiateur modifie peu les formes de spectre, 
si ce n’est un léger accroissement du rayonnement 
mou pour toutes les valeurs de 0 (fig. 18). Dans le 
cas d’une épaisseur de radiateur importante 
(fig. 19), pour 0 — 459, le spectre garde la même 
forme que pour des épaisseurs plus faibles, avec 
renforcement du rayonnement mou où apparaît 
déjà un maximum ; pour 0 — 600 et 6 — 759, le 
maximum est beaucoup plus accusé et tend à se 
déplacer dans le sens des énergies croissantes. 

L'évolution de la forme des spectres en fonction 
de 8 et de l’épaisseur de radiateur se justifie par la 
variation de l’émission de rayonnement de freinage 
et de l'importance relative de son absorption dans 
le radiateur. Pour un radiateur d'épaisseur mince 
croissante, la production de rayonnement dur varie 
peu et il est faiblement absorbé ; tandis que l’émis- 
sion de rayonnement mou reste importante pour 
chaque tranche du radiateur (sensiblement propor- 
tionnelle à l’épaisseur du radiateur) mais la self- 
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absorption peut devenir prédominante pour des 
épaisseurs assez grandes et pour des énergies assez 
faibles. 

Le maximum des spectres de radiateurs épais se 
situe aux environs de 150 — 200 keV, domaine 
intermédiaire entre les X et les y, ce qui montre 
l'intérêt pratique que présente l’utilisation du 
rayonnement de freinage [10]. 


Nous tenons à remercier M. Mare L’Allemand 
pour l’aide qu’il nous a apportée lors des mesures 
expérimentales et des longs calculs nécessités par 
notre travail. 

Manuscrit reçu le 23 janvier 1961. 
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VITESSE DE PROPAGATION DES ONDES THERMOSTATIQUES DANS LES LIQUIDES 


Par A. CARRELLI et M. MARINARO, 
Institut de Physique Expérimentale à Naples. 


Résumé. — Dans ce travail on recherche les conditions qui permettent de mettre en éviderce 
un paramètre, plus proprement une vitesse de propagation trouvée par Lucas dans les liquides. 

On propose une nouvelle méthode pour sa mesure, qui permet son évaluation pour tous les 
liquides et on rapporte les résultats obtenus pour deux liquides. 

Les différences que l’on remarque entre cette vitesse et celle du son peuvent être attribuées à 


des phénomènes de dispersion. 


Abstract, — In this work we look for the conditions which allows us to establish a parameter 
and more exactly a velocity which was found by Lucas in liquids. 

We propose a new method for this measurement, which allows us to measure this velocity in 
all liquids and we give the results which we have obtained for two liquids. 

The differences which we have noted in this velocity and in that of sound may be attributed 


to dispersion phenomena. 


D’après un raisonnement proposé par Lucas [1], 
on constate que, sous l’effet de la propagation de la 
chaleur dans un liquide, il apparaît une force 
agissant sur un flotteur qui y est plongé. Des me- 
sures expérimentales, faites par Lucas et répétées 
par A. Carrelli et E. Grossetti [2], ont confirmé 
l'existence de cette force qui peut, grâce à un rai- 
sonnement extrêmement simple, être mise en rela- 
tion avec une vitesse ® caractéristique. Si, en effet, 
nous indiquons respectivement par AV, o, € le 

volume, la masse spécifique, la chaleur spécifique 
- d’un liquide, un centimètre cube de celui-ci ab- 
sorbera, pour une variation AT de température, 
une quantité de chaleur : 


AQ = pcAVAT. (1) 


Cette énergie thermique se propage dans une 
- certaine direction z dans le fluide. Or, d’après les 
principes fondamentaux de la mécanique, on aura 
une variation AG par unité de volume de la quan- 
tité de mouvement liée à ce flux d’énergie et l’on 
peut admettre que : 

AG = AQJ®D (2) 

où Ÿ est une vitesse caractéristique du liquide. 
Si cette variation de quantité de mouvement se 
produit dans un temps At une force }, donnée par : 


FE 10 
LDAF dr 


te 
® dt 


agira sur l’unité de volume du fluide. 

Selon l’équation de la propagation de la chaleur 
on peut déterminer la valeur de d7/dt quand, par 
exemple, on suspend un flotteur dans un récipient 
contenant le liquide étudié et si l’on produit instan- 
tanément un refroidissement du liquide de bas en 


haut. Le calcul montre que le flotteur est soumis à 
une force F;(t) normale au fond du flotteur de 
section S, dirigée vers le bas et dont la valeur, en 
fonction du temps, est donnée par : 


Fil) = Vee 4 63) 
où z est la distance entre le flotteur et le fond du 
récipient et k la conductibilité thermique du liquide. 

Un premier point que nous avons voulu fixer 
dans cette recherche est l'influence de l’épaisseur 
du fond du flotteur sur la force agissant sur le 
flotteur lui-même. 

Il est évident qu’une force F apparaît, laquelle 
est due à la quantité de chaleur cédée par la base 
métallique du flotteur. Puisque le fond est métal- 
lique, on peut retenir que cette cession est due a 
toute la masse de métal ; donc cette quantité de 
chaleur et conséquemment l’effet, dépendent de la 
capacité thermique de la masse métallique du 
récipient. 

Si le fond dispose d’une grande quantité de cha- 
leur à transmettre au liquide, c’est-à-dire si 
me —= 9sS est très grand, le récipient pratiquement 
ne perd pas d'énergie et donc le flux de chaleur du 
récipient au liquide ne peut donner naissance à 
aucun effet, et la force n’apparaît pas. 

Nous avons avant tout cherché à vérifier cette 
déduction, c’est-à-dire, nous avons cherché à éta- 
blir que plus pc est grand, s et S ne subissant 
aucune modification, plus la valeur de Fest petite, 
c’est-à-dire, il faut prévoir que pour s et S cons- 
tants (en faisant varier l’épaisseur s d’autres mani- 
festations apparaissent, commeonle verraplusloin), 


les produits mcF doivent être presque constants. 
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Voici les résultats obtenus pour s — 0,3 mm qui 
comme on peut voir, sont satisfaisants. 


MéÉTAUx mc GR.CALORIE F DYNESs mcF 
Cuivre 20 A0 159 2,6.10— 
Aluminium 14,1.10—2 2,0 23410 
Argent 15,0.10—2 0,9 1,3.10—1 
Or 17,1.10—? 1,6 DA AO 


Il est à noter, cependant, que nous avons trouvé 
que la force apparaît lorsque l’on emploie aussi des 
fonds non métalliques. Ainsi, avec un fond de 
celluloïde ou de verre, par exemple, pour 
s — 0,70 mm, le produit mcF acquiert la valeur 
de 77,0 10—? pour le celluloïde et 7,7 10—2 pour le 
verre, c’est-à-dire valeurs qui dans un cas sont 
notablement plus petites. 

Un second point consiste à établir si le polissage 
du fond du flotteur a une influence sur la force, 
puisque la force F provient de la surface du fond 
et cette surface, pour produire l’effet demandé, doit 
être plane. Dans ce but, nous avons utilisé comme 
fond des petits disques de verre dépoli ou gravé, et 
l’on a remarqué que, toutes les autres conditions 
étant conservées, la force, avec ce fond non dépoli, 
subit une diminution importante. 

On a enfin procédé à l’élaboration d’une nouvelle 
méthode par la détermination du paramètre ® qui 
apparaît dans la seconde relation (2). Le but de 
cette détermination consiste essentiellement à 
chercher s’il existe un lien entre ce paramètre et la 
vitesse V du son dans le liquide à l’étude. (Un 
calcul très simple démontre plus particulièrement 
que ® devrait être égale à 4/3 V.) 

Il est à remarquer que lorsque l’on se sert de la 
méthode proposée par Lucas, la mesure de la 
force F est rendue difficile par le fait que l’on aen 
même temps une variation de la densité du liquide 
à cause du refroidissement subi par celui-ci ; cette 
variation de densité entraîne une variation de la 
poussée d’Archimède et, par conséquent, l’appa- 
rition d’une force FA opposée à Fr, et dont la 
valeur peut être beaucoup plus grande, rendant 
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impossible la mesure de F# dans de nombreux 
liquides. 

On peut surmonter cet inconvénient, en consi- 
dérant que lorsque l’épaisseur du flotteur est très 
grande (de l’ordre de 2 à 3 mm) l’effet Lucas, pour 
les raisons indiquées ci-dessus, disparaît. 

Prenons donc, comme sur la figure 1, deux 
flotteurs constitués par deux cylindres de verre 
de 5,5 em de hauteur et de 4 em de diamètre. 

La base de A est formée par un disque de laiton 
de 0,03 cm d’épaisseur. Celle de B par un disque 
de 0,3 cm d’épaisseur. Tous deux sont placés dans 
deux récipients de structure analogue et à la même 
distance du fond ; lorsqu'on produit un refroidis- 
sement, on a sur chacun d’eux une variation de 
poussée qui a pratiquement la mêmé valeur, tandis 
que, selon ce qui a été dit plus haut, la force de 
Lucas apparaît uniquement sur celui des deux 
flotteurs dont l’épaisseur est mince, cette force étant 
seule responsable du déséquilibre qui s’ensuit et 
pouvant être mesurée. Les mesures ont été effec- 
tuées à l’aide d’une balance très sensible, ayant une 
période d’oscillations d’environ 19 s. L’axe de rota- 
tion portait un petit miroir et sur un tambour on 
pouvait obtenir une trace qui donnait la mesure de 
la rotation. 

La force qui agit sur le fond du flotteur est 
donnée par la relation (3). C’est une fonction carac- 
téristique du temps. L’équation du mouvement de 
la balance est donc la suivante : 


L'intégrale de (4), compte tenu de l’expression 
théorique de F,(1) a été calculée numériquement 


et l’on a ainsi, dans le cas de l’eau, la courbe théo- 


€ secondes 


180 


EG, 


+ ota = Fu(i. (:) 


fi 
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rique représentée sur la figure 2. Il faut noter que 


la variation de « en fonction de #, obtenue par l’inté- 


gration de (4), c’est-à-dire en tenant compte de la 
période 7 — 19 s, et de l’amortissement de la 
balance utilisée, correspond pour des valeurs éle- 
vées de £ à la variation de la fonction : 


ft) = Vreli = e 
V£ 


—2#pc/4kt 


(voir fig. 2). Les enregistrements expérimentaux 
qui donnent la valeur de « en fonction du temps peu- 


._ vent être comparés avec /(t) plutôt qu'avec l’inté- 


 grale de l’équation (3) pour des valeurs élevées du 
temps (pratiquement pour & — 100 dans notre cas). 
| Ceci revient à négliger, pour & > 100, l'influence du 


mouvement oscillatoire de la balance, en Suppo- 
sant que le déplacement du bras est directement 
proportionnel à la force appliquée, c’est-à-dire, indi- 
quant, a étant une constante de proportionnalité 
entre force appliquée F£ et déviation : 


D = Veelk + 


En utilisant divers liquides et enregistrant les 
déviations « en fonction de t on a trouvé des varia- 
tions d’un liquide à un autre outre la variation due 
à la diversité de k, e, c pour les divers liquides ; et 
ceci démontre encore que cette vitesse ® existe 
effectivement. 

Les enregistrements expérimentaux pour deux 
liquides sont reportés sur la figure 3a et b en les 
comparant avec les graphiques de f(t) (sur la 
figure 4) on obtient la valeur de ®. Le tableau 
donne les résultats obtenus : 


& = aFx = a — Ra ee 


X 
LiQuipes pee pen 
sec 3. sec 
Alcool méthylique 2:8-10P 1,49.105 
Alcool amylique 1588 102 MANU 
Eau 21,0.105 2,0 .105 


Comme on peut voir, conformément à ce quia 
été obtenu précédemment, les valeurs de ® peuvent 
être considérées comme distinctes de celles de 4/3 V, 
Des recherches faites avec cette méthode et avec 
des installations instrumentales plus précises pour- 
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Alcool méthylique 
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ront confirmer d’autant mieux cette différence 

entre ® et V. 

La propagation dans un liquide d’un flux de cha- 
leur, est une propagation d’énergie et de quantité 
de mouvement ; la vitesse avec laquelle la quan- 
tité de mouvement se propage peut évidemment 
se comparer à la vitesse des ondes longitudinales 
ayant une longueur d’onde très petite. 

Les valeurs de cette vitesse sont connues pour 
des longueurs très grandes en comparaison avec 
celles des groupes moléculaires, et dans cette 
zone, V est indépendante de à, c’est-à-dire on n’y 
rencontre pas dispersion. Les valeurs de ® obte- 
nues avec cette méthode pour certains liquides, 
différentes des valeurs de V, feraient donc penser à 
une dispersion quand, comme on a dit, on utilise 
des ondes de À très petites. 


Manuscrit reçu le 10 janvier 1961. 
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ÉTUDE CINÉMATIQUE 
DE L’'ÉLECTRO-PRODUCTION DE MÉSONS x 


Par M. Paul KESSLER, 
Laboratoire de Physique Atomique du Collège de France. 


L'étude qui suit s’applique en fait non seulement à la 
production de mésons x par les électrons de haute 
énergie, mais également et de manière identique à leur 
production par d’autres particules, pourvu que celles-ci 
soient suffisamment relativistes, c’est-à-dire telles que 
leurs éne gies et impulsions puissent être confondues 
(en prenant c — 1). Notamment, elle concerne les expé- 
riences actuellement en préparation sur la diffusion 
élastique et inélastique de mésons y à très haute énergie 
(de l’ordre du GeV au moins). 


1. Raprel des données ceinématiques de la diffusion 
élastique (fig. 1). — Soit E et E’ respectivement l’im- 
pulsion (confondue en valeur absolue avec l'énergie) 


FE" 
0 
. 
M 
El 
F1G. 1. 


initiale et finale de la particule, 0 l’angle de diffusion, 
M la masse du noyau-cible, g l’impulsion transférée au 
noyau. £, M et 0 sont donnés. 
Les équations de conservation sont : 
OS CN | (a) 
E—E + M=Va + M, (b) 


soit encore, en élevant chaque membre des deux équa- 
tions au carré : 


(HE) L4RBECSin2 0/2 "7? 
(E— EL 9M(É —E'") L M°=g + M? 


En soustrayant membre à membre ces deux équa- 
tions, on obtient. : 


LEE’ sin? 0/2 — 2M(E — EF’) = 0, 
et en posant : 


y = [1 + 2(EJM) sin? 0/2], 


on à : 
E' = YE 
et : 
= 2MET =) EU ENT 


F 
{+ 
Ch 


5 
1 


2. Détermination du seuil de production — Soit pu 


la masse du méson 7. Les équations de conservation au 
seuil de production sont : 


E—E = q (a) 
E—E + M=Ve + M +u, (b) 

soit encore : | 
(E— EP +2 TME = q (a’) 


(E:— E°)° + 2M(E = E)= qe UV EEE 
(b°) 


En soustrayant ces deux équations membre à 
membre, on obtient : 
(2/y) ME’ —2ME = — Ju Va? + M? —u2 © 


Posons € — (u/M) et notons que e € 1. L’équation 


ci-dessus peut s’écrire : 


E’ — ÿE — ye Vqg? + M?— (yJ2) e? M. (c’) 


Le second membre de cette équation donne la dis- 


tance entre le seuil de production et le pic élastique. « 


Si l’on se borne à l’approximation d’ordre un en €, on 
vérifie facilement qu'il suffit de négliger le second 


terme du second membre, et de remplacer dans le il 


premier terme g? par qà (la valeur de g? dans la diffu- 
sion élastique). On a ainsi pour le seuil de production : 


EYE — yeVoù + M?. 


Si l’on veut pousser l’approximation jusqu’au second 
ordre, il faut se servir du système d'équations formé 
par (a’) et (c’), de la manière suivante : Dans (c'), on 
transfère au premier membre le terme en e?, et on 
élève ensuite les deux membres au carré. On obtient 
ainsi (en négligeant un terme en ef) : 


(E° 


VE)? + ye? M(E' — VERTE e?(g? + M). 


On'‘remplace ensuite g? par son expression en fonc- 
tion de Æ" donnée par (a'). On obtient ainsi une équa- 
tion du second degré en Æ”, dont la solution est (en 
négligeant les termes à partir de eÿ) : 


E'=YE— yeVqù + M? 


Y) (VE — Mi]. 


3. Corrélation entre le recul du noyau et l’impulsion 
du méson x créé ( fig. 2). — Ici, nous supposons donnés, 
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en plus de Æ, M et 0, l'énergie finale £' de la particule 
diffusée (£’ inférieur au seuil de production). Dans ces 
conditions, on peut poser : 


E—E|- a; E—E +M=—b. 


E’ 


Tee; 


Soit k l'impulsion du méson x créé, et + l’angle de k 
avec q. Notons que l’ensemble du processus (création 
du méson et recul du noyau) présente une symétrie 
cylindrique autour de l’axe portant le vecteur E — E. 
Les équations de conservation sont : 


a? — q? + 2kq cos o + k? (a) 
b = Va? + M? +VRk? + pe. (b) 


L’équation (b) permet de déterminer k# pour toute 
valeur de gq, l’équation (a) détermine ensuite ©. La 
courbe Æ(g) est limitée du fait des inégalités : 


g+k>a; g — k| < 4; 4 < max ; k < Kmax, 
où 


Imax — V(b —u)}— M?; Kmex — V{b — M}? — p?. 


4. Application numérique. — Nous considérons des 
électrons ou des mésons y d’une impulsion de 2 GeV/c, 
diffusés sous un angle de 309 sur des protons. Les 


données initiales sont donc les suivantes (toutes les 
énergies étant exprimées en GeV) : 


E—2; M—0,9%; 6— 300 


1) Diffusion élastique. On trouve : 
PR 1:55 


— 1,02. 
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2) Seuil de production. Avec p — 0,14, on trouve : 
E' = 1,88. a 


3) Corrélation entre g et k. Nous nous plaçons par 
exemple à une énergie finale £' — 1. On trouve pour 
les différentes constantes : 


a —= 1,24 : b —= 1,94 max — 1,94 : Fax —= 0,99, 


La courbe Æ(q) a l'allure indiquée sur la figure 3, où 
les droites correspondant aux différentes limites sont 
indiquées en pointillé. 

Le présent travail a été réalisé sous les auspices du 
Centre National de la Recherche Scientifique, et a 
bénéficié également de l’aide du Commissariat à 
l'Energie Atomique. 


Lettre reçue le 2 mars 1961. 


ÉMISSION AU COURS DU TEMPS 
DE QUELQUES RAIES D'IONS 
PRODUITS 
PAR UNE ÉTINCELLE CONDENSÉE SOUS VIDE 


par M. Jean-Louis Scxwog 
et Mile Germaine BALLOFFET, 


Laboratoire des Hautes Pressions (Bellevue) 


L'étude des variations en fonction du temps de l’in- 
tensité des raies émises par des étincelles est généra- 
lement limitée aux raies d’atomes neutre [1] ou d'ions 
de degré d’ionisation 1 ou 2 [2] [3] produites par des 
étincelles ordinaires ; on observe, dans l’ensemble, 
que le maximum d’émission d’une raie donnée se pro- 
duit d’autant plus tard que l’énergie requise pour l’ob- 
tention de cette raie est plus faible [2] [3]. Il nous 
a semblé intéressant de compléter les données exis- 
tantes en étendant ces mesures au spectre des étin- 
celles sous vide, dont on sait [4] qu’il est émis par des 
ions de degrés d’ionisation divers couvrant un domaine 
d'énergie compris entre 20 et 300 eV avec les puis- 
sances que nous utilisons normalement (conditions 
maxima : 8uF-25 kV). Rappelons ici que l’on calcule 
l’énergie nécessaire à l’obtention d’une raie en ajoutant 
les énergies d’ionisation intervenant jusqu’au degré 
observé et l’énergie d’excitation du niveau considéré. 

Les mesures (*) ont été faites à l’aide d’un spec- 
trographe à vide pour l’ultraviolet lointain [5], équipé 
d’un dispositif de réception photoélectrique compor- 
tant une série de fentes de sortie montées sur une pla- 
tine matérialisant le cercle de Rowland, et préalable- 
ment ajustées sur les raies étudiées. Une commande ex- 
térieure de déplacement permet de régler en cours de 
mesure la position des fentes à mieux que 1/20 mm. 
Les faisceaux lumineux sélectionnés par ces fentes sont 
réfléchis par des miroirs et dirigés sur des récepteurs 
photoélectriques commerciaux associés à des couches 


(*) Ce travail a été effectué en liaison avec le Service 
de Fusion Contrôlée du C. E. N. de Fontenay, auquel nous 
sommes redevables, en particulier, de la réalisation de 
l’appareillage électronique et haute tersion ; nous sommes 
heureux de remercier M. Hubert de l’assistance dont nous 
avons bénéficié. 
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fluorescentes (salicylate de Na). Chaque photomul- 
tiplicateur est relié à un oscillographe ; la charge dans 
le collecteur est réduite à 1 kQ ou 2 kQ afin de pas 
déformer le signal. 

Les étincelles étaient produites en utilisant le dispo- 
sitif comportant une électrode d’amorçage [6], l’élec- 
trode émettrice (anode) étant en aluminium et conte- 
nant, le cas échéant une poudre de chlorure d’ammo- 
nium (Cl, N, H), ou un nitrate (N et O), ou de l’alu- 
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mine (O seul). La capacité du circuit de décharge pou- 
vait prendre les valeurs 2-4-6-8 LF correspondant 
respectivement à des décharges de période : 10,5-12,5- 
14 et 15 ps, c’est-à-dire des selfs résiduelles de 1,4-0,8 
et 0,7 uH. 

Enfin l’oscillogramme du courant dans la décharge 
est obtenu en utilisant comme signal la tension pro-. 
duite aux bornes d’un shunt convenable par le cou- 
rant induit dans une ceinture de Rogowski. 


TABLEAU I 
ÉLÉMENT À TRANSITION ÉNERGIE 
Hydrogène 1248,07 00) 4s 2S,y2 — 2p ?P 40,15 
A 1 862,782 3e 29 Jp Ps 31,34 
Ne 1 238,800 DST De 178,28 
Ov 1 031,912 28 Sy» —2P Pal 305,69 
O+ 1 371,287 2s2p 1P9 —9p?1D, 208,96 
ER 813,002 bob 25 3h 3h, 327 


Le tableau I donne les différentes raies ayant fait 
l’objet de cette étude, avec leur longueur d’onde, la 
transition électronique et l’énergie correspondante. 
On a représenté (fig. 1) quelques oscillogrammes ca- 


Alrrr. 1 862 À avec LUF Ovr. 1 031 À avec 4uF 


Alrrr. 1 862 À avec 8 uF 


Clyrr. 813 À avec 4UF 
PrGad 


ractéristiques de ces raies en indiquant la valeur de la 
capacité utilisée, le balayage est toujours le même : 
2 ps par division. La courbe inférieure représente l’évo- 


lution du courant au cours de la décharge, la courbe 
supérieure l’émission de la raie indiquée. Le phéno- 
mène le plus marquant et le plus général pour les raies 
étudiées qui sont des raies de résonance d’ions alca- 
linoïdes [7] correspondant à des énergies extrêmement 
différentes (tableau I) est l’existence d’un maximum « 
d'émission au même instant, soit 3 ou 4 us après le 
début de la décharge. La reproductibilité de ce maxi- 
mum en position dans le temps est de l’ordre de 1 ps. 
Comme le premier maximum du courant ne se déplace 
que de 2,6 à 3,75 us environ pour les valeurs extrêmes 
de capacité que nous avons utilisées, il est difficile de 
préciser si le maximum d’émission a lieu exactement en 
même temps que le maximum du courant. Il semble 
cependant que le premier se déplace très peu dans le 
temps, lorsqu'on augmente la capacité. Par ailleurs, la 
durée de ce premier pic d'émission est d’autant plus 
longue que la capacité est plus grande. 

D'autre part, on peut noter l’existence d’un ou deux 
maxima secondaires souvent très importants, mais peu 
reproductibles en intensité ; cependant pour toutes 
les raies étudiées ces maxima se produisent toujours 
aux mêmes instants, soit : un peu après les second et 
troisième maxima du courant respectivement. Ces 
maxima sont très nets dans le cas de l’émission par des 
poudres (raies de Ny, Clyrr, Or et Lyman « de l’hydro- 


RTC: 


Ovr. 1 031 À {courbes du haut) 
Ov. 1 371 À (courbes du bas) 


À gauche, avec 2UF ; à droite, avec 8uF. 


N° 6 


gène). Lorsquon augmente la capacité de la décharge 

(à partir de 4 uF). A l'inverse de ce qu’on observe dans 
le cas du premier maximum d'émission, les maxima 
suivants correspondent à des instants qui dépendent 
de la période du courant. 

L'étude comparative d’une raie d’un ion non alca- 
linoïde (de degré d’ionisation plus faible d’une unité) 
a été faite avec la raie de Oÿ à 1 371 À en observant 
simultanément les oscillogrammes donnés par Oy 
et Oyr (fig. 2). Le maximum d'émission de O, se produit 
de plus en plus tard lorsqu'on augmente la capacité ; 
au bout de trois microsecondes avec 2uF et au bout 
de 5us avec 8 uF, c’est-à-dire toujours un peu après 
l'instant correspondant au premier maximum du cou- 
rant. Remarquons que, le maximum de Oy, ayant lieu 
toujours au bout de 3 ou 4 us-quelque soit la capacité 
utilisée, le maximum d’émission de Oy se produit soit 
avant, soit après celui de ©, selon la valeur de la 
capacité. 

On peut penser que, à partir du moment où l’on 
dispose, au cours de la montée de courant de décharge, 
d’une énergie instantanée suffisante, soit Æ:, pour 
ioniser les atomes évaporés de l’anode jusqu’au der- 
nier électron périphérique inclus, l’émission des raies 


de résonance des ions alcalinoïdes se produira par re- 


combinaison radiative des électrons avec les ions 
formés, cette émission ayant lieu au bout d’un temps 
qui ne dépend pratiquement pas de l’énergie de la 
décharge. L’autre processus d’émission, habituel, 
dans les décharges de type courant consiste en ionisa- 
tions successives puis excitation du niveau intéressé 
(sans dépasser ce niveau). Ici les deux processus vont 
probablement contribuer à l’émission, le premier pré- 
dominant peu à peu lorsqu'on augmente l’énergie 
électrique instantanée. Lorsque cette énergie devient 
suffisamment supérieure à Æ, pour permettre l’exci- 
tation des raies de l’ion de degré i + Î ayant perdu ses 
électrons périphériques, voire pour ioniser une fois de 
plus cet ion, l’émission de l’ion alcalinoïde (de degré i) 
se produira essentiellement par recombinaison ra- 
diative, d’où un maximum d’émission décalé vers les 
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temps plus longs par rapport au cas d’une excitation 
directe par impacts successifs. Ceci semble se produire 
avec l’aluminium qui émet assez facilement les raies 
de l’ion Al, et même Al, ; effectivement l’apparition 
du maximum d’émission des raies de Al, se produit 
un peu plus tard (à 5 4s du début de la décharge, donc 
après le maximum de courant) lorsqu'on augmente 
la capacité (à 8 4F). De même pour une valeur élevée 
de l’énergie instantanée, l'émission des ions de degré 
i — 1 (0y) sera retardée par rapport à celle des ions 
alcalinoïdes (0,1), étant donné le plus grand nombre 
de recombinaisons nécessaires avant d’atteindre le 
niveau de base de cet ion (fig. 2). 

A partie d’une valeur de la capacité de 4 UF, l'énergie 
de la seconde demi-oscillation de courant et générale- 
ment de la troisième est encore suffisante pour renou- 
veler le processus de ionisation jusqu’au dernier élec- 
tron périphérique, mais il faut maintenant l’intensité 
maximum du courant pour le produire, et par suite 
le pic d'émission a toujours lieu un peu après le maxi- 
mum de courant. 

En résumé, selon la valeur de l’énergie instantanée 
on aura donc soit une émission des ions de degré i — 1 
précédant celle des ions de degré 1, soit l’inverse, et 
dans ce cas l’émission des ions de degré : — 1 se produit 
lors du refroidissement du plasma. 


Lettre reçue le 3 mai 1961. 
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MESURE BALISTIQUE 
DES PETITS FLUX D'INDUCTION. 
DISPOSITIF MULTIPLICATEUR 
DE LA SENSIBILITÉ 


Par J. CAYREL, 


Laboratoire de Physique Générale 
de la Faculté des Sciences de Bordeaux, 


. Considérons un galvanomètre balistique ayant une 
résistance propre g — 340 () et une résistance critique 


PR, = 81 340 Q. Comparons les sensibilités au flux de 
ce balistique, quand il fonctionne au régime critique : 
R=g+rex = Re, et quand il fonctionne en flux- 
mètre : À = g, en supposant négligeable devant g la 
résistance de la bobine connectée au galvanomètre, 
dans laquelle on produit la variation de flux A, à 
mesurer. 
L’équation du mouvement du cadre étant : 


2 
KO" + (r, +5) DAC D 


Soient 0, et 0 les maxima des déviations angulaires 
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du cadre pour les deux régimes considérés. La compa- 
raison des sensibilités au flux 484 — 0. /A®,, et 
Sr — 0,/AP.. dans les deux cas donne : 


: G avec eau en VEN 
2 1—(F,[2VKC) 


SUITE 


ATEN 
Pour le galvanomètre (1) étudié on avait : 
G == LS et St [Se = 1585: 


Ainsi quand la résistance du circuit du balistique 
passe de 81 340 Q à 340 Q la sensibilité au flux est 
seulement multipliée par le facteur 1,83, alors que la 
quantité d'électricité traversant le galvanomètre est 
multipliée par R]g, soit par 239. La disparité entre 
l’augmentation de la déviation et celle de la quantité 
d'électricité s'explique du fait que la division par 240 
de la résistance du circuit du galvanomètre multiplie 
dans les mêmes proportions, et la quantité d’électri- 
cité induite génératrice de l’impulsion donnée au cadre, 
et l'intensité des courants induits par l’aimant du gal- 
vanomètre dans son propre cadre au cours de la rota- 
tion de ce dernier. Ces courants ont pour effet, confor- 
mément à la loi de Lenz, de freiner énergiquement le 
mouvement du cadre. On se rend compte de l’impor- 
tance de ce freinage dans le fonctionnement en flux- 
mètre par la lenteur extrême du retour au zéro du spot 
qui exige plus d’une demi-heure. 

Si la quantité d'électricité induite développée dans 
le circuit du balistique, quand il fonctionne en flux- 
mètre (À — g), traversait le cadre du même appareil 
fonctionnant en circuit ouvert, c’est-à-dire affranchi 
de tout freinage électrique, la déviation 6 du régime 
critique serait multipliée par 239 e, soit par 650. 

Des considérations qui précèdent, il résulte que si 
l’on pouvait produire très rapidement la quantité d’élec- 
tricité induite avant que le cadre du balistique ait bougé 
dans un circuit induit, dont la résistance se réduirait 
à la résistance g du galvanomètre, et immédiatement 
après, couper le circuit du galvanomètre (AR — co), 
la déviation du balistique pour une même variation 
de flux extérieur A®.. deviendrait 650 fois plus grande, 
et non plus 1,83 fois plus grande, qu’au régime cri- 
boue (R=TR)). 

Tel est le principe de l’appareil que nous avons 
imaginé. 

Il se compose de deux rupteurs K,, K, actionnés par 


(1) Pour beaucoup de galvanomètres on a sensiblement 
Adour Sorel MES 


les valeurs théorique et expérimentale de F tient à une 


deux roues à cames en Lucoflex montées sur le mê 
axe et dont la rotation est obtenue par la tension d° 
ressort agissant sur une manivelle solidaire de l 
Une gâchette déclenche la rotation des roues à cam: 
Le premier rupteur K, est intercalé dans le cire 
inducteur. Lorsqu'il fonctionne à l’époque f, le coure 
inducteur passe de Z à O0 dans un temps très cou 
Le second rupteur K, est inséré dans le circuit ind 
composé de la bobine induite et du galvanomèt 
Le rupteur K, doit fonctionner postérieurement | 
rupteur K, à l’époque {, + +. Le temps + doit être as. 
long pour que la f. e. m. induite dans le circuit g 
vanométre puisse y produire intégralement la quantitém 
d'électricité A®,, /g. RE 
Par ailleurs, il doit être assez court pour que le cadre 
du galvanomètre n’ait pratiquement pas bougé avant 
la rupture du circuit induit, qui libère le cadre du ba: 
listique des courants de freinage. Rien ne prouve que 
les deux conditions sus-énoncées soient pleinement 
compatibles. 5 
Le réglage de + s'obtient en décalant plus ou moins 
la roue à came du rupteur K, par rapport à celle du 
rupteur K;,, et aussi en réglant la tension du ressort 
qui actionne les roues. La mesure de + à l’oscilloscopex | 
nous a donné un temps de l’ordre de 107 s. 44 
La valeur expérimentale du facteur F, de multiplis 
cation de la sensibilité par rapport à celle du régime 
critique atteint 382, soit 59 % de la valeur théorique 2 
Rce]g = 650, $ 
Autrement dit, l’élongation maxima du spot, que 
nous obtenons avec notre dispositif, pour un courant 
inducteur de 2 mA, est de 130 mm sur une échelle de 
lecture placée à 1 mètre, alors que, sans le dispositifs 
on obtient une élongation maxima de 0,34 mm. Bier. 
entendu, cette dernière valeur est déduite d’une mesure 
avec courant inducteur de 100 mA qui donre uné 
élongation de 17 mm, d’où avec 2 mA au lieu de 100 mA 
la déviation 0,34 mm. 
Le problème qui nous reste à résoudre est de savoiïf 
si le facteur de multiplication de 380 peut être rappro= 
ché de la valeur théorique 650 par une réduction dem 
l'intervalle de temps + séparant le fonctionnement des 
deux rupteurs. Ou si, au contraire, la différence entre 


incompatibilité entre les conditions théoriques fixant. 
pour le galvanomètre considéré, les valeurs minimum 
et maximum de +. Même dans cette dernière hypothèse | 
une multiplication par 380 de la sensibilité au flux est 
assez grande pour être susceptible d'applications. % 
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